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SUMMARY 


PROBLEM 

Develop  a  computer  program  to  provide  the  capability  for  analysis  and  design  of  any 
conceivable  electromechanical  vibrator  by  implementing  a  generalized  method  based  on 
multiport  network  interconnection  techniques. 

RESULTS 

The  computer  program  SEADUCER  has  been  developed  and  applied  successfully  to 
several  problems  involving  steady-state  analysis  and  design  of  various  longitudinal  vibrators 
for  Navy  sonar  systems.  Calculations  from  mathemaii.  al  models  of  transducer  flexing  heads 
have  been  curve-fit  and  then  used  in  the  program  when  vibrational  modes  other  than  plane 
wave  have  been  known  to  exist.  Applications  other  than  longitudinal  vibrator  analysis  have 
also  been  made  successfully  using  the  multiport  interconnection  techniques  available  in  the 
program. 

SEADUCER  has  been  written  so  that  it  can  easily  be  installed  on  most  medium-  to 
large-scale  digital  computers.  The  program  is  being  used  at  other  installations  and  is  avail¬ 
able  to  government  and  private  organizations. 

RECOMMEN  D  ATION  S 

Develop  a  more  general  technique  to  solve  for  selected  individual  port  variables. 

Also,  implement  hybrid  interconnections  to  provide  a  more  powerful  tool.  Additional 
mathematical  models  for  element  types  should  be  added  as  the  need  arises.  Applications 
to  other  fields  such  as  electronics  should  be  explored. 
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INTRODUCTION 


The  purpose  of  this  document  is  to  provide  a  working  guide  to  using  the  SEADUCER 
computer  program.  SEADUCER  is  an  acronym  for  Steadystate  Evaluation  and  Analysis  of 
TransDUCERs.  The  techniques  used  by  SEADUCER  to  combine  sub-assemblies  of  trans¬ 
ducers  are  applicable  to  math  models  for  a  wide  variety  of  transducers.  However,  in  this 
document  SEADUCER  is  applied  primarily  to  math  models  for  longitudinal  electromechan¬ 
ical  vibrators.  To  a  large  extent,  SEADUCER  represents  the  state  of  the  art  both  in  the 
math  models  and  in  their  computer  implementation  for  analysis  and  design  of  this  type  of 
transducer.  These  computerized  techniques  have  been  developed  and  implemented  by  the 
Transducer  and  Array  Systems  Division,  Code  601,  of  the  Naval  Undersea  Center. 
SEADUCER,  originally  called  MOD  8,  is  the  latest  in  a  long  series  of  transducer  analysis 
and  design  computer  programs  developed  by  NUC  Code  601  and  is  significantly  different 
in  two  respects: 

(1)  the  mathematical  method  for  combining  the  elements  of  the  transducer  is  more 
flexible,  and  (2)  the  program  implementation  can  be  easily  adapted  for  use  on  most 
medium-  to  large-scale  digital  computers. 

For  transducer  analysis,  SEADUCER  has  as  its  fundamental  building  blocks  the 
mathematical  models  of  transducer  subassemblies.  While  some  of  these  models  have 
existed  for  years,  their  first  combination  into  a  transducer  analysis  computer  program 
(MOD  1 )  was  begun  in  1961  by  Dr.  Gordon  Martin  of  Code  601 .  Since  that  time  several 
individuals  have  made  a  number  of  improvements  in  math  models  for  sub-assemblies  and 
in  their  use  in  transducer  analysis  programs. 

The  development  of  SEADUCER  began  early  in  1968  when  it  became  apparent 
that  the  existing  MOD  7  program  was  not  capable  of  handling  several  transducer  configu¬ 
rations  of  importance.  To  overcome  this  problem,  H.  Ding  and  J.  Ward  of  Code  601 
described  the  method  which  forms  the  basis  for  combining  elements  in  SEADUCER. 

This  multiport  interconnection  technique  is  used  primarily  for  transducer  analysis  and 
design;  however,  it  can  be  used  within  the  SEADUCER  program  for  solution  of  other 
problems  in  electronics,  mechanics,  thermal  flow,  etc. 

With  the  generalized  multiport  concepts  available,  L.  McCleary  began  the  computer 
implementation  after  appropriate  computer  input  notation  for  the  network  interconnec¬ 
tions  was  selected.  The  programming  has  proceeded  continuously  since  that  time. 

The  basic  purpose  of  this  report,  as  stated  above,  is  to  provide  the  potential  user 
with  enough  information,  including  examples,  to  enable  him  to  prepare  typical  transducer 
application  computer  runs.  The  mathematical  equations  used  as  a  basis  for  the  programs 
are  given;  but,  in  general,  no  derivation  or  detailed  discussion  of  them  is  included,  par¬ 
ticularly  if  either  is  readily  available  in  other  literature. 


The  report  consists  of  three  major  sections  in  addition  to  the  Introduction.  The 
first.  Definitions  and  Conventions,  present’s  the  terms,  definitions,  and  conventions  used 
throughout  the  report.  The  second,  Theoretical  Description  of  Multiport  Interconnection 
Technique,  contains  a  mathematical  description  of  the  multiport  interconnection  tech¬ 
nique.  The  third.  General  Description  of  Computer  Program,  presents  the  basic  features 
of  the  SEADUCER  Program  for  both  programmers  and  non-programmers.  There  are  six 
appendices,  A  through  F.  Appendix  A  is  programmer  oriented  in  that  it  provides  a 
description  of  each  subroutine  in  the  program;  the  primary  common  storage  blocks  are 
also  described.  Appendix  B  defines  the  types  of  networks  implemented,  giving  the  input 
data  required,  sign  conventions,  and  mathematical  equations.  Appendices  C,  D,  and  E 
present  discussions  and  specific  examples  of  computer  runs  for  various  applications 
involving  longitudinal  electromechanical  vibrators.  Appendix  F,  which  is  an  active  band¬ 
pass  filter  analysis,  illustrates  an  application  in  a  field  other  than  transducer  analysis. 

A  complete  computer  program  listing  is  too  large  to  include  in  this  report;  however, 
copies  of  the  program  will  be  released  in  a  form  suitable  for  a  particular  installation  on 
receipt  by  NUC  Code  601  of  an  official  request. 

In  order  to  make  effective  use  of  rapidly  advancing  technology,  the  SEADUCER 
program  and  documentation  must  be  frequently  updated.  In  order  to  accomplish  this 
easily  and  efficiently,  Appendix  A  had  been  completely  computerized  using  a  type¬ 
setting  program  so  that  future  changes,  deletions,  and  insertions  can  be  made  quickly 
and  inexpensively.  Such  information  will  be  supplied  either  as  updates  to  this  report 
or  as  follow-on  reports  when  appropriate. 


DEFINITIONS  AND  CONVENTIONS 

The  following  definitions  and  explanatory  notes  are  stated  to  aid  clarity.  Figure  1 
illustrates  many  of  the  definitions  given  here.  “Voltage”  and  “current”  are  used  here  but 
the  same  principles  can  be  applied  to  a  variety  of  problems,  e.g.,  mechanical,  in  which 
case  the  preceding  terms  would  be  replaced  by  “force  and  velocity  respectively. 

A  port  is  defined  as  a  terminal-pair  wherein  the  current  at  one  terminal  must  at 
every  instant  be  equal  and  opposite  to  the  current  at  the  other  terminal  of  the  pair.  An 
accessible  port  is  a  terminal-pair  of  a  connected  array  of  networks,  whereas  a  network 
port  is  a  terminal-pair  of  a  single  network. 

A  multiport  network  is  defined  as  a  set  of  linear  equations  relating  current  and 
voltage  (or  other  analogous  variables)  at  each  of  the  one  or  more  ports. 

A  connection  must  join  together  two  or  more  network  ports  with  a  single  acces¬ 
sible  port.  These  ports  are  said  to  be  involved  in  the  connection. 

A  series  connection  is  a  connection  wherein  all  involved  ports  have  a  common  cur¬ 
rent;  i.e.,  a  common-current  loop  is  formed  by  the  connection. 

A  parallel  connection  is  a  connection  wherein  all  involved  ports  have  a  common 
voltage;  i.e.,  a  common-voltage  node  pair  is  formed  by  the  connection. 

A  connection  involving  no  more  than  two  network  ports  can  be  treated  either  as 
a  series  or  parallel  connection,  depending  on  the  type  of  accessible  port  involved. 

A  series  port  is  a  port  involved  in  a  series  connection. 

A  parallel  port  is  a  port  involved  in  a  parallel  connection. 


An  interconnection  is  defined  as  a  cpnnected  array  of  a  number  of  multiport  net¬ 
works  containing  at  least  one  connection.  An  interconnection  has  one  or  more  accessible 
ports  and  is  itself  a  multiport  network  which  relates  current  and  voltage  (or  other  analo¬ 
gous  variables)  at  each  of  the  one  or  more  accessible  ports. 

A  series  interconnection  is  composed  exclusively  of  series  connections  and  is  called 

Type  S. 

A  parallel  interconnection  is  composed  exclusively  of  parallel  connections  and  is 
called  Type  P. 

A  hybrid  or  mixed  interconnection  may  contain  both  series  and  parallel  connec¬ 
tions. 

Compaction  is  the  process  of  removing  undesired  ports  having  trivial  boundary  con¬ 
ditions  from  an  interconnection,  leaving  only  those  accessible  ports  desired  in  the  final 
relationship.  This  is  accomplished  by  placing  the  proper  short-  or  open-circuit  boundary 
condition  on  each  undesired  port.  Decompaction ,  on  the  other  hand,  is  the  process  of 
determining  the  voltages  and  currents  at  each  port  of  every  multiport  network  in  the 
interconnection,  subject  to  the  necessary  and  sufficient  boundary  conditions  applied  to  the 
accessible  ports  of  the  interconnection. 

Several  different  matrix  forms  are  used  in  describing  multiport  networks.  The 
Z-form  relates  voltage  to  current,  i.e., 

[E]  =  [ZJ  [I] 

where  [E]  is  the  voltage  vector,  [Z]  is  impedance  matrix,  and  [11  is  the  current  vector. 

The  Y-form  relates  current  to  voltage,  or 

[I]  -  [ Y]  [E]  , 

where  [I]  and  [E]  are  vectors  as  before,  and  [Y]  is  the  admittance  matrix.  Another  useful 
form  which  is  defined  for  even-order  networks  only  is  the  A-form.  which  maps  output 
port  variables  into  the  input  port  variables,  i.e.,  for  a  two-port  network. 


The  naming  convention  adopted  for  network  ports  is 
±  ±  i  X  n  Yk 


where 


i  is  the  port  number  (arbitrarily  numbered  from  1 . n  for  each  network). 

X  and  Y  are  any  pair  of  letters  used  to  uniquely  name  the  network  or  intercon¬ 
nection, 

n  is  the  number  of  ports  possessed  by  network  XnY,  and 

k  is  the  piece  number  (used  only  for  certain  types  of  networks  which  allow  auto¬ 
matic  cascading  of  pieces  with  the  same  name  -  normally  2%  2  matrices  in  the  A- 
form). 

Note:  A  section  is  either  an  input  network  which  has  no  piece  numbers  or 
a  network  which  is  composed  of  cascaded  pieces. 


Series  Connection  Parallel  Connection 


Series  Interconnection  Parallel  Interconnection 


figure  1.  Illustration  of  Definitions.  In  many  cases  an  accessible  port  is  introduced  artificially  in  order 
to  satisfy  the  definition  of  a  connection  and  then  is  later  terminated  with  a  trivial  boundary  condition 
(opened  or  shorted). 
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The  signs  preceding  the  name  are  used  to  define  the  sign  convention  prescribed 
by  the  user  and  represent  the  sign  of  the  dependent  variable  (or  force,  for  series)  and 
independent  variable  (or  velocity,  for  series),  respectively.  For  examply,  consider 
the  network,  A2B.  shown  in  figure  2.  The  first  accessible  port  is  1 A2B,  and  the  second 
is  2A2B.  The  connection  joins  two  pieces,  A2B1  and  A2B2. 

Several  different  network  sign  conventions  may  be  used  simultaneously  with  an 
interconnection  provided  that  the  sign  convention  ascribed  to  a  given  network  is  consis¬ 
tent  with  the  given  network  equations. 


Figure  2.  Network  Illustrating  Naming  Convention. 


THEORETICAL  DESCRIPTION  OF  MULTIPORT 
INTERCONNECTION  TECHNIQUE 

The  technique  described  in  this  section  is  simply  a  transformation  of  variables  proc¬ 
ess.  The  voltages  and  currents  at  each  network  port  (old  variables)  are  transformed  into 
the  voltages  and  currents  at  each  accessible  port  (new  variables)  by  means  of  two  transfor¬ 
mation  matrices  determined  by  the  interconnection  topology  and  by  Kirchoff  s  laws.  The 
determination  of  these  two  transformation  matrices  and  subsequent  operations  are  much 
simpler  if  the  appropriate  connection  definition  (series  or  parallel)  is  used  wherever  there 
is  a  choice. 

The  following  development  assumes  a  hybrid  or  mixed  form  of  interconnection  for 
th 

generality.  Consider  the  b  interconnection  of  N  multiport  networks  where  each  network 

has  pn  ports  (n=l,2 . r,...,N).  Each  network  may  be  written  in  the  following  form  relating 

network  port  variables: 

[W[]  =  [Kfj]  [Uj]  (i,j=l,2,3,...,pr), 

where  [Wp  and  [Uj]  are  column  vectors  of  pr  components,  [Kjj]  is  a  square  matrix  of  pr 
X  pr  components,  and  and  U^  are  the  two  port  variables  at  the  a^2  port  of  the  r^ 
network. 


Proper  formulation  of  [W?']  and  [Uj]  implies  the  following  conditions  for  the  r 
network: 


th 


if  ‘a’  is  a  series  network  port,  -  E^  and  -  1^; 
if  ‘a’  is  a  parallel  network  port,  =  1^  and  =  E^, 


where  E  denotes  port  voltage  and  I  denotes  port  current. 

These  “proper  formulation’  rules  are  required  to  allow  an  automatic  application  of 
short-  and  open-circuit  boundary  conditions  to  undesired  accessible  ports  during  the 
compaction  process. 

Since  each  network  matrix  equation  is  independent  of  the  interconnection,  the 
entire  set  of  interconnected  multiport  networks  may  be  arrayed  in  a  single  matrix  equation, 
once  the  individual  network  equations  are  in  proper  form. 


(W‘) 

(Wf) 


(Kj)  (  0  )  ...  (  0  ) 


(  0  )  (Kg) 


(  0  )  (  0  ) 


(  0  ) 


<) 


(Uj1) 

(Uj2) 

'A  ■  ■ 

L - 

(i,j=l,2....,pn), 


or,  expressed  simply,  [W?1]  -  [K^]  [U-1]  (n-1 ,2,...  ,N)(i,j-l  ,2,...  ,pn).  where  [W-  ]  and 

^  N  ^  ^  a 

[Un]  are  column  vectors  of  2  pn  components  (each  a  voltage  or  current)  and  [K--]  is  a 
J  n=l  J 

N  N 

square  matrix  of  2  pn  X  I  pn  components. 

n=l  11=1  tu  . 

Consider  the  relationship  of  accessible  port  variables  for  the  b  interconnection 

(pb  in  number): 


[w£]  =  tKk£]  lU£]  (k,£=l,2,...,pb). 


where  [W^]  and  [Ug  ]  are  column  vectors  of  pb  components,  each  a  voltage  or  current,  and 
[K^g]  is  a  pb  X  pb  square  matrix. 

From  previous  definitions  and  conventions,  the  pb  accessible  ports  include,  at  a 
minimum,  a  series  accessible  port  in  every  series  connection,  a  parallel  accessible  port  in 
every  parallel  connection,  and  an  unaffected  accessible  port  for  every  unaffected  network 
port  (which  may  be  specified  as  either  series  or  parallel  connection).  If  additional  accessible 
ports  are  introduced,  they  will  lead,  at  best,  to  formation  of  hybrid  connections  necessi¬ 
tating  further  work. 


Again,  proper  formulation  of  [Wb]  and  |  Ub]  implies  that 

if  a  is  a  series  port,  Wb  =  Eb  and  Ub  =  Ib  : 
if  a  is  a  parallel  port,  Wb  =  Ib  and  Ub  =  Eb, 


where  E  denotes  port  voltage  and  1  denotes  port  current. 

In  general  (for  the  hybrid  or  mixed  interconnection),  the  W  and  U  matrices  contain 
both  E  and  I  elements  mixed  according  to  the  occurrence  of  series  and  parallel  connections 
and  unaffected  ports  in  the  hybrid  interconnection.  If  all  connections  are  series  or  all  con¬ 
nections  are  parallel,  the  W  and  U  matrices  will  not  be  mixed. 

The  transformation  matrices  are  defined  as  follows: 


iwbi  =  [cbi  twH] 


tuH]  =  [Db] 


fU 


[Wb]  =  [C^wf1]  =  [Cb]  [Ky]  [Uj1] 


(subscripts  as  before) 


[Cb][K"][Dbl[Ub] 


which  implies 


'KkC 


[Cb]  [K|]]  [Db] 


The  f Cb ]  and  [Db]  connection  or  transformation  matrices  are  determined  by  apply¬ 
ing  Kirchoffs  laws  to  the  interconnection  topology  and  will  contain  only  integer  elements 

,  N  h  N 

equal  to  zero  or  ±1.  [Cb]  has  pb  rows  and  £_  pn  columns,  and  [ DD]  has  2  pn  rows  and 

pb  columns.  n=1  n_1 

If  there  exists  a  connection  between  two  or  more  network  ports  without  at  least  one 

involved  accessible  port  (contrary  to  definition),  a  row  or  column  of  the  [Cb]  and  [Db] 

matrices  will  contain  all  zeros  and  will  cause  an  invalid  variable  transformation. 

Consider  the  matrix  equation  for  all  accessible  ports  for  the  b^  interconnection 

following  transformation  of  variables: 


[Wb]  =  [Kbg]  [Ub]  (k,B=l,2,...,pb). 


Many  of  the  accessible  ports  present  during  transformation  of  variables  may  have  trivial 
boundary  conditions  (i.e.,  open-  or  short-circuit)  and  this  may  not  be  desired  in  the  final 
relationship.  Compaction  of  the  interconnection  is  accomplished  as  follows: 


Let 


a  -  number  of  desired  accessible  ports, 

1 3  =  number  of  undesired  accessible  ports  (to  be  shorted  if  series 
and  opened  if  parallel  in  form). 


Therefore, 


Pb  =  a  +/?• 

The  matrix  equation  after  transformation  of  variables  must  be  partitioned  (by  the  method 
described  later  in  this  section).  The  rearranged  and  partitioned  equation  may  be  written  as 
(dropping  the  interconnection-naming  superscript,  “b”) 

'[Wa]l  [" [Kac]  [Kad]l  r [Uc3l  (a,c=  1 ,2, ...  ,a) 

[wb]J  =  [[KbC]  tKbdJ  |_[UdlJ  (b,d=l. 2, .../)). 

The  undesired  ports  may  be  removed  by  setting  [W^]  =  [0] .  This  effectively  shorts 
a  series  port  (E  -»  0)  and  opens  a  parallel  port  (I  -►  0)  which  generally  represents  the  trivial 
boundary  conditions.  Therefore,  setting  [Wb]  =  [0]  and  rearranging  terms,  we  find  that 

[Ud]  =  [Kbcl[Uc], 

So 

[Wa]  =  [[Kacl  -  [Kadj  [Kbd]->  [ Kbc l]  iuc], 

or 

[Wa]  =  [Hac]  [Uc], 

where 

IHac]  =  lKac]  -  £Kad]  [Kbd]-'  f Kbc ] 

This  is  the  generalized  compacted  interconnection  equation  relating  voltages  and  currents 
(or  analogous  quantities)  at  the  a  desired  accessible  ports  of  the  N  interconnected  multi- 
port  networks. 

To  illustrate  the  decompaction  process,  consider  the  generalized  compacted  inter¬ 
connection  equation  for  the  b^  interconnection: 

[Wa]  =  [Hjc]  [uj?]  (a,c=l,2,...,a). 

The  applied  boundary  conditions  are  such  that  [Uc]  can  be  evaluated  as  a  function  of  the 
known  terminating  impedances,  voltages,  and  currents,  i.e., 

(It  is  necessary  and  sufficient  that  a  independent  boundary  conditions 
be  specified.) 
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where  W,  W/U,  and  U  are  the  necessary  boundary  conditions  at  the  accessible  ports  ot  the 
compacted  interconnection.  From  the  previous  development  of  the  compacted  intercon¬ 
nection  equation,  [U^]  may  be  computed: 


K 


be 


1  [Ub| 


So,  [Ub] 


may  be  computed: 


Also.  [Uj1]  may  be  computed: 


[Uj1]  -  [Db]  [Ub], 


(a,c  =  1 ,2....  ,a) 
(b,d=l  .2 . (3). 


(K=l,2 . p5). 


Decompaction  allows  computation  of  all  network  port  variables,  subject  to  boun¬ 
dary  conditions,  as  follows: 


[Uj1]  =  [Db] 


'[Ucbl 

.[-iKbdl_1  lKbc>  IUcfl 


(Subscripts  as  before) 


[Wj1]  =  [Kg]  [Uj1]. 


PARTITIONING  METHOD 

A  partitioning  method  for  rearrangement  of  variables  can  be  used  to  convert  any 
network  equation  into  a  more  desirable  form  than  might  occur  as  a  matter  of  course.  Con¬ 
sider  the  following: 

[W^]  =  [K^t]  [U^1  (a, 7=1 ,2,...  ,pb). 


or  simply. 


Where  certain  elements  of  [W]  are  voltages  (at  series  ports)  and  others  are  currents  (at 
parallel  ports).  Rearrangement  and  partitioning  yields 


where 


and 


*[Wsf 

’[Hj]  [H2f 

‘[usr 

.  IWpl. 

[H3]  [H4] 

_[uP]_ 

[Ws]  and  [ Ug]  exist  at 

series  ports 

[UpJ  and  [Wp]  exist  at  parallel  ports. 


i.e.,  the  series  port  variables  and  parallel  port  variables  have  been  grouped  in  each  of  the 
vectors.  Although  similar  manipulation  will  produce  various  forms  of  this  matrix  equation, 
the  relationship  in  the  Z-form. 


[wsr 

"rz, ]  [z2] 

"  [Us] " 

[Up] 

[z3]  [Z4] 

_[Wp]_ 

[E]  =  [Z]  [I] , 
is  found  by  solving 


[Ws]  =  [H,  ]  [Us]  +  [H2]  [Up] 
[Wp]  -  [H3]  [UsI  +  [H4]  [Up] 


for 


[Up]  =  [H4]-1  [Wp]  -  [H4]-1  [H3]  [Us], 


We  find  that 


[Ws]  =  [H2][H4]-,[Wp]  + 


[Hj] 


[H2]  [h4i -1  [H3] 


[Us]. 


Therefore. 


[I H j  ]  -  [H2]  [H4]-]  [H3]]  [Ht]  [ H4 ] —  1 
.  -  [H4]-1  [H3]  [H4]-1 
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The  resultant  Z-form  relates  the  voltages  at  the  network  ports  to  the  currents  at  the  network 
ports. 

In  general,  this  partitioning  method  is  used  in  converting  individual  multiport  equa¬ 
tions  into  required  form,  and  is  used  in  compaction.  Also,  if  it  is  necessary  to  short  a  parallel 
port  (H  —  0)  or  open  a  series  port  (I  -*  0),  this  method  can  be  used  prior  to  compaction  to 
arrange  the  elements  of  [W]  appropriately.  Finally,  in  setting  boundary  conditions,  this 
method  is  used  to  solve  for  [Uc]  from  the  actual  boundary  conditions  applied  to  the  desired 
accessible  ports. 

The  current  computer  program  will  not  handle  hybrid  interconnections  (i.e.,  mix¬ 
tures  of  series  and  parallel  interconnections);  therefore,  only  two  types  of  interconnections 
will  be  discussed. 

Type  S 


This  interconnection  is  composed  exclusively  of  series  connections.  Type  S  requires 
the  use  of  a  Z-form  description  of  the  individual  networks,  which  is  a  common  form  tor  net¬ 
work  description.  This  type  occurs  most  frequently  in  transducer  models.  Figure  3  shows 
an  example  of  a  series  interconnection  of  three  multiport  networks  having  two  desired 
accessible  ports  where  loop  4  involves  only  two  network  ports  and,  therefore,  may  be 
treated  as  either  a  parallel  or  series  connection.  Figure  4  shows  the  interconnection  with 
all  four  accessible  ports,  one  being  the  unaffected  network  port  and  one  involved  in  each 
of  the  three  series  connections. 


Figure  3.  Series  Interconnection. 
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Figure  4.  Series  Interconnection  With  Sign  Convention. 

The  following  equations  describe  each  of  the  networks  shown  in  the  figure: 


The  Z-form  was  chosen  since  every  port  involved  in  the  connection  is  a  series  port  and  the 
description  of  the  ports  should  therefore  use  current  as  the  independent  variable,  i.e., 

[Wf]  =  [Ej1]. 
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[Uj1]=[lj1]. 

The  entire  set  of  network  equations  can  be  written 

[W?1]  =  [K?j]  [Uj1]  (n=A,B,C). 

or 

[Ej1]  =  [Z-j]  [Ij1]  (for  series  ports). 

Thus, 
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Next,  we  write  the  relationships  between  voltages  and  currents  at  the  four  accessible 
ports  of  the  interconnection  (naming  it  “D”): 

lWk]  =  [KkC]  [Uf]  (k.C  =  1 ,2,3,4). 

Since  all  four  accessible  ports  are  series  ports, 

|wk]  =  |EkK 

fJ?i  = 


and 
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The  transformation  of  variables  is  accomplished  by  means  of  the  [CD]  and  [DD]  matrices, 
which  can  be  determined  by  inspection  of  figure  4  (for  sign  conventions,  etc.)  and  by  referral 
to  Kirchoffs  laws. 
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Finally,  using  the  expression  from  the  general  description 

[KkC]  =  [C°]  fKij]  [D°^ 

we  find  that 
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IKk(I 


0  0  0  0  0  0  0  1 

0  0  0  10  0  1  0 

0  1-10  0  10  0 

10  0  0  10  0  0 


Since  the  boundary  conditions  at  accessible  ports  3  and  4  are  trivial  short  circuits,  the  matrix 
may  be  compacted.  This  compaction  will  not  be  done  in  detail  but  will  be  sketched  as  follows 
(no  rearrangement  is  necessary  in  this  instance): 


Recall  that  if 


[Kac]  [Kad] 
[Kbc]  [Kbd] 


a. c  It  desired 

b. d  It  undesired. 
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and 


[Wb]  =0 


then 


[wa]  = 


IKac3  -  [Kad]  [Kbd]  [Kbc] 


[Uc]  -  [HD][UcJ 


Therefore, 


hd  hd 
H1 1  H12 

itD  ttD 
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which  will  yield  through  simple  matrix  multiplication  and  addition  a  final  result  in  the  form 


fcl 

n 

hd 

H11 

HD 

H12 

i 

_ j 

c2 

H21 

hd 

Qri 

_ i 

where  the  elements  of  the  H-matrix  are  known  in  terms  of  the  elements  of  the  Z-matrices  for 
the  three  networks  which  were  given  originally. 

Type  P 


This  interconnection  is  composed  exclusively  of  parallel  connections.  The  example 
is  included  here  to  illustrate  that  connections  involving  only  two  network  ports  may  be 
treated  either  as  series  or  parallel.  Consider  the  same  three  networks  used  in  the  Type  S 
example,  but  connected  in  parallel  as  shown  in  figure  5.  Node-pair  4  is  introduced  instead 
of  Loop  4  as  before,  and  is  treated  as  a  parallel  connection  here.  The  Y-form  will  be  used 
for  this  interconnection,  since  every  port  involved  is  a  parallel  port.  In  view  of  the  require¬ 
ments  of  proper  formulation  discussed  previously,  voltage  must  be  the  independent  variable. 
Thus, 

[W"]  -  [If] 

[ Uf  ]  =  [Ej1]  (for  parallel  ports). 
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Figure  5.  Parallel  Interconnection. 


If  the  Z-form  is  known  for  each  network,  the  following  conversion  from  the  Z-form  to  the 
Y-form  is  required: 

We  know  that  [E]  =  [Z]  [I], 

[I]  =  [Z] -1  [E] . 

and 

[1]  =  [Y]  [E] . 

Therefore, 

[Y]  =  [zr1. 

In  particular,  for  A,  B,  and  C, 

[  I A  j  =  [Ya]  [EA] ,  where  [YA]  =  [ZA]_1 ; 

[IB]  =  [YB]  [EB],  where  [YB]  =  [ZB]_1; 

and 

[IC]  -  [ YC]  [EC] ,  where  [YC]  =  [ZC]_1. 
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It  should  be  noted  that  the  unaffected  port  can  also  be  considered  as  a  series  port,  as 
is  done  for  Type  S  interconnection,  in  which  case  the  matrix  equation  for  network  C 
would  take  the  following  hybrid  K-form  (not  yet  implemented  in  the  computer  program 
package): 
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where  [KE]  is  determined  from  [ZE]  by  means  of  the  partitioning  method  discussed  pre¬ 
viously.  However,  this  unaffected  network  port  is  treated  as  a  parallel  port  in  this  example 
to  simplify  the  problem  and  to  allow  the  use  of  the  simpler  Y-form  notation. 

The  entire  set  of  network  equations  can  be  written  as 
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^4 

Again,  if  we  name  the  interconnection  “Network  #D’\  the  relationship  between  the  voltages 
and  currents  at  the  four  accessible  ports  is  described  by  the  following  equations: 

[Wk]  =  [KkC]  [U?]  Ck,C=  1 ,2.3,4), 
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where 


tWk]  =[Ik]’ 

[U£D]  =  [Ejjp] , 

and  [Kke1  =  [YkC]' 

Thus,  [l£]  =[Y^][eJ)]. 

The  transformation  of  variables  is  again  accomplished  by  means  of  the  [CD]  and  [D°] 
matrices, which  are  determined  from  the  diagram  and  Kirchoff’s  laws. 

[l£]  =  [CD]  [I?];  [Ej1]  =  [DD][EgD] 


20 


Finally, 
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Since  ports  3  and  4  are  trivially  open  circuited,  the  matrix  equation  can  again  be  compacted 
to  yield  this  expression: 


n 

X 

JO 

ed 

1 - - - 

2 

where  [H^]  is  now  in  the  Y-form. 

The  fact  that  only  series  and  parallel  interconnections  are  considered  here  should  not 
be  construed  by  the  reader  to  mean  that  the  theoretical  methods  presented  are  restricted  to 
these  types  of  connections.  The  theory  is  general  and  valid  for  hybrid  interconnections: 
however,  the  computer  program  will,  at  present,  handle  only  series  interconnections  or 
parallel  interconnections. 


GENERAL  DESCRIPTION  OF  COMPUTER  PROGRAM 

The  purpose  of  this  section  is  to  provide  insight  into  the  basic  features  of  the  com¬ 
puter  program  for  both  programmers  and  non-programmers.  The  description  is  organized 
into  the  following  sections: 

( 1 )  SEADUCER  Computer  Programming  Considerations, 

(2)  Network  Types, 

(3)  Interconnection  Table, 

(4)  Catalog, 

(5)  Boundary  Condition  Specification, 

(6)  Control  Program, 

(1)  Program  Conversion  from  Double  Precision  to  Single  Precision,  and 
(8)  Summary. 
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SEADUCER  COMPUTER  PROGRAMMING  CONSIDERATIONS 

Basic  design  goals  for  the  programming  of  SEADUCER  are  as  follows: 

1 .  Provide  for  computerized  analysis  and  design  of  transducers  based  on  multiport 
network  interconnection  techniques. 

2.  Produce  the  program  in  a  form  as  near  standard  as  possible  so  that  it  can  easily 
be  implemented  on  most  medium-  to  large-scale  digital  computers. 

3.  Provide  for  automatic  run-by-run  allocation  of  the  core  storage  required  for 
storing  all  input  data  as  well  as  computed  variables  and  matrices  for  any  trans¬ 
ducer  configuration. 

4.  Provide  an  extremely  straightforward  method  for  converting  the  program  from 
double  precision  to  single  precision. 

5.  Allow  for  complete  program  flexibility  by  requiring  for  each  application  the 
preparation  of  a  main  routine,  written  in  the  programming  language,  which 
explicitly  controls  the  flow  of  input,  computation,  and  output. 

6.  Provide  for  output  options  to  be  controlled  independently  of  computations 
performed.  These  options  are  for  printout  as  well  as  output  via  tape  or  cards 
to  be  processed  with  a  separate  plotting  package. 

In  order  to  implement  these  basic  goals,  the  following  measures  have  been  taken: 

1.  Even  though  the  multiport  technique  has  been  implemented  primarily  for 
transducer  applications,  a  great  deal  of  care  has  been  taken  to  allow  the 
implementation  of  other  applications. 

2.  ASA  FORTRAN  has  been  followed  as  faithfully  as  possible  throughout  the  body 
of  the  program.  The  primary  exception  is  the  use  of  the  FLD  function,  which 
must  be  provided  for  bit  handling. 

3.  By  setting  Blank  Common  to  an  appropriate  size  within  the  main  routine,  the 
user  is  able  to  set  aside  only  as  much  core  storage  as  is  required  for  a  particular 
run.  This  is  especially  advantageous  when  using  computer  systems  which  charge 
for  the  amount  of  core  used.  It  also  eliminates  the  problem  of  dimensioning 
large  areas  of  core  that  are  not  used  in  all  applications. 

4.  The  primary  mode  of  the  program  is  double  precision.  However,  a  consistent 
convention  has  been  used  throughout  the  program  making  it  simple  to  create 
a  single  precision  version.  This  allows  an  installation  to  maintain  both  single 
and  double  precision  versions,  if  necessary.  Where  20  to  23  bits  of  precision 
are  available  in  single  precision,  many  applications  will  require  double;  whiie 
for  48  bits  in  single  precision,  it  is  unlikely  that  double  would  be  required. 
Further  analysis  of  the  various  precisions  available  will  need  to  be  performed 
for  particular  applications. 


5.  Since  no  package  of  application-oriented  subroutines  can  ever  contain  every 
conceivable  option  and  variation  required  in  a  fast-moving  research  and 
development  environment,  SEADUCER  requires  that  for  each  application  a 
FORTRAN  main  routine  be  written  (referred  to  as  the  Control  Routine).  For 
the  typical  transducer  application  such  as  that  described  in  Appendix  C,  the 
Control  Routine  involves  several  subroutine  calls  to  load  the  network  configu¬ 
ration,  several  subroutine  calls  to  load  the  input  data,  and  several  calls  for  the 
desired  logic  and  output.  Only  in  cases  where  a  large  volume  of  repetitive 
data  is  required  are  actual  data  cards  read  in  for  input  parameters.  In  general, 
it  is  desirable  to  have  the  input  data  loade  d  within  the  main  routine  since 
only  in  this  way  can  values  be  directly  computed  as  they  are  loaded  into  the 
program.  This  method  also  allows  complete  flexibility  with  respect  to  the 
order  and  number  of  input  combinations  without  having  to  settle  for  options 
preset  within  the  program. 

The  success  of  this  approach  depends  upon  having  a  method  for  the  particular 
computer  system  which  allows  Tie  storage  (preferably  in  source  and  relocatable 
code)  of  the  entire  program  on  magnetic  tape  (or  disk,  drum,  etc.).  This  storage 
must  be  in  a  fosm  such  that  for  a  particular  run,  the  user  need  only  keypunch 
the  Control  Routine  and  system  control  cards  required  for  automatically 
compiling  the  Control  Routine  as  well  as  loading  and  executing  the  desired 
portion  of  the  program.  Here  is  where  sophisticated  computer  system  loaders 
are  highly  desirable,  since  they  can  be  used  to  load  only  the  sections  of 
SEADUCER  needed  for  a  particular  run.  Presently,  the  entire  program  cannot 
be  loaded  into  a  32K  word  machine  at  one  time,  although  with  proper  selec¬ 
tive  loading  techniques,  the  same  computer  can  be  used  for  many  of  the  runs 
of  interest.  When  this  limit  has  been  reached  in  a  particular  computer,  the 
only  solution  is  tc  have  a  system  for  dynamic  program  segmentation  (i.e., 
where  sections  of  logic  are  automatically  rolled  into  and  out  of  core  during 
the  course  of  the  runs).  Also  of  significant  benefit  to  the  user  is  a  computer 
operating  system  which  offers  the  capabilities  just  mentioned  as  well  as  the 
ability  to  selectively  modify  and  recompile  any  routine  previously  stored 
away,  and  to  use  the  newly  compiled  version  in  the  same  run. 

6.  In  order  to  inhibit  vast  amounts  of  unnecessary  printout,  most  SEADUCER 
subroutines  have  a  print  option  in  a  common  block  which  can  be  set.  Since 
SEADUCER  is  so  large  and  since  plotting  hardware  and  software  are  so  non¬ 
standard  among  computer  installations,  the  information  which  needs  to  be 
plotted  from  SEADUCER  is  output  on  cards  or  magnetic  tape  in  a  form  con¬ 
venient  for  loading  into  a  separate  plotting  program. 

NETWORK  TYPES  (GENERAL  DESCRIPTION) 

All  networks  input  by  the  user  must  have  a  special  type  number  associated  with 
them.  The  routine  associated  with  this  particular  type  number  is  preprogrammed  to 
compute  the  elements  of  a  specific  form  of  the  matrix  (generally  Z-form)  as  a  function  of 
specified  independent  variables  (frequency  in  all  applications  of  transducer  analysis  and 


23 


design).  The  computations  are  in  general  a  function  of  ( 1 )  the  independent  variable,  (2)  the 
sign  convention  assumed  for  the  network,  (3)  the  number  of  ports,  and  (4)  the  port  num¬ 
bering  originally  specified. 

A  library  of  the  network  types  encountered  in  transducer  analysis  presently  exists. 

It  should  be  noted  that  this  library  is  continually  expanding  as  different  networks  are  used, 
and  may  be  expanded  to  include  other  matrix  forms  or  other  physical  problems.  Specific 
types  already  implemented  are  described  in  Appendix  B. 

The  user  can  specify  types  of  networks  other  than  those  available  in  the  existing 
library  with  the  following  specifications: 

(1)  network  equations  (matrix  form;  element  equations), 

(2)  port  sign  conventions  and  numbering,  and 

(3)  definitions  of  the  input  constants  needed  to  compute  the  network 
matrix  elements. 

If  the  network  under  consideration  is  actually  composed  of  a  number  of  two-port  net¬ 
works  cascaded  together,  the  user  may  take  advantage  of  a  short-cut  technique  used  widely 
for  transducer  problems.  Each  piece  is  simply  considered  as  a  network  having  the  same 
name  as  the  composite  except  for  an  additional  number  of  the  right  side  of  the  name  which 
designates  piece  number.  The  program  automatically  cascades  the  A-form  matrices  repre¬ 
senting  each  piece  of  the  same  name  together  to  acquire  the  matrix  for  the  composite 
network,  called  a  section. 

INTERCONNECTION  TABLE 

All  interconnections  of  linear  networks  may  be  considered  to  be  single  networks 
themselves  having  similar  properties,  i.e.,  sign  conventions,  port  numbering,  boundary  con¬ 
ditions,  and  names.  The  accessible  ports  are  the  ports  of  the  interconnection,  and  network 
ports  are  the  ports  of  the  component  networks.  All  connections  involve  just  one  accessible 
port  and  at  least  one  network  port.  The  interconnection  table  is  a  description  of  the  inter¬ 
connection  for  the  computer  program  and  contains  sufficient  information  to  allow  com¬ 
plete  non-ambiguous  reconstruction  of  the  interconnection. 

Basically,  for  each  connection,  the  following  information  is  required  to  describe 
the  interconnection: 

( 1 )  connection  type  (i.e.,  series  or  parallel), 

(2)  sign  convention  at  the  connection, 

(3)  accessible  port  name,  and 

(4)  names  of  all  involved  network  ports. 

For  example,  consider  the  interconnection  shown  in  figure  6. 


1D4D 


The  information  contained  in  these  equations  can  then  be  transferred  into  an  intercon¬ 
nection  table  as  follows: 

S  ++ 1  D4D,++3C3C 
S  ++2D4D,  ++2B2B,  ++2C3C 
S  ++3D4D,  ++1  A2A,  ++1C3C 
S  ++4D4D,  -+2A2A,  ++1B2B 

where  S  denotes  series  connection,  and  signs  are  for  voltage  and  current  (or  force  and 
velocity,  etc.),  respectively. 
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The  interconnection  table  (or  port  tabic)  is  a  basic  input  into  the  computer  pro¬ 
gram.  Each  interconnection  must  be  described  via  its  own  interconnection  table.  This 
information  is  stored  and  defines  the  system  of  linear  equations  that  is  to  be  solved  for 
the  program. 

A  powerful  feature  of  Seaducer  is  its  provision  for  nesting  of  interconnection 
tables.  Therefore,  by  judiciously  choosing  interconnection  tables,  the  user  can  substruc¬ 
ture  his  problem  so  that  portions  of  the  entire  network  which  involve  changes  in  the 
physic„l  characteristics  of  components  can  be  studied  in  depth  without  recomputing  the 
entire  network.  (See  Appendices  C  and  D.) 

CATALOG 

The  program  utilizes  a  master  storage  area  to  retain  all  input  data  as  well  as  com¬ 
puted  variables  and  matrices  for  the  individual  and  combined  networks.  As  each  new 
network  is  encountered,  this  area  is  dynamically  allocated  in  the  sense  that  pointers  (sub¬ 
scripts)  are  set  and  stored  away  in  the  catalog.  The  catalog  stores  the  following: 

(1)  network  name, 

(2)  matrix  size, 

(3)  piece  number  (for  cascaded  pieces), 

(4)  total  number  of  pieces  (for  a  section  composed  of  cascaded  pieces), 

(5)  network  type, 

(6)  total  amount  of  storage  reserved  for  basic  input  values  and 
permanently  stored  computed  quantities,  and 

(7)  location  in  master  storage  area  for  network  with  above  properties. 

For  example,  the  catalog  for  the  network  shown  in  figure  6  might  be 


Name 

Matrix 

size 

Piece 

no. 

Total 

number  of 
pieces 

Network 

type 

Data 

storage 

Location  in 
master  storag 

1 

DD 

4 

0 

0 

0 

0 

1 

AA 

2 

0 

1 

0 

0 

65 

3 

AA 

2 

1 

1 

1 

7 

81 

4 

BB 

0 

0 

1 

0 

0 

104 

5 

BB 

2 

1 

1 

1 

7 

120 

6 

CC 

3 

1 

0 

1 1 

34 

143 

It  should  be  pointed  out  that  for  each  network  port,  twice  the  amount  of  storage  needed  is 
allocated  to  allow  the  user  storage  for  different  forms  (Z-form,  A-fonp,  etc.)  of  the  matrix 
other  than  that  provided  by  the  calculation  of  the  “ network  type ”  subprogram.  In  summary, 
the  catalog  is  used  to  handle  the  tedious  bookkeeping  task  involved  in  the  storage  of  the  data 
and  matrices  of  each  network  involved  in  the  interconnection. 


BOUNDARY  CONDITION  SPECIFICATION 


Boundary  conditions  must  be  specified  when  the  process  of  compaction,  or  decom¬ 
paction,  is  to  be  carried  out.  Recall  that  compaction,  as  described  under  Theoretical 
Description  of  Multiport  Interconnection  Technique,  is  the  process  of  removing  undesired 
ports  from  the  interconnection  through  application  of  trivial  boundary  conditions  (open 
circuit  or  short  circuit).  The  compaction,  as  done  automatically  by  the  program,  always 
assumes  zero  left-hand  side  in  the  equations  relating  variables  at  the  undesired  ports.  This 
effectively  shorts  a  series  port  and  opens  a  parallel  port.  Compaction  is  called  automati¬ 
cally  for  interconnections  of  the  same  name  which  have  fewer  ports  in  a  suceeding  inter¬ 
connection  table.  For  example,  if  an  interconnection  table  exists  for  A6B,  and  the  name 
A4B  appeared  in  a  later  interconnection  table,  then  A6B  would  automatically  be  reduced 
to  A4B  assuming  ports  5  and  6  to  be  the  undesired  ports. 

The  decompaction  process  is  the  means  of  determining  voltages  and  currents  at 
each  port  of  every  multiport  network  in  the  interconnection  subject  to  sufficient  boun¬ 
dary  conditions  specified  by  the  user.  The  generalized  decompaction  process  is  currently 
not  available,  although  the  force  and  velocity  at  the  edge  of  each  piece  of  a  3-3  mode 
ceramic  stack  can  be  found,  and  decompaction  can  be  done  for  all  cascaded  two-port 
networks. 

CONTROL  PROGRAM 

The  program  is  run  by  a  user-supplied  control  routine  that  contains  the  prob'em 
specification.  This  control  routine  defines  the  interconnection  tables,  defines  the  network 
type  for  each  network  name  appearing  in  an  interconnection  table,  and  specifies  the  cal¬ 
culations  and  output  to  be  performed.  The  control  program  is  a  FORTRAN  routine 
supplied  by  the  user  which  allows  added  generality  for  the  solution  of  the  user  s  specific 
problems.  Essentially,  a  library  of  subroutines  is  provided  to  perform  various  operations 
pertinent  to  the  solution  of  any  problem  which  can  be  represented  by  an  interconnec¬ 
tion  of  multiport  networks.  Examples  of  the  control  routines  used  in  the  solution  of 
several  different  types  of  problems  are  presented  in  the  appendices.  Also  included  in 
Appendix  A  are  write-ups  of  all  subroutines  available  with  the  present  program. 

The  first  step  to  be  taken  in  every  control  routine  is  to  define  the  interconnection 
topology,  or  port  table.  This  is  done  by  calling  a  subroutine,  LDPORT,  once  for  each 
accessible  port,  by  specifying  the  names  of  the  network  ports  involved  in  the  connection, 
and  by  specifying  the  connection’s  sign  conventions  and  whether  it  is  series  or  parallel. 

Next,  each  network  port  is  specified  by  calling  the  appropriate  (corresponding 
to  network  type)  routine  that  stores  in  the  master  storage  area  information  that  will 
allow  the  computation  of  the  impedance  (or  similar  form)  matrix  for  the  network.  The 
calculation  of  the  impedance  matrix  actually  takes  place  in  a  separate  routine  if  the 
matrix  is  a  function  of  an  independent  variable  (e.g.,  frequency).  In  other  words,  every 
network  type  has  a  routine  associated  with  it  that  stores  in  the  master  storage  area  per¬ 
tinent  data  about  the  network,  and  calculates  the  impedance  matrix  if  the  matrix  is  not 
a  function  of  another  independent  variable.  However,  if  the  impedance  matrix  is  a 
function  of  an  independent  variable,  the  calculation  ot  the  matrix  takes  place  within 
another  routine  that  is  called  during  the  cycling  of  the  independent  variable. 


The  impedance  matrix  for  the  entire  network  is  found  by  effectively  pre-  and  post- 
multiplying  by  the  [C]  and  [D]  matrices  as  described  under  Theoretical  Description  of 
Multiport  Interconnection  Technique.  This  operation  is  performed  in  subroutine  BLDKKL. 
Actual  storage  and  matrix-multiplication  is  not  required  due  to  the  nature  of  the  [C]  and 
\D]  matrices. 

The  steps  described  above  are  the  basic  steps  that  are  performed  in  almost  all  control 
routines.  As  mentioned  previously,  the  user  can  perform  any  calculations  that  are  necessary 
and  supply  all  print  options. 

PROGRAM  CONVERSION  FROM  DOUBLE  PRECISION  TO  SINGLE 
PRECISION 

A  consistent  programming  convention  has  been  used  throughout  the  computer  pro¬ 
gram  in  order  to  facilitate  the  straightforward  conversion  from  double  precision  to  single 
precision.  The  following  two  types  of  changes  are  required: 

1 .  Cards  to  be  deleted  for  the  single  precision  version  contain  the  name 
of  the  routine  starting  in  Column  73. 

2.  Cards  to  be  changed  contain  the  name  of  the  routine  starting  in 
Column  73,  immediately  followed  by  an  *. 

Example 

Double  Precision 

SUBROUTINE  DCACSH  (DINR,  DINI,  DOUTR,  DOUTI) 

DOUBLE  PRECISION  DINR,  DINI,  DOUTR,  DOUTI,  DTOUTR,  DTOUTI  DCACSH 
DOUBLE  PRECISION  DCANG 1 ,  DSQRT,  DLOG  DCACSH 


DOUTR  =  DLOG  (DSQRT(DTOUTR*DTOUTR+DTOUTI*DTOUTI))  DCACSH* 


Single  Precision 

SUBROUTINE  DCACSH  (DINR,  DINI,  DOUTR,  DOUTI) 

delete 

delete 


DOUTR  =  ALOG(SQRT(DTOUTR*DTOUTR+DTOUTI*DTOUTI)) 


SUMMARY 


This  section  has  attempted  to  provide  a  description  of  the  computer  program  in 
terms  that  are  useful  for  non-programmers  as  well  as  programmers.  Due  to  the  general 
nature  of  the  program  package,  a  clear-cut  description  of  its  use  is  not  possible  without 
first  defining  a  problem.  For  this  reason,  several  examples  are  included  in  the  appendices 
to  show  the  user  how  to  construct  a  control  routine  for  the  solution  of  a  particular 
problem.  The  use  of  a  FORTRAN  control  routine  for  each  different  data  set  allows  the 
user  the  full  generality  of  the  programming  system  for  physical  simulation  of  his  specific 
problem. 


BIBLIOGRAPHY 


Notes  for  Seminar  on  Transducer  Array  Analysis  and  Evaluation ,  Transducer  Division, 
Naval  Undersea  Warfare  Center,  Intra-Division  Memo  No.  D604-160,  June,  1965. 

H.  J.  Carlin  and  A.  B.  Giordano,  Network  Theory:  An  Introduction  to  Reciprocal 
and  Nonreciprocal  Circuits ,  Prentice-Hall,  Lnglewood  Cliffs,  N.J.,  1964,  Chapter  3. 


APi’EHDIX  A 
Program  Description 


Introduction 

Presented  in  this  appendix  are  write-ups  describing  the  primary 
common  storage  blocks  and  all  subroutines  in  the  conpucer  program.  The 
primary  common  blocks  are  given  in  order  to  provide  easy  reference#  and  to 
avoid  having  to  describe  a  common  block  several  times.  The  subroutine 
write-ups  describe  the  basic  function^  or  the  routines,  and  define  the 
variables  in  the  parameter  lists- 


Packed  and  Unpacked  Network  Names 

Netvork  names  appear  throughout  the  program  in  two  forms,  referred  to 
as  packed  and  unpacked .  A  network  name  in  the  packed  form  is  one  that  has 
not  been  decoded  by  subroutine  UNPACK.  The  name  still  contains  port 
number,  total  number  of  ports  for  the  network,  and  the  piece  number  if 
applicable,  as  well  as  a  terminating  $.  An  unpa eked  name  is  simply  the 
two  letter  network  designation  assigned  by  the  user. 


A  -  1 


Description  of  the  Primary  Common  Blocks 


COMMON  /  /  DATA ( Si ze=MXDATA ) 

Blank  common  is  used  exclusively  as  the  master  storage  area  discussed 
in  Section  4.4. 

COMMON  /  DTA  TOT  /  MX  DATA 

MX  DATA  -  total  storage  dimensioned  for  master  storage  area 

COMMON  /  DTA  SUB  /  I  DTA  SB 

I  DTA  SB  -  counter  for  the  nvuiber  of  data  cells  which  have  been 
allocated  in  the  master  storage  area 

COMMON  /  CPT  WRD  /  N  BYT  P  W,  N  B  ?  BYT 

N  BYT  P  W  -  number  of  bytes  (1  byte=l  Hollerith  character)  per 
word  for  the  particular  machine  being  used 

N  3  P  BYT  -  number  of  bits  per  byte 

COMMON  /  PRT  OFF  /  NO  PRT 

NO  PRT  -  print  flag 

=  0  indicates  print 
^  1  suppresses  print 

COMMON  /  BASIC  8  /  D  F,  D  OMEGA,  D  PI 
D  F  -  frequency 

D  OMEGA  -  2 *pi* frequency 

D  PI  -  3 . 141b92653b89 
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CU-VV.-K  /  CATALG  /  NAME(IOO),  iiTR  SIZ(IOO),  I  PC (100), 

N  PC  ■PO'KIOU)  ,  I  TYPE (100)  ,  II  DT  GTR(IOO),  LOC(IOO),  I  CI..G  3D 

This  area  is  used  for  the  raster  storage  area  catalog. 


NAME 
MTR  SIZ 
I  PC 

N  PC  TOT 

I  TYPE 
N  DT  STR 

LOC 

I  CLG  SB 


-  network  name (unpacked) 

-  matrix  size 

-  piece  number  (0  if  not  applicable) 

-  total  number  of  pieces  with  this  network  name(0  if 
not  applicable) 

-  network  type(0  if  not  applicable) 

-  amount  of  data  storage  reauired  for  this  network  in 
addition  to  the  matrix  storage 

-  index  into  master  storage  where  information 
pertinent  to  this  particular  network  is  located 

-  total  number  of  entries  that  have  been  made  into  the 
catalog. 
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COMMON  /  PRT  T3L  /  NET  SP(100),  I  SIGN  C(100),  I  SIGN  D(100), 

J  PORT (100),  N  PORTS (100) ,  NET  NAM (100), 

MUM  APO(IOO),  NUM  NPO(IOO),  N  AHEAD(IOO), 

N  BACK (100),  LC  J  ROW(IOO),  N  PRT  SB 

This  area  is  used  for  the  interconnection  (or  port)  table- 

NET  SP  -  flaa  denoting  either  series  or  parallel  connection 

I  SIGN  C  -  sign  associated  with  voltage(or  force) 

I  SIGN  D  -  sign  associated  with  current  (or  velocity) 

J  PORT  -  the  network  port  number 

N  PORTS  -  total  number  of  ports  associated  with  this 
particular  network 

NEr  NAM  -  network  name (unpacked) 

NUM  APO  -  accessible  port  order  number-  Accessible  ports  are 
ordered  (NUM  APO  =  J  PORT  )  ,  however,  for  network 
ports  this  is  the  accessible  port  number  to  which 
the  network  port  is  connected 

NUM  NPO  -  network  port  order  number (negative  for  accessible 
ports) 

N  AHEAD  -  a  subscript  that  is  a  pointer  to  the  entry  in  the 

table  which  has  the  next  higher  network  port  order 
for  network  ports,  and  next  higher  accessible  port 
order  for  accessible  ports 

N  BACK  -  subscript  that  is  a  pointer  to  the  entry  that  has 

next  lower  network  port  order  for  network  ports,  and 
next  lower  accessible  port  order  for  accessible 
ports 

LC  J  ROW  -  starting  location  in  master  storage  for  port 
information 

N  PRT  SB  -  mmber  of  entries  currently  in  the  port  table- 
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Cu-i  iO!\ /  CA  2BY2  /I  AllR  ,1 A21R  ,1 A12R  ,1 A22R ,  I  AllI  ,1 A21I  ,IA12I  ,IA22I 


! 


Indices  for  referencing  conplex  quantities  of  a  2  by  2  A-form  matrix 

CCM40N/  CZ  2BY2  /IZ11R,IZ21R,1Z12R,IZ22R,  IZ11I  ,IZ21I  ,IZ12I  ,IZ22I 
Indices  for  referencing  conplex  ouantities  of  a  2  by  2  Z-form  matrix 

COrlMON/  CY  2BY2  /I YllR ,1 Y21R ,1 Y12R ,1 Y22R ,  IY11I  ,IY21I  ,IY121 ,1 Y22I 
Indices  for  referencing  conplex  quantities  of  a  2  by  2  Y-form  natrix 

CCMMON  /  CZ  3BY3  / 

IZ11R,  IZ21R,  IZ31R,  IZ12R,  IZ22R,  IZ32R,  IZ13R,  IZ23R,  IZ33R, 
IZ11I,  IZ21I,  IZ31I,  IZ12I ,  IZ22I ,  IZ32I ,  IZ13I,  IZ23I,  IZ33I 

Indices  for  referencing  conplex  quantities  of  a  3  by  3  Z-form  matrix 


CU»N  /  CY  3BY3  / 

I  YllR ,  IY21R,  I Y31R ,  IY12R,  IY22R,  IY32R,  1Y13R,  IY23R,  IY33R, 
IY11I t  IY21I,  IY31I,  I Y12I ,  IY22I,  IY32I,  IY13I,  IY23I,  IY33I 

Indices  for  referencing  conplex  quantities  of  a  3  by  3  Y-form  matrix 
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Subroutine  Descriptions 


Subroutine  descriptions  have  been  grouped  according  to  the 
functions  performed#  and  are  preceeded  by  an  alphabetical  index  to 
allow  easy  access  if  the  routine  name  is  known-  The  subroutines  have 
also  been  classified  to  aid  the  user  in  determining  priorities  in 
studying  tlie  write-ups-  The  classifications  are  as  follows- 

*  -  routines  for  which  the  purpose  should  be  understood 

by  the  user  for  operation  of  the  program , 

**  —  routines  which  do  not  reguire  detailed  understanding 

by  the  user# 

***  -  routines  that  are  used  internally  in  the  program, 

and  will  probably  not  be  directly  needed  by  the 
user- 
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Subroutine 


Page 


ALOCAT - 

BLDKKL - 

BLDMTK - 

BLKDTA - 

CASCDS - 

CHKSTR - 

CMPGAM - 

CMPZS  - 

DACMUL - 

DCACSH - 

DCANG1 - 

OCMDIV - 

DOMUL - 

DCPCT2 - 

OCRCIP - 

DCSHCH - 

DCSNCS - 

DECOM8 - 

DLETRC - 


- A_  26 

- A-110 

- A_  84 

- A_  89 

- A_  85 

- A_  76 

- A_  74 
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DJHCH  • 
F11FN  ■ 
FNDFDP- 
FNDFZV- 
FNDKRC' 
FNDLCJ' 
FNDMTR1 
FNDNAM' 
FNDNLT' 
FNDFRM' 
FNDTL  ■ 
FI1DTMS- 
FSRCH  ■ 
HEAD24- 
LA2BY2' 
LDAM2S- 
LD3C0L- 
LDDRIM- 
LDDT1  • 
LDDT11- 
LDDT14- 
LDDT15- 
LDDT2  ■ 
LDDT20* 


-A-  88 
-A-12U 
-A-126 
-A-113 
-A- 106 
-A-  18 
-A-  24 
-A-  25 
-A-  21 
-A- 12 5 
-A-  71 
-A-  70 
-A- 118 
-A-  94 
-A-  62 
-A-116 
-A- 127 
-A-  67 
-A-  28 
-A-  39 
-A-  43 
-A-  45 
-A-  30 
-A-  49 
-A-  51 


LDDT21 
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LDDT24 


A-  5  b 

LDDT25 - A_  57 

LDDT26 - A“  39 

- - 32 

- - 33 

LDDT7 - 34 

- - 37 

LDDT9 - 38 

LDFREQ - 16 

LDIMTR - 75 

- - 12 

LDS311 - 68 

LDZRVR - A- 114 

LOCMTR - 22 

LOCNAM - 23 

LSI  41 - - 69 

- - 65 

LY3BY3 - 66 

LZ2BY2 - 63 

LZ3BY3 - 64 

MAXCHK - - 

iUNCHK - A- 108 

MTR1-10V - 72 
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MTRi'lUI.' 


•A-  73 


PAGE  - 

PuLEO  - 

PRBFIT - 

PKCTLG - 

PREPZV - 

PRf-lARI - 

PRPTBL - 

PRT4BN - 

PRTMTR - 

PRTPRil - 

RDCPRT - 

RLIMMG - 

SKPLIN - 

SRCHBC - 

SRMIDB - 

STRLIN - 

TIME  - 

TP1415 - 

TYPCAL - 
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LDPOFT  * 


LDPORT  loads  into  the  interconnection  table  the  description  of  the 
network  ports  comprising  a  particular  accessible  port.  For  calls  to  this 
routine,  the  accessible  ports  should  be  defined  in  consecutive  and 
ascending  order. 

calling  sequence: 

CALL  LDPQRT(  IH  SP,  IPORT) 


input : 


IH  SP  -  flag  denoting  series  (S)  or  parallel  (P) 

interconnection  (left  justified  Hollerith  word). 

IPORT  -  a  string  of  characters  (terminated  by  $)  which 

describes  the  accessible  port  by  name,  and  each 
network  port  involved  in  the  connection. 

example: 

CALL  LDPQRT( *S* ,  ,8K8D,-+2D3A,1B2E,$* ) 

The  eighth  accessible  port  of  K8D  is  series  connected  to  the  second 
port  of  D3A  and  the  first  port  of  B?E. 

printout : 

Each  call  to  LDPORT  prints  the  information  being  loaded. 
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BLDMER  * 


BLDMER  constructs  the  matrix  relating  variables  at  the  accessible 
ports  of  a  network  for  which  an  interconnection  table  has  already  been 
loaded.  This  is  performed  by  automatically  chaining  successive  calls  to 
BLEKKL  and  RDCPRT.  This  matrix  may  be  of  size  less  than  or  equal  to  that 
size  as  loaded  in  the  port  table. 

calling  sequence; 

CALL  BLDMTR(LH  ANA) 

input: 

LH  ANA  -  packed  name  terminated  by  $  (two  computer  words),  of 
the  interconnected  network. 
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HLDKKL  ** 


BLDKKL  initiates  the  computation  of  the  individual  network  matrices 
(through  CASCDE),  and  using  these  matrices  computes  the  matrix  relating 
variables  at  the  accessible  ports  by  pre-  and  post-multiplying  by  the  C 
D  matrices  which  are  determined  by  Kirchoff's  laws  (without  actual 
storage  of  C  and  D  ). 

calling  sequence: 

CALL  BLDKKL (LH  ANA) 


LH  ANA  -  packed  name  terminated  by  $  (two  computer  words),  of 
the  interconnection  network. 


input: 


RDCPRT  ** 


RDCPRT  reduces  a  network  by  the  compaction  process  described  in 
Section  3. 

Bote:  after  reducing  AnA  to  AmA  do  not  reduce  AkA  to  AmA  for  k 
greater  than  m. 

calling  sequence : 

CALL  RDCPRT(LH  ANA,  LH  AMA) 

input: 

LH  ANA  -  packed  name  terminated  by  $  (two  computer  words),  of 
the  original  network. 

LH  AMA  -  packed  name  terminated  by  $  (two  computer  words),  of 
the  new  smaller  network. 
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LDFREX}  * 


LDFRB5  loads  the  frequency  into  a  connon  area  named  BASICS,  as  well 
as  the  corresponding  value  of  omega  (2.  x  pi  x  frequency). 

calling  sequence : 

CALL  LDra£Q(DDF) 


input: 


DDF  -  frequency 


printout : 

Value  of  frequency  (controlled  by  COMMON  /  PRT  OPF  /  NO  PRT) 
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UPPTBL  *** 


UPPTBL  updates  the  interconnection  table  with  network  port  order 
numbers,  and  forward  and  backward  pointers  to  network  ports. 

calling  sequence: 

CALL  UPPTBL 
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FNDLCJ  *** 


FUDLCJ  obtains  the  location  in  master  storage  for  all  matrices  needed 
for  a  particular  network,  and  allocates  storage  for  all  networks  for  which 
storage  has  not  been  set  aside.  Additional  storage  is  also  set  aside  for 
networks  that  will  automatically  be  reduced  in  size. 

calling  sequence : 

CALL  FNDLCJ (N  FIRST,  N  LAST) 

input: 

N  FIRST  -  entry  in  interconnection  table  to  begin  processing. 

N  LAST  -  entry  in  interconnection  table  to  end  processing. 
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ALOCAT  ** 


ALOCAT  allocates  storage  in  the  master  storage  area  by  inserting  a 
new  entry  into  the  catalog. 

calling  sequence: 

CALL  ALOCAT  (NAM  MTR,  M  NEW,  NMT  TOT) 

input: 

NAM  MTR  -  unpacked  name  of  network  to  be  added  to  catalog. 

M  NEW  -  size  of  matrix. 

NMT  TOT  -  total,  number  of  storage  locations  required. 
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CHKSTR  ** 


CHKSTR  checks  to  insure  that  the  current  total  master  storage  area 
has  not  exceeded  the  allowable  storage  limit,  MXDATA. 

calling  sequence : 

CALL  CHC  STR(I  DTA  CH) 

input: 

I  DTA  CH  -  counter  which  is  compared  to  total  allowable  storage 
limit,  MX  DATA. 
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FNDNET  ** 


FNDNET  searches  the  interconnection  table  for  a  particular 
interconnection  denoted  by  input  name,  and  returns  the  first  and  last 
entries  in  the  port  table  for  that  name. 

calling  sequence : 

CALL  FND  NET(LH  ANA,  N  FIRST,  N  LAST) 


input : 


LH  ANA  -  packed  name  terminated  by  $  (two  computer  words),  of 
the  interconnection  network. 


output ; 

N  FIRST  -  first  entry  for  this  interconnection. 
N  LAST  -  last  entry  for  this  interconnection. 
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LOCMTR  ** 


LOCMER  returns  the  location  in  master  storage  of  the  network 
specified  through  input  name  and  matrix  size. 

calling  sequence : 

I  =  L0CMTR(I  NAME,  I  MTR  SZ) 

input : 

I  NAME  -  network  name  (unpacked). 

I  MTR  SZ  -  matrix  size. 

output; 

I  -  master  storage  location  (  0  if  not  found). 
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LOCNAM  ** 


LOCNAM  returns  the  location  in  master  storage  of  the  first  network 
specified  by  input  name.  No  check  is  made  on  matrix  size. 

calling  sequence: 

I  =  LOCNAM(  INAME ) 

input: 

I  NAME  -  network  name  (unpacked) . 

output: 

I  -  master  storage  location  (  0  if  not  found). 
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FNDMTR  ** 


FNDMTR  returns  to  the  calling  program  the  entry  number  in  the  master 
storage  area  catalog  that  agrees  with  the  input  name  and  matrix  size. 

calling  sequence: 

CAUL  FNDMTR  (I  NAME,  I  MTR  SZ,  I  SUB) 


input : 

I  NAME 

-  network  name  (unpacked). 

I  MTR  SZ 

-  matrix  size. 

output : 

I  SUB 

-  subscript  in  master  storage  area  catalog  (  0  if  not 
found). 
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FNDNAM  ** 


FNDNAM  searches  the  catalog  for  the  first  occurence  of  an  input  name 
and  returns  the  entry  number  in  the  catalog.  Note  that  no  check  is  made 
on  matrix  size. 

calling  sequence: 

CALL  FNDNAM( I  NAME,  I  SUB) 

input: 

I  NAME  -  network  name  (unpacked). 

output : 

I  SUB  -  subscript  in  master  storage  area  catalog  (  0  if  not 
found). 
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CASCDE 


CASCDE  computes  the  individual  (typed)  network  matrices  through  a 
routine  called  TTPCAL,  and  cascdes  (or  multiplies)  all  2  by  2  matrices 
that  are  flagged  appropriately  and  stores  the  product  matrix  into  the 
primary  storage  area  for  that  section. 

calling  sequence: 

CALL  CASCDE(ISUB) 


input: 

I  SUB  -  entry  in  master  storage  catalog  for  the  product 
matrix. 

note:  An  additional  input  parameter  is  supplied  via  COMMDN/C  FUG  AZ/N  ZAY 

N  ZAY  -  flag  indicating  form  of  cascaded  matrix  desired. 

0  indicates  convert  A  to  Z  form  for  primary  area 

1  indicates  leave  in  A  form 

2  indicates  convert  A  to  Y  form  for  primary  area 
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TYPCAL  *** 


TYPCAL  checks  the  network  type  of  the  entry  in  the  catalog  specified 
by  input,  and  branches  to  the  proper  routine  to  calculate  the  matrix 
describing  the  network. 

calling  sequence: 

CALL  TYPCAL( IC  SUB) 

input: 

IC  SUB  -  catalog  entry  number  of  the  network  for  which  matrix 
computations  are  to  be  performed. 
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LDOTl  * 


LDOT1  loads  data  for  a  type  1  network  which  represents  a  straight 
piece. 

calling  sequence: 

CALL  LDHTl(LH  ANAI,  RHO,  AREA,  AID,  AOD,  PIECL,  CR,  Cl,  CM) 

input : 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

RHO  -  density. 

AREA  -  cross-sectional  area. 

AID  -  inner  diameter. 

AOD  -  outer  diameter. 

PIECL  -  length  of  single  piece. 

CR  -  real  part  of  velocity  of  sound  in  medium. 

Cl  -  imaginary  part  of  velocity  of  sound. 

CM  -  multiplier  for  complex  velocity  of  sound. 

note:  If  Cl  or  CM  is  non  zero,  then  the  missing  parameter  will  be  computed 
from  CM  =  CI/CR. 

printout : 

Values  of  input  data  ( controlled  by  COMMON  /  PRT  OFF  /  NO  PRT ) . 
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TYPE1  **+ 


TYPE1  computes  the  matrix  for  a  type  1  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPEl(NAMFR,  NNOTA) 

input: 

NAMTR  -  master  storage  location  for  matrix. 

MNDTA  -  master  storage  location  for  data. 


A  -  29 


LDDTT2  * 


LDDT2  loads  data  for  a  type  2  network  which  represents  a  conical 
piece. 

calling  sequence: 

CALL  LDDT2(LH  ANAI,RH0,  A1,AID1, A0D1,  A2,AID2,A0D2,  PIECL,  CR,CI,CM) 

input: 


LH  ANAI 

- 

piece  name  (packed)  terminated  by  $  (two  words). 

RHO 

- 

density. 

A1 

- 

cross-sectional  area  of  end  1. 

AUDI 

- 

inner  diameter  of  end  1. 

A0D1 

- 

outer  diameter  of  end  1. 

A2 

- 

cross-sectional  area  of  end  2. 

AID2 

- 

inner  diameter  of  end  2. 

AOD2 

- 

outer  diameter  of  end  2. 

PIECL 

- 

length  of  single  piece. 

CR 

- 

real  part  of  velocity  of  sound  in  medium. 

Cl 

- 

imaginary  part  of  velocity  of  sound. 

CM 

_ 

multiplier  for  complex  velocity  of  sound. 

note:  If  Cl  or  CM  is  non  zero,  then  the  missing  parameter  will  be  computed 
from  CM  =  CI/CR. 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  PRT  OFF  /  NO  FRT). 
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TYPE?  *** 


TYPE?  computes  the  matrix  for  a  type  2  network  as  a  function  of 
frequency. 

calling  sequence: 

CALL  TYPE?(NAMTR,  NNDTA) 

input: 

NAMIR  -  master  storage  location  for  matrix. 

NNDTA  -  master  storage  location  for  data. 
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LDDTU  * 


LDOTU  loads  data  for  a  type  4  network  which  represents  an  ideal 
transformer. 

calling  sequence: 

CALL  LDDT4(LH  ANAI,  SEC  TRN) 

input: 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

SEC  TRN  -  secondary  to  primary  turns  ratio. 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  PRT  OFF  /  NO  FRT). 
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LDDT5  * 


UDDT5  loads  data  for  a  type  5  network  which  is  an  arbitrary  matrix. 
calling  sequence : 

CALL  LDDT5(LH  ANAI,  ZR,  ZI) 

input : 

LH  ANAI  -  section  name  (packer!)  terminated  by  $  (two  words). 
ZR  -  real  part  of  Zij  array. 

ZI  -  imag  part  of  Zij  array. 

printout : 

Values  of  input  data  ( controlled  by  COMMON  j  PRT  OFF  /  NO  PRT ) . 
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LDDT7  * 


LDDT7  loads  data  for  a  type  7  network  which  is  a  rational  polynomial 
fit  to  an  impedance  or  mobility  matrix. 

calling  sequence: 

CALL  LDDT7(LH  ANAI,  FACTOR,  IFCRM) 

input: 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

FACTOR  -  scale  factor  for  all  coefficients. 

IFQRM  -  flag  to  determine  matrix  form, 

=  0  Z  form  (impedance) 

>  1  Y  form  (admittance)  *=  inverse  of  Z  form. 

data  card  input  for  polynomial  fit  coefficients : 

Card  1  (  FORMAT(3I5,  l6AL)  ) 

N  CRD  N  -  order  of  numerator  polynomial. 

N  GRD  D  -  order  of  denominator  polynomial. 

N  SYM  KJ  -  symmetry  flag  for  loading  matrix  coefficients. 

=  -1  indicates  negative  symmetry  across  diaginal. 

=  0  indicates  no  symmetry. 

=  +1  indicates  synmetry  across  diagonal. 

note:  if  synmetry  exists  then  only  coefficients  for 
elements  in  upper  triangle  are  needed  and  should  be 
ordered  by  columns  as  (l,l),  (l,2),  (2,2), 

(1,3),...,  (N,N).  However,  if  no  synmetry  exists 

then  the  order  should  be  (l,l),  (2,l),  ...»  (N,l), 

(l,2),  (2,2),  .... 

IHEAD  -  l6  word  alphanumeric  heading  to  identify 
coefficients. 

Cards  2  -  NQRDN+1  (  P0RMAT(E17.10,  IX,  E17.10)  ) 

DCR  -  real  part  of  numerator  coefficient  in  ascending 

powers  of  the  independent  variable. 

-  Imaginary  part  of  numerator  coefficient  in  ascending 
powers  of  the  independent  variable. 
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DC  I 


Cards  H0RDR*2  -  NORDN+NORDD+1  (  P0RMAT(E17.1O,  IX,  E17.10)  ) 

DCR  -  real  part  of  denominator  coefficient  in  ascending 

powers  of  the  independent  variable. 

DC  I  -  imaginary  part  of  denominator  coefficient  in 

ascending  powers  of  the  independent  variable. 


printout : 

Values  of  input  data  (controlled  by  C0MM3N  /  PRT  OFF  /  NO  PRT). 
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TTPE7  *** 


TYPE7  computes  the  matrix  for  a  type  7  network  as  a  function  of 
frequency. 

calling  sequence: 

CALL  TYPE7(NAMTR,  NNDTA) 

input : 

NAMTR  -  master  storage  location  for  matrix. 

NNDTA  -  master  storage  location  for  data. 
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LDDT8  * 


LDDT8  loads  data  for  a  type  8  network  which  is  a  voltage  to  voltage 
converter . 


calling  sequence: 

CALL  LDDT8(LH  ANAI,  ZIR,  Zn,  ZGR,  ZGI,  ALPR,  AIPI) 

input : 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 
ZIR  -  real  part  of  amplifier  input  impedance  (ohms). 

ZH  -  imag  part  of  amplifier  input  impedance  (ohms). 

ZGR  -  real  part  of  amplifier  output  impedance  (ohms). 

ZGI  -  imag  part  of  amplifier  output  impedance  (ohms). 

ALFR  -  real  part  of  amplification  factor. 

ALPI  -  imag  part  of  amplification  factor. 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  PRT  OFF  /  NO  FRT). 
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LDET9  * 


LDDT9  loads  data  for  a  "type  9  network  which  is  a  voltage  to  current 
converter. 

culling  sequence : 

CALL  LDDT9(LH  ANAI,  ZIR,  ZII,  ZCR,  ZOI,  BETAR,  BETAl) 

input: 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

ZIR  -  real  part  of  amplifier  input  impedance  (ohms). 

ZII  -  i™»g  part  of  amplifier  input  impedance  ( ohms ) . 

ZQR  -  real  part  of  amplifier  output  impedance  (ohms). 

ZOI  -  im»g  part  of  amplifier  output  impedance  (ohms). 

BETAR  -  real  part  of  voltage  to  current  gain  (mhos). 

BETAl  -  imag  part  of  voltage  to  current  gain  (mhos). 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  PRT  OFF  /  NO  PRT). 
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LDDT11  * 


LDDT11  loads  data  for  a  type  11  network  which  is  a  3-3  mode  ceramic 
stack. 

culling  sequence: 

CALL  LDtDTll(LH  ANAI,  NPCS,  RHO,  AREA,  AID,  AOD,  PIECL, 
DE33TR,DE33TI,DE33TM,  DG33R,DG33I,DG33M,  DS33DR,DS33DI,DS33DM) 

input: 

LH  ANAI  -  section  name( packed)  terminated  by  $  (two  words). 
NPCS  -  total  number  of  pieces  for  this  ceramic  section. 

RHO  -  density. 

AREA  -  cross-sectional  area. 

AID  -  inner  diameter. 

AOD  -  outer  diameter. 

PIECL  -  length  of  single  piece. 

DE33TR  -  real  part  of  E33T  (always  double  precision). 

DE33TI  -  irnag  part  of  E33T  (always  double  precision). 

DE33IM  -  E33T  multiplier  (always  double  precision). 

DG33R  -  real  part  of  G33R  (always  double  precision). 

DG33I  -  imag  part  of  G33I  (always  double  precision). 

DG33M  -  G33  multiplier  (always  double  precision). 

DS33DR  -  real  part  of  S33D  (always  double  precision). 

DS33DI  -  imag  part  of  S33D  (always  double  precision). 

DS33DM  -  S33D  multiplier  (always  double  precision). 

note:  For  ceramic  parameters  if  either  I  or  M  is  zero,  the  missing 
parameter  will  be  computed  by  M—  I/R. 
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Values  of  input  data  (controlled  by  CQMMaSf/PRT  OFF/NO  PRT) . 
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TYPE11  *** 


TfPEll  computes  the  matrix  for  a  type  11  network  as  a  function  of 
frequency. 

calling  sequence; 

CALL  TYPEll(NAMTR,  NNDTA) 

input; 

NAMTR  -  master  storage  location  for  matrix. 

NNDTA  -  master  storage  location  for  data. 
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TllSBR  *** 


T11SBR  performs  operations  required  by  type  11,  and  stores  into 
c cannon  areas  used  in  type  11  calculations. 

calling  sequence : 


CALL  TllSBR 


LDDT14  * 


LDDTlU  loads  data  for  a  type  14  network  which  is  a  3-1  mode  ceramic 

bar. 

calling  sequence: 

CALL  LDDT14(LH  ANAI,  RHO,  THICK,  FOIL  W,  PIEC  L, 
DE33TO,DE33TI,DE330M,  DG31R,DG3H,DG31M,  DS11ER,DS11EI,DS11EM) 

input : 

LH  ANAI  -  section  name( packed)  terminated  by  $  (two  words). 

RHO  -  density. 

THICK  -  thickness. 

FOIL  W  -  foil  width. 

PIECL  -  length  of  single  piece. 

DE33TR  -  real  part  of  E33T  (always  double  precision). 

DE33TI  -  imag  part  of  E33T  (always  double  precision). 

DE33TW  -  E33T  multiplier  (always  double  precision). 

DG3IR  -  real  part  of  G3IR  (always  double  precision). 

DG31I  -  imag  part  of  G31I  (always  double  precision). 

DG31M  -  G31  multiplier  (always  double  precision). 

DS11ER  -  real  part  of  SHE  (always  double  precision). 

DS11EI  -  imag  part  of  S11E  (always  double  precision). 

DSHEM  -  SHE  multiplier  (always  double  precision). 

note:  For  ceramic  parameters  if  either  I  or  M  is  zero,  the  missing 
parameter  will  be  computed  by  M=-  I/R. 

printout: 

Values  of  input  data  (controlled  by  COMMON  j  PRT  OFF  f  NO  PRT). 
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TfPElt  *** 


TYPElU  computes  the  matrix  for  a  type  It  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TfPElL(NZMTR,  NNDTA) 

input: 

NZMTR  -  master  storage  location  for  matrix. 

NNDTA  -  master  storage  location  for  data. 


A 


44 


LDDT15  * 


LDDT15  loads  data  for  a  type  15  network  which  is  a  5-1  mode  ceramic 

tube. 

calling  sequence: 

CALL  LDDT15(LH  ANAI,  RHO,  R  IN  D,  R  OUT  D,  PIEC  L, 
DE35rm,DE33TI,DE33TM,  DG31R,DG31I,DG31M,  DS11ER,DS11EI,DS11EM) 

input: 

LH  ANAI  -  section  name( packed)  terminated  by  $  (two  words). 

RHO  -  density. 

R  IN  D  -  inner  diameter. 

R  OUT  D  -  outer  diameter. 

PIECL  -  length  of  single  piece. 

DE33TR  -  real  part  of  E33T  (always  double  precision). 

DE33TI  -  imag  part  of  E33T  (always  double  precision). 

DE33IM  -  E33T  multiplier  (always  double  precision). 

DG31R  -  real  part  of  G3IR  (always  double  precision). 

DG3H  -  imag  part  of  G3II  (always  double  precision). 

DG3IM  -  G3I  multiplier  (always  double  precision). 

DS11ER  -  real  part  of  SHE  (always  double  precision). 

DS11EI  -  imag  part  of  SHE  (always  double  precision). 

DS11EM  -  SHE  multiplier  (always  double  precision). 

note:  For  ceramic  parameters  if  either  I  or  M  is  zero,  the  missing 
parameter  will  be  computed  by  M=-I/R. 

printout : 

Values  of  input  data  (controHed  by  CQMMDN  /  PRT  OFF  /  NO  PRT). 


A  -  L5 


TYPE15  *** 


TYPE15  computes  the  matrix  for  a  type  15  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPE15(NZMTB,  NNDTA) 

input : 

NZMPR  -  master  storage  location  for  matrix. 

NNDTA  -  master  storage  location  for  data. 


-  46 


A 


T1U15S  *** 


T1415S  performs  operations  required  by  type  lU  and  type  15,  and 
stores  into  common  areas  used  in  type  lU  and  type  15  calculations. 

calling  sequence; 

CALL  T1415S 


A  -  1+7 


TP1U15  *** 


TPl4l5  performs  operations  required  by  type  lU  and  type  15  networks. 
calling  sequence: 

CALL  H>l4l5 


LDDT20  * 


LDDT20  loads  data  for  a  type  20  network  which  is  a  linear 
transformer. 

calling  sequence : 

CALL  LDDT20(LH  ANAI,  RP,  HIP,  R0,  HLO,  RS,  HLS,  PHI) 

input; 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

RP  -  copper  loss  in  primary. 

HLP  -  winding  Inductance  in  primary. 

RO  -  core  loss. 

HLO  -  magnetizing  Inductance. 

RS  -  copper  loss  in  secondary. 

HT-S  -  winding  inductance  in  secondary. 

PHI  -  secondary  to  primary  turns  ratio. 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  PRT  OPT  /  NO  PRT). 


A  -  U9 


TYPE20  *** 


TYPE20  computes  the  matrix  for  a  type  20  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPE20(  NAMER ,  HMDTA) 

input: 

NAMTR  -  master  storage  location  for  matrix. 

NNDTA  -  master  storage  location  for  data. 
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IX0T21  * 


UDT21  loads  data  for  a  type  21  network  which  is  a  series  impedance. 
CftlUng  sequence ; 

CALL  LDDT21(LH  ANAI,  RS,  HLS,  CS) 

input: 

T.H  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

RS  -  resistance  (Ohms). 

HLS  -  inductance  (Ry). 

CS  -  capacitance  (Fy). 

printout : 

Values  of  input  data  (controlled  by  COWON  /  PRT  OFT  /  NO  FRT). 


A  -  51 


TYPE21  *** 


TTPE21  computes  the  matrix  for  a  type  21  network  as  a  function  of 
frequency. 

calling  sequence: 

CALL  TYPE21(NAMER,  NHOTA) 

input: 

NAMTR  -  master  storage  location  for  matrix. 

NNOTA  -  master  storage  location  for  data. 
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IDDT22  * 


LDDT22  loads  data  for  a  type  22  network  which  is  a  parallel 
impedance . 

calling  sequence : 

CALL  LDDT22(LH  ANAI,  RS,  HLS,  CS) 

input: 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

RS  -  resistance  (Ohms). 

HLS  -  inductance  (Hjr). 

CS  -  capacitance  (py). 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  FRT  OFF  /  NO  PRT). 


A  -  53 


T5fPE22  *** 


TJTPE22  computes  the  matrix  for  a  type  22  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPE22(N  MPR,  N  OTA) 

input: 

N  MTR  -  master  storage  location  for  matrix. 

N  DTA  -  master  storage  location  for  data. 


A  -  5^ 


LEDT24  * 


LDOTpU  loads  data  for  a  type  24  network  which  is  a  series 
admittance. 

calling  sequence: 

CALL  LDDT24(LH  ANAI,  RP,  HLP,  CP) 

input: 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

RP  -  resistance  (Ohms). 

HLP  -  inductance  (Ry). 

CP  -  capacitance  (Fd). 

printout : 

Values  of  input  data  ( controlled  by  CQMMDN  /  PRT  OFF  /  NO  PRT). 


A  -  55 


TYPE2U  *** 


TYPE2U  computes  the  matrix  for  a  type  24  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPE24(RAMTR,  RKDTA) 

input: 

NAMPR  -  master  storage  location  formatrix. 

NKOTA  -  master  storage  location  for  data. 
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LDET25  * 


LDDT25  loads  data  for  a  type  25  network  which  is  a  parallel 
admittance. 

calling  sequence : 

CALL  LDDT25 ( LH  ANAI,  RP,  HLP,  CP) 

input: 

LH  ANAI  -  piece  name  (packed)  terminated  by  $  (two  words). 

RP  -  resistance  (Ohms). 

HLP  -  inductance  (Ify) . 

CP  -  capacitance  (Fd). 

printout : 

Values  of  input  data  (controlled  by  COMMON  /  PRT  OFT  /  NO  PRT). 


A  -  57 


TYPE25  *** 


TYPE25  computes  the  matrix  for  a  type  25  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPE25(»  MR,  N  DTA) 

input; 

N  MER  -  master  storage  location  for  matrix. 

N  OTA  -  master  storage  location  for  data. 
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LDDT26  * 


LDDTp6  loads  data  for  a  type  26  network  which  is  a  differential  input 
operational  amplifier. 

calling  sequence ; 

CALL  LDnT26(LH  ANAI,  R1,C1,  R2,C2,  RD,CD,  R3,HL3, 

0MEGA1,  0MEXJA2,  0MEXJA3 ) 

input ; 

LH  ANAI  -  section  name( packed)  terminated  hy  $  (two  words). 

R1  -  input  resistance  (Ohms). 

Cl  -  input  capacitance  (farads). 

R2  -  input  resistance  (Ohms). 

C2  -  input  capacitance  (farads). 

RD  -  differential  input  resistance  (Ohms). 

CD  -  differential  input  capacitance  (farads). 

R3  -  output  resistance  (Ohms). 

C3  -  output  capacitance  (farads). 

HL3  -  output  inductance  (henries). 

AO  -  dc  open  loop  gain. 

OMEGA1  -  first  pole  frequency  of  open  loop  transfer  function 

(l/sec). 

OMEGA2  -  second  pole  frequency  of  open  loop  transfer  function 
(l/sec). 

0MBSA3  -  third  pole  frequency  of  open  loop  transfer  function 
(l/sec). 

printout: 

Values  of  input  data  (controlled  by  CCMM3N  /  PRT  OFT  /  NO  PRT). 


A  -  59 


TYPE26  *** 


TYPE26  computes  the  matrix  for  a  type  2 6  network  as  a  function  of 
frequency. 

calling  sequence : 

CALL  TYPE26(W  MER,  N  DTA) 

input: 

N  MHR  -  master  storage  location  for  matrix. 

N  DTA  -  master  storage  location  for  data. 


A 


60 


AIDOD  ** 


AIDOD  calculates  either  inner  diameter,  outer  diameter,  or  area  given 
the  proper  combination  of  the  other  two. 


calling  sequence: 

CALL  AID0D( A, 

AID,  A0D) 

input/output  ( see  tables  below ) : 

A 

-  area 

AID 

-  inner  diameter 

A0D 

-  outer  diameter 

The  following  table  shows  the  input  quantities  and  output  quantities 
for  every  combination  of  input. 


_ input _ 

Area  ID  OD 

~T~  “0  ~0 

0X0 
0  0  X 

X  X  0 

X  0  X 

0  X  X 

XXX 
0  0  0 


where  X  -  non- zero  input 
C  -  calculated 


output 

Area  TO  GD 

~T~  ~0  ~0 

(insufficient  data) 
COX 
X  X  C 

X  C  X 

C  X  X 

XXX 
(insufficient  data) 


caution:  a  zero  input  should  be  passed  through  a  variable  name,  and  not  as 
a  literal  in  the  subroutine  call  since  the  value  may  be  changed  in  the 
subroutine. 


The  following  relationships  are  used  in  the  computation  and  are  valid 
for  types  1,  2,  and  11. 

A  *=  (pi/U  )(0D**2  -  H>**2) 

ID  -  SQRT(0D**2  -  UA/pi) 

0D  «=  SQJRT (U A/pi  -  ID**2) 


A  -  6l 


LA2BY2  ** 

LA2BJT2  loads  indices  for  master  storage  for  the  A-form  of  a  2  port 
network  (2x2  matrix)  into  cannon  area  CA2BY2. 

calling  sequence : 

CALL  LA2BST2 

input: 

N  MIR  -  master  storage  location  for  matrix 


A  -  62 


LZ2BY2  ** 


LZ2Bf?  loads  indices  for  master  storage  for  the  Z-fora  of  a  2  port 
network  (2x2  matrix)  into  c cannon  area  CZ2B3T2. 

calling  sequence : 

CALL  LZ2BY2 

input: 

N  MTR  -  master  storage  location  for  matrix 


A  -  63 


LZ3BY3  ** 


LZjBY^  loads  indices  for  master  storage  for  the  Z-form  of  a  3  port 
network  (3x3  matrix)  into  comnon  area  CZ3BY3. 

calling  sequence : 

CALL  LZ3BY3 

input: 

N  MTR  -  master  storage  location  for  matrix 


A  -  6k 


LY2BJf2  ** 


LY2BY2  loads  indices  for  master  storage  for  the  Y-form  of  a  2  port 
network  (2x2  matrix)  into  cannon  area  CY2BST2. 

calling  sequence: 

CALL  LY2BY2 

input : 

N  MTR  -  master  storage  location  for  matrix 


A  -  65 


LY3®C3  ** 


LX3BY3  loads  indices  for  master  storage  for  the  Y-form  of  a  3  port 
network  (3x3  matrix)  into  common  area  CY3BY3. 

cnlllng  sequence ; 

CALL  IY3BY3 

Input: 

N  MTR  -  master  storage  location  for  matrix 


A  -  66 


LDDRIM  *** 


LDDRIM  loads  into  the  master  storage  area  the  real,  imaginary,  and 
multiplier  supplied  on  input.  Note  that  only  the  real(R)  and  either  the 
imaginary(I)  or  multiplier(M)  or  both  need  be  supplied  and  the  missing 
parameter ( input  as  zero)  will  be  computed  by  M*-I/R.  The  parameters  are 
stored  in  sequence  using  the  present  value  of  IOTASB  which  is  the  total 
number  of  storage  cells  used  in  master  storage. 

calling  sequence : 

CALL  LDDRIM(DR,  DI,  DM) 


input : 


DR 


real  part. 


DI 

input /output : 

DM 


-  imaginary  part. 


-  multiplier  (DM  =  -DI/ER ) . 
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LDSB11  *** 

LDSB11  sets  up  a  group  of  indices  that  enables  easy  referencing  of 
network  type  11  (3-5  mode  ceramic  stack)  data.  These  indices  are  stored 
in  cannon  block  C0MN11. 

calling  sequence : 

CALL  LDSBll(NDTA) 

input: 

RDTA  -  master  storage  location  where  data  starts. 
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LS1U15  *** 


LS1U15  sets  up  a  group  of  Indices  that  enables  easy  referencing  of 
network  type  lU(j-l  mode  ceramic  bar)  and  type  15(3-1  mode  ceramic  tube) 
data.  These  indices  are  stored  into  common  block  CM1415. 

calling  sequence : 

CALL  LS1415(NDTA) 

input; 

MDTA  -  master  storage  location  where  data  starts. 


A  -  69 


FNOTMS  ** 


PTTOTMS  determines  the  total  mass  of  all  type  1,  type?,  and  typell 
networks. 

calling  sequence ; 

CALL  FNOTMS 

printout; 

mass  in  MKS  and  English  units  (controlled  by  CCMON  /  PRT  OFF  /  NO 
PRT). 


A  -  70 


FNDTL  ** 


FNOTL  determines  the  total  length  of  all  networks  specified  by  input 
name  for  types  1,  2,  11,  14,  and  15. 

calling  sequence : 

CALL  FMDTL(L  NAMES,  N  NAMES) 

input ; 

L  NAMES  -  an  array  containing  the  names (unpacked)  of  sections 
to  be  used  in  total  length  computation. 

N  NAMES  -  number  of  names  in  the  array  L  NAMES. 

printout : 

value  of  length  (  controlled  by  C0M40N  /  PRT  OFT  /  NO  PRT). 


A  -  71 


MTRMOV  ** 


MPRMOV  moves  a  matrix  of  given  size  from  one  area  in  master  storage 
to  another  area  in  master  storage. 

calling  sequence: 

CALL  MPRMCf\r(NHEW,  NOLD,  MERSIZ) 


input: 


KNEW  -  master  storage 

moved. 

HOLD  -  master  storage 

hfPRSIZ  -  matrix  size. 


location  where  matrix  is  to  be 

location  where  matrix  is  stored. 
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KTRMLJL  ** 

MTRMUL  multiplies  two  matrices  stored  in  master  storage. 
calling  sequence: 

CALL  MTRMUL(LCHEW,  LCOLDA,  LCOLDB,  MA,  HA,  MB,  NB) 
where  (matrix  HEW)  =  (matrix  A)  x  (matrix  B) 

input: 

LCHEW  -  master  storage  location  where  product  is  to  be 
stored. 

LCOLDA  -  master  storage  location  where  matrix  A  is  stored. 

LCOLDB  -  master  storage  location  where  matrix  B  is  stored. 

MA  -  number  of  rows  of  matrix  A. 

HA  -  number  of  columns  of  matrix  A. 

MB  -  number  of  rows  of  matrix  B. 

NB  -  number  of  columns  of  matrix  B. 

note:  LCNEW  cannot  coincide  with  either  LCOLDA  or  LCOLDB,  and  HA  must 
equal  MB. 
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DLE3RC  ** 


DLETOC  deletes  specified  rows  and/or  columns  from  a  matrix  and 
performs  a  scalar  multiplication  on  the  matrix. 


calling  sequence; 

CALL  DLEERC(NN  OLD  R,  NN  OLD  I,  M  OLD,  K  OLD,  I  OUT,  J  OUT, 
NI  OUT,  NJ  OUT,  N  NEW  R,  N  NEW  I,  SCALAR) 


input: 


NN  OLD  R 

NN  OLD  I 

M  OLD 
N  OLD 
I  OOT 
J  OOT 
NI  OOT 
NJ  OOT 
N  NEW  R 

N  NEW  I 

SCALAR 


-  master  storage  location  of  real  part  of  input 
matrix. 

-  master  storage  location  of  imaginary  part  of  input 
matrix. 

-  number  of  rows  in  input  matrix. 

-  number  of  columns  in  input  matrix. 

-  array  containing  indices  of  rows  to  be  deleted. 

-  array  containing  indices  of  columns  to  be  deleted. 

-  number  of  rows  to  be  deleted. 

-  number  of  columns  to  be  deleted. 

-  master  storage  location  of  real  part  of  output 
matrix. 

-  master  storage  location  of  imaginary  part  of  output 
matrix. 

-  scale  factor  for  output  matrix. 
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LDIMER  ** 


LDIMIH  loads  an  identity  matrix  of  specified  size  into  the  desired 
master  storage  area. 

calling  sequence: 

CALL  LD3MPR(ISUB,  MSIZ) 

input : 

ISUB  -  master  storage  location  for  unit  matrix. 

MSIZ  -  matrix  size. 
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DBCQMS  ** 


DECOMB  solves  the  system  of  equations  AX  =  B  for  the  vector  X  (INVFG 
=  0)  or  computes  the  inverse  of  the  matrix  A  (INVFG  =1). 

calling  sequence; 

CALL  DECCM8(DA,NA,  DAI,NAI,  DAINV,NAINV,  DAINVI,NAINVI, 

NSIZE,  MDFA,  INVFG) 


input: 


output: 


input: 


output: 


DA 

NA 

DAI 

NAI 

DAXNV 

NAINV 

DAINVI 


input: 

NAINVI 

NSIZE 

MDFA 


INVFG 


array  containing  the  real  part  of  matrix  A. 
starting  subscript  in  DA. 
array  containg  imaginary  part  of  matrix  A. 
starting  subscript  in  DAI. 

array  for  storing  the  real  part  of  A  inverse. 

subscript  in  DAINV  to  begin  storing. 

array  for  storing  Imaginary  part  of  A  inverse. 

subscript  into  DAINVI  to  begin  storing. 

size  of  matrix  A. 

matrix  decomposition  flag, 

=  0,  decompose  matrix. 

=  1,  matrix  has  already  been  decomposed. 

either  solve  for  vector  X  (INVFG  =  0),  or  compute 
inverse  of  A  (INVFG  =1). 
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ATOZ  ** 


ATOZ  converts  a  2x2  matrix  in  the  A- form  to  the  Z-form. 
calling  sequence : 

CALL  ATOZ(N  AMPR,  N  ZMIR) 

input: 

N  AMIR  -  master  storage  location  for  A-form. 

N  ZMTR  -  man  ter  storage  location  for  Z-form. 

transformation : 

Zll  =  All  /  A21 

Z12  =  A12  -  All  A22  /  A21 

Z21  =  1  /  A21 

Z22  =  -A22  /  A21 
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ZTOA  ** 


ZTOA  converts  a  2x2  matrix  in  the  Z-form  to  the  A- form. 
calling  sequence : 

CALL  ZTOA(N  ZMTR,  N  AMTR ) 

input : 

N  ZMTR  -  master  storage  location  for  Z-form. 

N  AMUR  -  master  storage  location  for  A-form. 

transformation : 

All  =  Zll  /  Z21 
-  A12  «  Z12  -  Zll  Z22  /  Z21 
A21  =  1  /  Z21 
A22  =  -Z22  /  Z21 
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ATOT  ** 


ATOT  converts  a  2x2  matrix  in  the  A- form  to  the  Y-form. 
calling  sequence ; 

CALL  ATOT(N  AMTR,  N  YMTR ) 

input: 

N  AMTR  -  master  storage  location  for  A- form. 

N  YMTR  -  master  storage  location  for  Y-form. 

transformation : 

Yll  =  A22  /  A12 

Y12  =  (A12  A21  -  All  A22)  /  A12 
Y21  =  1  /  A12 
Y22  =  -All  /  A12 
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YTOA  ** 


YTOA  converts  a  2x2  matrix  in  the  Y-form  to  the  A- form. 
calling  sequence : 

CALL  YTOA(N  YMTR,  N  AMTR) 

input; 

N  YMTR  -  master  storage  location  for  Y-form. 

N  AMTR  -  master  storage  location  for  A-form. 

transformation ; 

All  =  -Y22  /  Y21 
A12  ■  1  /  Y21 

A21  =  (Y12  Y21  -  Yll  Y12)  /  Y21 
A22  =  Yll  /  Y21 
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DACMUL  ** 


DACMJL  performs  a  complex  accumulative  multiplication. 
calling  sequence: 

CALL  DACMUL(DA,  DB,  DC,  DD,  DE,  DF) 

input: 

DA  -  real  part  of  first  input  number. 

DB  -  imag  part  of  first  input  number. 

DC  -  real  part  of  second  input  number. 

DD  -  imag  part  of  second  input  number. 

input /output : 

DE  -  real  part  of  resultant. 

DF  -  imag  part  of  resultant, 

the  following  steps  are  performed: 

DE  =  DE  +  DA  *  DC  -  DB  *  DD 

DF  =  DF  +  DA  *  DD  +  DB  *  DC 

note:  resultant  variables  may  occupy  the  same  locations  as  either  input. 
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DCMMUL  ** 


DCMMUL  performs  a  complex  multiplication. 
calling  sequence: 

CALL  DCMMUL(DA,  DB,  DC,  DD,  DE,  DF) 

input: 

Da  -  real  part  of  multiplicand. 

DB  -  iniag  part  of  multiplicand. 

DC  -  real  part  of  multiplier. 

DD  -  imag  part  of  multiplier. 

output : 

DE  -  real  part  of  product. 

DF  -  •’’•nag  part  of  product, 

note:  the  result  may  occupy  the  same  locations  as  either  input. 


A  -  8? 


DCMDIV  ** 


DCMDIV  performs  a  complex  division. 
calling  sequence; 

CALL  DCMDIV  (DA,  DB,  DC,  DD,  DE,  DF) 

input; 

DA  -  real  part  of  dividend. 

DB  -  imag  part  of  dividend. 

DC  -  real  part  of  divisor. 

DD  -  imag  part  of  divisor. 

output ; 

DE  -  real  part  of  quotient. 

DF  -  imag  part  of  quotient, 

note;  the  result  may  occupy  the  same  locations  as  either  input. 


A  -  83 


DCRCIP  ** 


DCRCIP  computes  the  reciprocal  of  a  complex  number. 
calling  sequence; 

CALL  DCRC3P(DINR,  DINI,  DOUTR,  DOUTl) 

input: 

DINR  -  real  part  of  input. 

DINI  -  imag  part  of  input. 

output : 

DOUER  -  real  part  of  output. 

DOUTI  -  -tmag  part  of  output, 

note:  the  result  may  occupy  the  same  locations  as  either  input. 


A  -  8U 


DC  SQRT  ** 


DC SQRT  computes  the  square  root  of  a  complex  argument. 
calling  sequence: 

CALL  DC3QRT(DINR,  DINI,  DCSQTR,  DCSQTI) 

Input: 

DINR  -  real  part  of  argument. 

DINI  -  imag  part  of  argument. 

output : 

DCSQTR  -  real  part  of  square  root. 

DCSQTI  -  Imag  part  of  square  root, 

note:  the  result  may  occupy  the  same  locations  as  the  input. 


A  -  85 


DCANGl  ** 


DCANG1  returns  the  phase  (in  radians)  of  a  complex  number  adjusted  so 
that  the  result  appears  between  -pi  and  pi. 

calling  sequence : 

ANG  =  DCANGi(DX,  DY) 

where:  AUG  =  ARC  TAN(DY/DX) . 
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DCACSH  ** 


DCACSH  calculates  the  complex  inverse  hyperbolic  cosine,  Y  *  ARC 
COSH(Z). 

where:  real  Y  =  ln( |Z  +  SQRT(Z**?-l) I ) 
imag  Y  =  ang(Z  +  SQRT(Z**?-l) ) 

CALLING  SEQUENCE: 

CALL  DCACSH(DINR,  DINI,  DOUTR,  DOUTI) 

input: 

DINR  -  real  part  of  input,  Z. 

DINI  -  imag  part  of  input,  Z. 

output : 

DOUTR  -  real  part  of  output,  Y. 

DOUTI  -  imag  part  of  output,  Y. 

note:  the  result  may  occupy  the  same  locations  as  the  input. 
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DSHCH  ** 


DSHCH  computes  the  hyperbolic  sine  and  hyperbolic  cosine  of  an 
argument. 

calling  sequence ; 

CALL  DSHCH(DARG,  DSINH,  DCOSH) 

input: 

DARG  -  input  argument. 

output: 

DSINH  -  hyperbolic  sine  of  DARG. 

DCOSH  -  hyperbolic  cosine  of  DARG. 
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DCSNCS  ** 


DCSNCS  computes  the  sine  and  cosine  of  a  complex  argument. 
calling  sequence; 

CAUL  DCSNCS(DARGR,DARGI,  DSINR,DSINI,  DC0SR,DC0SI) 

input; 

DARGR  -  real  part  of  input  argument. 

DARGI  -  imag  part  of  input  argument. 

output ; 

DSINR  -  real  part  of  sine. 

DSINI  -  imag  part  of  sine. 

DCOSR  -  real  part  of  cosine. 

DCOSI  -  imag  part  of  cosine, 

note:  either  result  may  occupy  the  same  locations  as  the  input. 
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DCSHCH  ** 


DCSHCH  computes  the  hyperbolic  sine  and  cosine  of  a  complex 
argument. 

calling  sequence : 

CALL  DCSHCH(DARGR,DARGI,  DSDfHR ,DSINHI,  DCOSHR,DCOSHl) 

input: 

DARGR  -  real  part  of  input  argument. 

DARGI  -  imag  part  of  input  argument. 

output: 

DSINHR  -  real  part  of  hyperbolic  sine. 

DSINHI  -  imag  part  of  hyperbolic  sine. 

DCOSHR  -  real  part  of  hyperbolic  cosine. 

DCOSHI  -  imag  part  of  hyperbolic  cosine, 

note:  either  result  may  occupy  the  same  locations  as  the  input. 
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PRTMTR  ** 


PRTMTR  prints  a  square  matrix  that  is  in  master  storage. 
calling  sequence: 

CALL  PRTMTR ( N  MTR,  MTR  SIZ) 

input: 

N  MTR  -  master  storage  location  of  matrix. 

MTR  SIZ  -  matrix  size. 

printout : 

see  above. 
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FRTMBN  ** 


PRTMBN  prints  an  M  by  N  matrix  that  is  in  master  storage. 
calling  sequence: 

CALL  PRTMBN(N  MER,  HER  M,  MER  N) 

input : 

N  MTR  -  master  storage  location  of  matrix. 

MTR  M  -  number  of  rows  in  matrix. 

MER  N  -  number  of  columns  in  matrix. 

printout : 

see  above. 
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FRMARI  ** 


PRMARI  prints  the  magnitude,  angle  (either  in  degrees  or  radians), 
real,  and  imaginary  parts  of  a  complex  number. 

calling  sequence: 

CALL  PRMARI( IHEAD,  IANGLE,  DREAL,  DIMAG) 

input: 

IHEAD  -  flag  that  determines  whether  a  heading  should  be 
printed  (=0,  no  heading). 

IANGLE  -  flag  that  specifies  either  degrees  ( =5hdeg )  or 
radians. 

DREAL  -  real  part  of  input  number. 

DIMAG  -  imag  part  of  input  number. 

printout : 

see  above. 
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HEADP4  ** 


HSAD2U  prints  a  given  string  of  2k  Hollerith  characters.  The  first 
character  need  not  be  blank  as  the  output  format  skips  the  first  print 
position. 

calling  sequence : 

CALL  HEAD?L(I  HEAD) 

input : 

I  HEAD  -  an  array  containing  the  character  string  to  be 
printed. 

printout ; 

see  above. 


A  -  94 


PAGE  ** 


PAGE  initiates  a  page  eject  on  the  output  form. 
calling  sequence : 

CALL  PAGE 
printout : 

see  above. 
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STOLIN  ** 


STRLIN  prints  a  line  consisting  of  80  asterisks  on  the  output  form. 
cnlUng  sequence : 

CALL  STRLIN 
printout : 

see  above. 
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SKPLIN  ** 


SKPLIN  causes  a  desired  number  or  lines  to  be  skipped  on  the  output 

form. 

calling  sequence : 

CALL  SKPLIN ( NUMLIN ) 

INPUT: 

NUMLIN  -  number  or  lines  to  be  skipped. 

printout : 

see  above. 
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FRPTBL  ** 


PRPTBL  prints  the  interconnection  (or  port)  table. 
calling  sequence: 

CALL  PRPTBL 
printout : 

see  above. 
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PRCTLG  ** 


PRCTLG  prints  the  master  storage  area  cataxog. 
caxxing  sequence : 

CALL  PRCTLG 
printout ; 

see  above. 
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UNPACK  *** 


UNPACK  decodes  the  input  string  of  Hollerith  characters  (until  it 
reaches  a  $)  to  determine  such  items  as  sign  conventions,  accessible  port 
number,  network  name,  number  of  ports,  and  piece  number  if  applicable. 

calling  sequence : 

CALL  UNPACK(  IPOiRT,  NIPORT,  MCOUNT,  N,  ISIGNC,  ISXGND, 

JPORT,  NPORTS,  NETNAM,  IPC,  NCHAR ) 

INPUT: 

IPORT  -  array  containing  string  of  Hollerith  characters. 

NIPORT  -  maximum  number  of  words  in  array  IPORT. 

MCOUNT  -  counter  for  actual  number  of  words  used  in  IPORT. 

N  -  counter  for  number  or  bytes  m  1P0RT(I). 

output : 

ISIGNC  -  sign  of  C-matnx  entry  (dependent  variable). 

ISIGND  -  sign  oi  D-matnx  entry  (independent  variable). 

JPCSRT  -  accessible  port  number. 

NPORTS  -  network  port  number. 

NETNAM  -  network  name. 

IPC  -  piece  number  11  applicable. 

NCHAR  -  counter  ror  number  oi  characters  in  NETNAM. 

note:  tnis  routine  calls  the  FID  function  available  in  FORTRAN  V  on  the 
Univac  1106  and  1108.  For  other  systems  it  will  need  to  be  provided. 
Although  more  general  on  the  Univac  systems,  the  FLD  function  is  used  here 
only  to  extract  the  left  most  byte  from  a  word. 

calling  sequence : 

ICRAR( IF,  IL,  IWCRD) 


input/output 


input: 


A  -  100 


IF 


EL 

IWCRD 


-  bit  in  IWQRD  to  begin  extracting  (==0  for  this 
case) . 

-  bit  in  IWQRD  tc  end  extracting  ( =nuaber  of  bits  per 
byte  for  this  case. 

-  word  from  which  bits  are  to  be  extracted. 
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BLKDTA  * 


BLKDTA  is  a  BLOCK  DATA  routine  used  only  for  initializing  certain 
auantities  in  COMMON. 

Note:  COMMON  /  CPT  WRD  /  N  BYT  P  W,  N  B  P  BYT  must  be  set  properly 
here  for  each  particular  computer . 

(See  Page  A-^) 
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TIME  ** 


TIME  uses  the  subroutine  LKCIKS  (net  provided  in  this  package)  to 
determine  the  present  clock  time  and  prints  the  elapsed  time  from  the  last 
call  to  TIME,  as  well  as  the  total  time  since  CALL  IKCLKS  (another  routine 
not  provided  which  initializes  time  in  LKCLKS  to  0). 

calling  sequence : 

CALL  TIME 

note:  LKCLKS  returns  the  current  clock  time  in  seconds. 
calling  sequence : 

CALL  LKCLKS(T) 


INPUT: 


T  -  clock  time  (seconds) 


printout : 

see  above. 
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Z1T0Z2  * 


Z1T0Z2  finds  the  impedance  (Z?)  seen  at  output  port  given  the  A- form 
of  the  network  matrix  and  the  impedance  (Zl)  at  input. 

calling  sequence; 

CALL  Z1T0Z?( NAME,  NPC  Z,  DZ1R,DZ1I,  DZ2R,DZ?l) 

INPUT; 

NAME  -  network  name(  unpacked ) . 

NPCZ  -  piece  number  (0  for  entire  section). 

DZIR  -  real  part  of  input  impedance. 

DZ1I  -  imag  part  of  input  impedance. 

output ; 

DZ2R  -  real  part  of  output  impedance. 

DZ2I  -  imag  part  of  output  impedance. 

where : 


Z?  =  (A22*Zl-A12)/(-A?l*Zl+All) 


DCPCT?  ** 


DCPCT?  computes  the  force  and  velocity  on  the  left  of  a  piece  (FI, VI) 
given  the  force  and  velocity  on  the  right  ( Fr,Vr )  and  the  Px2  A  matrix. 

calling  sequence ; 

CALL  DCPCTp( NAME,  N  PC  FV,  DPRR,DFRI,  DVRR,DVRI, 

DFLR,DFLI,  DVLR,DVLI) 

INPOT: 

NAME  -  network  name( unpacked ) . 

N  PC  FV  -  piece  number  (0  for  entire  section). 

DFRR  -  real  part  of  Fr. 

DFRI  -  imag  part  of  Fr. 

DVRR  -  real  part  of  Vr. 

DVR I  -  imag  part  of  Vr. 

output; 

DFLR  -  real  part  of  FI. 

DFLI  -  imag  part  of  FI. 

DVLR  -  real  part  of  VI. 

DVLI  -  imag  part  of  VI. 

where : 


FI  =  All*Fr  +  A12*Vr 


VI  =  A2l*Fr  +  A22*Vr 
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FNDKRC  ** 


FNDKRC  computes  the  right  handed  decompaction  matrix  for  a  3-3  mode 
ceramic  stack  (network  type  11).  (See  Appendix  D.  ) 

calling  sequence : 

CALL  FNDKRC ( IZ  C3S,  IN  ZC3,  IKR3C) 

input : 

IZ  C3S  -  master  storage  location  for  C3S  matrix. 

IN  ZC3  -  master  storage  location  for  single  piece  matrix. 

IKR3C  -  master  storage  location  for  storing  right  handed 

decompaction  matrix. 


SRCHBC  ** 


SRCHBC  performs  the  Search  process  as  described  in  Appendix  D. 
calling  sequence: 

CALL  SRCHBC(IA5S,  IBUC,  IC}S,  IR3C,  IMPS,  ISPS) 

input: 


IA5S 

-  master 

storage 

location 

for 

matrix 

A5S. 

IBUC 

-  master 

storage 

location 

for 

matrix 

BUC. 

IC}S 

-  master 

storage 

location 

for 

matrix 

C^S. 

IP3C 

-  master 

storage 

location 

for 

matrix  R3C. 

IMPS 

-  master 

storage 

location 

for 

matrix 

MPS. 

IS2S 

-  master 

storage 

location 

for 

matrix 

S?S. 
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MINCHK  *** 

MINCHK  compares  two  input  numbers  to  find  the  minimum.  If  a  new 
minimum  is  found,  then  a  new  identification  subscript  is  stored,  and  a 
flag  is  set. 

calling  sequence: 

CALL  MINCHK (DNEW,  DMIN,  INEW,  IMIN,  MINFND) 

INPUT: 

DNEW  -  new  value. 

input /output : 

IMIN  -  minimum. 

input : 

INEW  -  identification  subscript  corresponding  to  DNEW. 

input /output : 

IMIN  -  identification  subscript  corresponding  to  DMIN. 

output : 

MINFND  -  flag  to  indicate  that  new  minimum  was  found. 

logic : 


If  (DNEW  .LT.  DMIN)  then  DM3N=DNEW,  BfEN-ONEW,  MINFND=1 
otherwise  MINFND  =  0 
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MAXCHK  *** 


MAXCHK  compares  two  input  numbers  to  find  the  maximum.  If  a  new 
maximum  is  found,  then  a  new  identification  subscript  is  stored,  and  a 
flag  is  set. 

calling  sequence : 

CALL  MAXCHK (DNEW,  DMAX,  INEW,  IMAX,  MAXFND) 

input: 

DNEW  -  new  value. 

input/output : 

DMAX  -  maximum. 

input : 

INEW  -  identification  subscript  corresponding  to  DNEW. 

input /output : 

IMAX  -  identification  subscript  corresponding  to  IMAX. 

output : 

MAXFND  -  flag  to  indicate  that  new  maximum  was  found. 

logic : 


If  (DNEW  .GT.  DMAX)  then  DMAX=DNEW,  IMAX=INEW,  MAXFND=1 
otherwise  MAXFND  =  0 
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CMPZS  ** 


CMPZS  is  used  in  conjunction  with  the  SEARCH  procedure  (Appendix  D) 
to  compute  and  store  the  self  impedance.  This  routine  should  be  modified 
by  the  user  when  a  non-zero  self  impedance  is  required. 

calling  sequence: 

CALL  CMPZS 

printout: 

As  written  this  routine  prints  ZSGLF=0  only  once  in  a  run. 
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PREPZV  * 


PREPZV  reads  from  cards  or  card  images  on  magnetic  tape  a  heading 
card  followed  by  an  array  of  radiation  impedances.  The  array  is  scaled  by 
desired  scale  factors  and  written  onto  an  auxiliary  storage  unit  for  use 
later  in  the  calculations. 

calling  sequence: 

CALL  PR EPZV( NARRAY,  IFARAY,  INZVUN,  IOUTUN,  FSCALE,ZSC ALE, VSCALE) 

input : 

NARRAY  -  array  number. 

IFARAY  -  frequency  at  which  the  radiation  impedances  were 
generated. 

INZVUN  -  logical  input  unit  number. 

I017TUN  -  logical  output  unit  number  (not  equal  to  INZVUN ) . 

FSCALE  -  scale  factor  for  IFARAY. 

ZSCALE  -  scale  factor  for  radiation  impedances. 

VSCALE  -  scale  factor  for  velocities  associated  with  the 

impedances. 

heading  card  input  format: 

PORMAT(AL,A?,  IE,  15,  i6aL) 

columns  ?-  6  contain  ARRAY  (Hollerith  string) 

columns  7-10  contain  array  number 

columns  11-15  contain  number  of  impedances  to  read  in 
columns  16-79  contain  heading  information 
data  card  input  format: 

FCRMAT(?I5,  15,  515,  2E15.8,  2E15.8) 

SEE  LDZRVR  for  arguments 
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printout: 

prints  the  input  card  information. 


A  -  112 


FNDFZV  *** 


FNDFZV  reads  the  array  of  radiation  impedances  that  has  previously 
been  read  in  by  FPEPZV  for  the  purpose  of  checking  for  proper  array  number 
and  frequency. 

calling  sequence : 

CALL  FNDFZV ( IF  SWP) 


input : 


IF  SWP  -  not  presently  used  by  the  routine.  However  it  is 

the  value  of  the  frequency  (from  the  Search  process 
frequency  sweep)  in  the  event  that  the  user  wants  to 
change  the  array  of  radiation  impedances  as  a 
function  of  frequency. 
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LDZRVR  *** 


LDZRVR  reads  a  card  image  containing  radiation  impedance  information 
used  in  the  Search  process. 

calling  sequence : 

CALL  LDZRVR 

modified  card  image  format: 

FaRMAT(?lM5,5l3,  ?(ElU.7,E13.7)) 


11 

- 

row  number  of  element. 

J1 

- 

column  number  of  element. 

IF 

- 

frequency  for  radiation  impedance. 

ICODE 

- 

identification  code  supplied  by  user 

IB  ETAS 

- 

azimuth 

IGAMAS 

- 

depression  elevation. 

ZR 

- 

real  part  of  radiation  impedance. 

ZI 

- 

-traag  part  of  radiation  impedance. 

VR 

- 

real  part  of  velocity. 

VI 

imag  part  of  velocity. 
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SRMIDB  * 


SRMIDB  performs  the  single  reactor  midband  logic  as  described  in  the 
Search  process  of  Appendix  D. 

calling  sequence : 

CALL  SRMIDB( DFP ,  L  MX  MN,  L  EC  IC,  L  FLG  SP,  DLRM,  DLQK) 


input : 


DFP 


-  frequency 


L  MX  MN  -  flag  specifying  either  to  maximize  or  minimize. 
=  3HMAX  to  maximize 
ft  3HMAX  to  minimize 


L  EC  IC  -  specifies  whether  I  Zee  I  or  IZicl  is  to  be 
extremized. 

=  ?HEC  to  extremize  I Zee  I 
^  PHEC  to  extremize  IZicl 

L  FLG  SP  -  specifies  either  series  or  parallel  reactor. 
=  1HS  for  series 
ft  1HS  for  parallel 


DLBM 


-  series  resistance  of  the  reactor 


DLQM 


-  reactor  quality  factor 


printout : 

values  of  input  and  output  (controlled  by  COMMON  /  PRT  OFF  /  NO 
PRT). 
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LDAM?S  ** 


LDAM2S  computes  and  stores  the  A-form  of  the  M?S  network (  matching 
network)  used  in  the  Search  procedure  (Appendix  D). 

calling  sequence : 

CALL  LDAM2S 

printout : 


A  M?S  matrix  (controlled  by  COMMON  /  PRT  OFT  /  NO  PRT). 


CMPGAM  ** 


CMPGAM  is  used  in  the  Search  procedure  discussed  in  Appendix  D  to 
compute  (l)  voltage-velocity  control  impedance,  Zee,  (?)  current- velocity 
control  impedance,  Zic,  (3)  GAMMAec  =  Zself  +  Zic,  (U)  GAMMAic  =  Zself  + 
Zic  where  Zself  is  the  self  impedance  of  the  element. 

calling  sequence: 

CALL  CMPGAM( IN  A  A?S) 


input : 


IN  A  A?S  -  master  storage  location  of  the  A-form  matrix  for  the 
A?S  network. 


printout: 

A  S?R  matrix  and  GAMMA  ec,  GAMMA  ic  (controlled  by  COMMON  /  PRT  OFF  / 
NO  PRT). 
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FSRCH  * 

FSRCH  searches  out  performance  variables  for  a  transducer  over  a 
specified  frequency  range  with  an  input  array  of  of  radiation  impedances. 
The  method  used  is  discussed  in  Appendix  D. 

calling  sequence; 

CALL  FSRCH( IF  FRST,  IF  LAST,  IF  DEUT, 

I  OUT  OP,  I  UNIT  T,  IHLBLL,  IHLBL?) 

input : 

IF  FRST  -  initial  frequency. 

IF  LAST  -  maximum  frequency. 

IF  DELT  -  frequency  increment. 

I  OUT  OP  -  output  options  as  discussed  in  Appendix  D. 

I  UNIT  T  -  logical  unit  number  for  output  to  be  saved  (=0  if 

not  desired). 

THTBi i  _  first  word  of  a  label  for  output. 

IHLBL?  -  second  word  of  a  label  for  output. 

printout ; 

a  page  for  each  frequency  and  its  associated  performance  variables. 
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RLDMJ  ** 


RLBMJ  provides  the  functional  relationship  of  the  dependent  variable 
used  primarily  by  subroutine  FMFN.  The  routine  should  be  changed  for 
different  functional  relationships.  As  programmed  the  routine  finds  the 
admittance  Yc  of  the  ED  network  for  Run  80006  in  Appendix  E. 

calling  sequence: 

CALL  RLIMM3(DX,  DRL,  DIM,  DM3) 

input: 

DX  -  independent  variable. 

output: 

URL  -  real  part  of  dependent  variable. 

DIM  -  imag  part  of  dependent  variable. 

DM3  -  magnitude  of  dependent  variable. 

printout: 

above  values  (controlled  by  COMMON  /  C  FMN  PR  /  I  FMN  PR 

where  I  FMN  PR  =  0  for  no  print). 
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FMFN  * 


FMFN  computes  the  maximum  and/or  minimum  of  a  function(  evaluated 
RLIMMj)  given  an  initial  guess.  The  routine  utilizes  a  pole-z 
technique  which  is  further  refined  by  a  parabolic  fit. 

calling  sequence: 

CALL  FMFN(DFMOES,  DXM,  DYMM,  DXN,  DYNN) 

input: 

DFM3ES  -  initial  guess  for  independent  variable,  DXM. 

output : 

DXM  -  value  of  independent  variable  at  maximum  of 

dependent  variable. 

DYMM  -  relative  maximum  of  dependent  variable  in  region 

DFKJES. 

DXN  -  value  of  independent  variable  at  minimum  of 

dependent  variable. 

DYNN  -  relative  minimum  of  dependent  variable  in  region 

DFM3ES. 

additional  input  options  can  be  set  via 

COMMON  /  CMN  FMN  /  MIN  SKP,  MAX  SKP,  N1  POL  0 
MIN  SKP  -  flag 

=  1  will  not  find  minimum  of  the  function. 

/  1  will  find  minimum  of  the  function. 

MAX  SKP  -  flag 

=  1  will  not  find  maximum  of  the  function. 

/  1  will  find  maximum  of  the  function. 

N1  POL  0  -  flag  controlling  pole-zero  technique. 

=  1  apply  pole-zero  once 
/  1  apply  pole  zero  twice 

note:  all  three  flags  are  set  to  0  as  a  default  condition  in  BLK  DTA. 
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POLEO  *** 


POLEO  is  used  to  find  the  frequency  at  which  a  relative  maximum  and  a 
relative  minimum  of  the  complex  admittance  function  occurs  given  three 
frequencies  for  the  function  to  be  evaluated.  This  routine  has  been 
designed  to  apply  particularily  to  piezoelectric  electromechanical 
systems.  The  admittance,  Yc,  of  such  systems  with  widely  separated 
resonant  behavior  can  be  approximated  at  frequencies  well  below  the  second 
resonance  by  Y(f)  =  f(f-b)/(f-c).  (note:  for  magnetostrictive  systems  a 
similar  type  of  treatment  will  apply,  but  Y  must  denote  complex  impedance 
instead  of  complex  admittance.  ) 

calling  sequence: 

CALL  PCLEO(  DX1 ,  DX2,  DX3,  DXM,  DXK) 


input : 

DX1  -  first  value  of  the  independent  variable. 

DX?  -  second  value  of  the  independent  variable. 

DX3  -  third  value  of  the  independent  variable. 

output: 

DXM  -  value  of  the  independent  variable  at  which  the 

dependent  variable  is  a  maximum. 

DXN  -  value  of  the  independent  variable  at  which  the 

dependent  variable  is  a  minimum. 


method : 

The  technique  was  developed  from  the  idea  that  the  frequency  poles 
and  zeros  determine  the  admittance  as  a  function  of  frequency. 

Y(f)  =  (A(f-al)(f-a2)(f-a3)...)/((f-bl)(f-b2)(f-b3)...) 
where  al,a2,a3, . . . ,bl,b2, . . .  are  complex  constants. 

For  this  case  Y(f)=0  if  f=0.  Therefore,  al  is  set  to  0.  It  may  also 
be  assumed  that  only  a  small  portion  of  the  complex  frequency  affects  the 
function  so  that  only  one  pole  and  one  zero  need  be  considered,  i.e. 

Y(f)  =  f(f-b)/(B(f-c)) 
or 

Y(f)/f  =  B  +  B(c-b)/(f-c) 
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Y(f)/f  =  B  +  D/(f-c). 

Given  the  value  of  the  admittance  function,  Y(f ),  at  three 
frequencies  (fl,f?,f3)  allows  a  solution  for  the  constants  B,  D,  and  c. 
The  solutions  are 

B  =  -D/(fl-c)  +  Y1 

D  =  (  Yl-Y?)/( 1/ ( f 1-c )-l/( f ?-c ) ) 

c  =  (fl(f?-f3)Yl  +  f?(f5-fl)Y?  +  f3(fi-f2)Y3)  / 

((f?-f3)Yl  +  (f3-fl)Y2  +  (fl-f?)Y3) 

b  =  c  -  D/B 

Therefore,  the  value  of  b  gives  the  frequency  at  which  a  minimum 
occurs,  and  c  gives  the  frequency  at  which  the  maximum  occurs . 
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PRBFIT  *** 


PRBFIT  performs  a  parabolic  fit  to  the  dependent  variable  at  the 
input  value  of  the  independent  variable  to  determine  either  a  maximum  or 
minimum  point.  The  dependent  variable  functional  relationship  is 
described  through  RLBM3. 

calling  sequence: 

CALL  PRBFIT(DX,  DYM3) 

input /out put: 

DX  -  independent  variable  input  value  which  will  be 

modified  to  give  a  value  closer  to  the  actual 
maximum  or  minimum. 

output : 

DYMS  -  maximum  or  minimum  value  of  the  dependent  variable 

corresponding  to  DX. 

method: 

Assume  that  a  parabolic  fit  to  the  function  y  =  f(x)  is  desired  at 
the  point  xm.  Define 

eps  a  xm  *  l.E-4 


xl  =  xm  -  eps 


x2  =  xm 

x3  =  xm  +  eps 

and  determine  f(xl),  f(x2),  f(x3). 

Calculate 

eps  prime  =  (eps/2)( lf(xl) i  -  If(x3)l)  / 

(lf(xl)l  -  2 1 f ( x2 ) I  +  If(x3)l) 

and 

xm  *  x2  +  eps  prime 

If  llf(xm)l  -  I f ( x2 ) I  I  /  lf(xm)l  Is  less  than  the  error  tolerance,  then 
stop.  Otherwise  do  the  following  replacements  and  cycle  again: 
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eps  =  eps  prime 
xl  =  x? 
x?  =  xm 
x3  =  xm  +  eps. 
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FNDPRM  * 


FNDPRM  iterates  on  the  given  ceramic  parameters  until  agreement  with 
experimental  results ( i .  e . , FM,  FN,  YM,  YN,  CLO,  CLOM)  is  obtained  within  a 
certain  tolerance.  (See  Appendix  E.  ) 

calling  sequence: 

CALL  FNDFRM(C  NAME,  FM,  FN,  YM,  YN,  CLO,  CLOM,  FLO,  MAX) 


input : 


C  NAME  -  name  of  the  ceramic  section. 


FM 

FN 

YM 

YN 

CLO 

CLOM 

FLO 

MAX 


printout : 


-  frequency  where  admittance  is  a  maximum. 

-  frequency  where  admittance  is  a  minimum. 

-  maximum  admittance. 

-  minimum  admittance. 

-  the  electrical  input  parallel  capacitance  at  the 
frequency  FLO. 

-  the  dissipation  factor  at  frequency  FLO. 

-  frequency  which  is  below  FM  and  near  the  lower 
frequency  limit  of  the  actual  transducer  operating 
frequency. 

-  max  number  of  iterations. 


see  writeup  for  subroutine  PRTPRM, 


FNDFDP  *** 


FTOFDP  computes  the  values  to  compare  with  experimental  data  in  the 
Find  Parameters  procedure  described  in  Appendix  E. 

calling  sequence; 

CALL  FNDFDP 
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LDBCOL  *** 


LDBCOL  loads  a  column  of  the  Find  Parameters  matrix.  (See  Appendix 

E.  ) 

calling  sequence ; 

CALL  LDBCOL(J  COL,  N  CAP) 

input: 

J  COL  -  column  number  to  be  loaded. 

N  CAP  -  total  number  of  rows  in  matrix. 
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PRTPRM  *** 


PRTPRM  prints  the  output  from  the  Find  Parameters  procedure,  (see 
Appendix  E.  ) 

calling  sequence: 

CALL  PRTPRM 

printout : 

prints  at  each  iteration  the  experimental  values,  the  new  values 
vising  the  newly  computed  parameters,  and  the  relative  errors. 
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APPENDIX  B 


IMPLEMENTED  NETWORK  TYPES 


Table  1  summarizes  the  networks  presently  implemented  in  the  program.  For  each 
network  the  input  data,  sign  conventions,  and  equations  programmed  are  included  in  the 
descriptions  which  follow: 
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TABLE  1 .  Summary  of  Network  Types, 


(*  indicates  frequency  dependence) 


Type  Description  Input  Data 

1  *  Straight  Piece  NAME.RHO.A.ID.OD.L.Cr.Ci.Cm 

2  *  Conical  Piece  NAME.RHO.A1  .ID1  ,ODl  .A2.ID2. 

OD2,L.Cr.Ci,Cm 

4  Ideal  Transformer  NAME. TURNS  RATIO 

5  Arbitrary  n  by  n  Matrix  NAME.ZIJr.ZIJi 

7  *  Rational  Polynomial  in  Frequency  NAME.F ACTOR, IFORM.NSYMf  G. 

coefficients 


8  Voltage  to  Voltage  Converter 

Q  Voltage  to  Current  Converter 

1  1  *  3-3  Mode  Ceramic  Stack 

14  *  3-1  Mode  Ceramic  Thin  Bar 

15  *  3-1  Mode  Ceramic  Tube 


NAME. ZI.ZO. ALPHA 

NAME.ZI.ZO.BETA 

NAME.N.RHO.A.ID.OD.L. 

E33Tr.E33Ti,E33Tm, 

G33r.G33i.G33m. 

S33Dr,S33Di,S33Dm 

NAME.RHO.T.W.L. 

E33Tr,E33Ti.E33Tm, 

G31r.G31i.G31m, 

SI  1  Er.Sl  1  Ei.Sl  1  Em 

NAME.RHO.ID.OD.L, 

E33Tr,E33Ti,E33Tm, 

G31r.G31i.G31m, 

SI  1  Er.Sl  1  Ei.Sl  1  Em 
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*  Linear  Transformer 


NAME.RP.LP.RO.LO.RS.LS.PEII 


21  *  Series  Impedance 

22  *  Parallel  Impedance 

24  *  Series  Admittance 

25  *  Parallel  Admittance 

26  *  Linear  Model  of  a  Differential 

Input  Operational  Amplifier 


NAME.RS.  LS.CS 
NAME, RS, LS.CS 
NAME.RP.LP.CP 
NAME.RP.LP.CP 

NAME.R1  .Cl  .R2.C2.RD.CD,R3,L3, 
AO.  OMEGA  1 ,0MEGA2.0MEGA3 
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TYPE  1 :  STRAIGHT  PIECE  (non-ceramic) 


Input  Parameters: 

") 

A  =  Area  (m“) 

L  =  Length  (m) 

c  =  Sound  Velocity  =  cr  (l+jcm)  (m/sec) 
p  =  Density  (kg/m  ) 

OD  =  Outer  Diameter  (m)  1 


when  applicable 


ID  =  Inner  Diameter  (m) 
f  =  Frequency  (Hz) 
cc  =  2rrf 


(1)  (2) 


L 


Computed  Parameters: 
k  =  co/c 
A 

A 


11  =  A-n  =  cos  (kL) 

Zll 

=  -jpcAcot  (kL) 

.j  n  -  jpcAsin  (kL) 

Zl2 

=  -Z-q  =  jpcAcsc  (kL) 

2]  =  j(sin(kL))/(pcA) 

Z-»2 

=  jpcAcot  (kL) 
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Note  that 
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TYPE  2:  CONICAL  PIECE  (non-ceramic) 


Input  Parameters: 

Aj  =  Area  at  ( 1  )  (m~) 

A 

A->  =  Area  at  (  2)  (m~ ) 


L  =  Length  (nt) 

c  =  Sound  Velocity  =  cr  ( 1+jcm)  (m/sec) 


p  =  Density  (kg/m^) 

ODj  =  Outer  Diameter  at  (1 )  <m) 
ODt  =  Outer  Diameter  at  (2)  <m  ) 
ID]  =  Inner  Diameter  at  ( 1  )  (m) 
ID-i  =  Inner  Diameter  at  (2)  (m) 
f  =  Frequency  (Hz) 
co  =  27rf 


when  applicable 


(1) 


(2) 


vi  0Dr 


Computed  Parameters: 


k  =  co/c 

A]j  =  x/Aj/At  [cos(kL)  -  sin  (kL)/(kL)]  +  sin  (kL)/kL 

A]  i  =  jpc[  1  /( kL )  ( sin( k L )/( k£ )  +  cos(kL))  (  -  CAt  )"  +  v^l  A  2  SU1  (^L)] 

A21  =  j  sin  (kL )/(pc ^/A ]  At) 

An  =  /4,/A]  [cos  (kL)  -  sin  (k L )/( kL ) ]  +  sin  (kL)./(kL) 
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Z]  ]  =  - jpcAj  [cot  (kL)  -  1  /(kL)]  -jpcAT^n/CkL)] 
Z  i  2  =  -  Z->  j  =  jpc  \/Aj  A  -i  esc  ( k  L ) 


Ztt  =  jpcA2  [cot  (kL)  -  1  /(kL)]  +jpcv/A]  A->  [  1  /(kL)] 
Note  that 


"pr 

N 

1 _ 

zi  2 

1 

f2 

Cl 

N 

_ 1 

Z22 

1 

< 

tj 

i _ 

1 

T1 

< 

1 _ 

A1  2 

1 - 

r  1 

U, 

1 _ 

1 - 

< 

1 _ 

_a21 

a22 

V2J 

TYPE 

2 


V.  F. 
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TYPE  4:  IDEAL  TRANSFORMER 


Input  Parameters: 

N  =  Secondary /Primary  Turns  Ratio 


Computed  Parameters: 

Aj  j  =  1/N 

Aj  n  =  0 

A21  =  0 
An  =  N 

If  N=1 .  [A]  =  Unit  Matrix, 
where 


V 

A1  1 

> 

to 

_ 1 

v 

A21 

A~>n 

h 

NOTE:  THE  Z-FORM  DOES  NOT  EXIST 
THE  Y-FORM  DOES  NOT  EXIST. 
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TYPE  5 :  Arbitrary  n  by  n  Matrix 

Input  Parameters:  List  of  complex  numbers  for  elements  of  an  n  by  n  matrix. 
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TYPE  7:  RATIONAL  POLYNOMIAL  IN  FREQUENCY 
Input  Parameters: 

FACTOR  =  scale  factor  for  input  coefficients 
IFORM  =  flag  for  form  of  input  matrix 
NSYMFG  =  flag  for  symmetry'  of  input  matrix 

Cnjj  =  complex  coefficients  for  numerator  polynomial  (n=0,  .  .  .  .  N) 
dmjj  =  complex  coefficients  for  denominator  polynomial  (m=0.  1  .  .  . 
N  =  order  of  numerator  polynomial 
M  =  order  of  denominator  polynomial 
f  =  frequency  (Hz) 

where: 

-  i.J  =  1.2 . K 

K  =  the  number  of  network  ports 

Computed  Parameters: 

N 

y  c  ■  *  fn 

n=0 

2,1  =  A - 

1  dn,y*fra 

m=0 

Note:  For  definitions  of  IFORM  and  NSYMFG,  see  Appendix  A 
subroutine  writeup  for  LD  DT  7. 


M) 
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TYPE  S :  VOLTAGE  TO  VOLTAGE  CONVERTER  (can  be  used  us  Amplifier  mm 


Input  Parameters: 

Zj  =  Input  impedance  of  amplifier  (ohms) 

Zq  =  Output  impedance  (ohms) 
a  -  Voltage  gain 

Note:  All  three  parameters  are  complex  and  a  function  of  frequency  . 


Computed  Parameters: 


A11 

=  1/a 

N 

II 

_N 

Aj  -> 

=  Z  Q/a 

Zj  2  =  o 

At  j 

=  l/(aZj) 

Z-)  j  —  aZj 

A22 

=  ZQ/(aZi) 

z22  =  -ZG 

Note  that 


V 

> 

> 

i  j 

^ _ J 

- 1 

tn 

Ji. 

]  Ann 

L'-J 

hi 

Z11  Z12 

V 

i — 

tn 

Zoj  Ztt 

It 

Y  -  Z"1 
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TYPE  9:  VOLTAGE  TO  CURRENT  CONVERTER 
(can  be  used  as  Norton  Amplifier  Model) 


Input  Parameters: 

Zj  =  Input  impedance  of  amplifier  (ohms) 

Zq  =  Output  impedance  (ohms) 
ft  =  Voltage  to  current  gain  (transconductance.  Mhos) 
All  three  parameters  are  complex  and  a  function  of  frequency. 


Computed  Parameters: 

Aj  j  =  1/(/3Zq) 

Aj  2  =  1/(3 

At  j  =  1  / (jSZj  Zq) 
A22  =  l/(j3Zj) 
Note  that 


"El" 

A1 1 

> 

t  j 

_ i 

Jl„ 

_A21 

A“J 

'El" 

i — 

Cl 

nT 

e2 

_Z21 

Zn 

Y  =  Z-1 


Z1  2  =  0 
Z2]  =  0Zj  Zq 
Z22  =  -Zq 


I  A  I  f  1 
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TYPE  11:  3-3  MODE  CERAMIC  STACK 


Input  Parameters: 

p  -  Density  (kg/m  ) 

A  =  Area  of  poled  surface 
L  =  Piece  length  (m) 

N  =  Number  of  pieces 
OD  =  Outer  diameter 
ID  =  Inner  diameter 
Ej3  =  E33TR  (1  - jE33TM) 
g33  =  G33R  (1  -jG33M) 

s?3  =  S33DR  (1  - jS33DM) 
f  =  Frequency  (Hz) 

<jO  =  27rf 


(m“) 


when  applicable 


Ceramic  stack  parameters 
in  RMKSAa  Units 


Computed  Parameters: 

c  =  /(ps33) 

k  =  cc/c 

cq  =  Permittivity  of  vacuum  (8.8542  X  10_^“) 

P33  =  ^§33  ^“/s33  +  1  ^eOE33^ 

C0  =  A/O^f) 

$  =  g33C0/s33 
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Zq  =  <L/(jcoC0) 

z0  =  zo'N 

=  jpcAcot  (kL)  +  4>zq  =  -Zp 

c  c 

Zp  =  jpcAcsc  (kL)  +  4>Zq  = -Z-, j 

a  =  cosh-*  (Z^/Zp) 

Zq  =  -Zp  sinh  (a) 

K33  =  g33 

Zp  =  ZQCOth  (No)  +  4>Zq 
Zp  =  -Zposch  (Na)  -  4>Zq 

Q 

Z]  3  -  z0 
Z]]  =  -Zp 

Z22  =  -Z1  1 
Z23  =  Z1 3 
Z3 1  =  Z13 
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TYPE  14:  3-1  MODE  CERAMIC  THIN  BAR 


Input  Parameters: 

p  =  density  (kg/m:' ) 
t  =  thickness  (m) 
w  =  foil  width  (m ) 

L  =  length  (m) 

eJ3  =  E33TR  (1  -  jE33TM ) 
g31  =  G31R  (1  -jG31M) 

sfj  =  SI  1ER  (1  - jSl  1  EM) 

f  =  frequency  (Hz) 
cu  =  27rf 


COMPUTED  PARAMETERS: 


k  -  co/c 

e0  =  8.8542  X  1CT12 

D31  =  E33  X  eO  X  £31 
(k3])2  =  g3j  X  D3]/sEj 

Amech  *  X  w 
Aelec  =  E  X  w 
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mass  -  p  X  Amech  X  L 


C0  =  (-^Y2)  X  £o  X  E33  X  "  -  K31 1 

•»  =  \-lcx  /  I->3 1  /(  'Vi  XL) 
z33  =  l/(j^C0) 

z0  =  $  x  Z33 

zfj  -  -jpC  X  Amech/tan  (k  X  L)  -  -Z^ 
zf2  =  jpC  X  Amech/sin  (k  X  L)  =  -Z§j 
Z  ]  1  =  Z^j  +  ^zq  =  _Z22 

Zj2  =  z]2  ~  ®z0  ~  ~Z21 
Z1 3  =  zO  =  Z31  =  Z23 

Z32  =  “Z1 3 
Note  that 
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TYPE  15:  3-1  MODE  CERAMIC  TUBE 


Input  Parameters: 

p  =  density  (kg/m  .) 

ID  =  inner  diameter  (m) 
OD  =  outer  diameter  (m) 

L  =  length  (m)  (all  foiled) 

eJ3  =  E33TR  (1  -jE33TM) 
g3]  =  G31R  (1  -jG31M) 

sE  =  S11ER  (1  -  jSl  1  EM) 

f  =  frequency 
cu  =  27rf 


COMPUTED  PARAMETERS: 

1 

c  = - 

p  X  Sj  i 

k  =  to/c 

e0  =  8.8542  X  1(T12 

D31  =  E33  X  eO  X  S31 
(K31)-  =  g31  X  D31/sEj 

Amech  =  (OD  +  ID)  (0D  "  ID)  X  */4 
Aelec  =  (0D  +  ID)  x  L  X  tt/2 
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mass  =  p  X  X  L 


Cq  =  27rL/Cn  (OD/ID)  X  eQ  X  Ej3  X  (1  -K3]  ) 

*  =  AelecXD3l/'(sn  X  L) 

Z33  =  l/jojC0) 
zO  =  *  X  Z33 

Z^j  -  jpc  X  Amech/tan  (k  X  L)  =  -Z22 
Z^-,  -  jpc  X  Amech/sin  (k  X  L) =  -Z?, 

Z1 1  =  zn  +  <l>zO  =  _Z22 

Z1 2  =  z^2  ~  $zO  =  ~Z21 
Z1 3  =  zO  =  Z31  =  Z23 

Z32  =  -Z1 3 
Note  that 


3 
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TYPE  20:  LINEAR  TRANSFORMER 


Input  Parameters: 

4>  =  Turns  ratio;  secondary /primary 
Rp.Rs  =  Copper  losses  in  primary  and  secondary  (ohms) 

Rq  =  Core  loss  (ohms) 

Lp,L§  =  Winding  inductances  in  primary  and  secondary  (H) 
Lq  =  Magnetizing  inductance  CH) 

Zp  =  Primary  impedance  -  Rp +jcoLp  (ohms) 

Yn  =  Shunt  Admittance  =— —  + — - -  (mhos) 

0  Ro  j"L0 

Zg  =  Secondary  impedance  =  Rg +jcoLg  (ohms) 

f  =  Frequency  (Hz) 
co  =  27rf 
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9  + 


Computed  Parameters: 


Zp  —  Rp  +  jcoLp 

Zq  =  1  /Y0 

ZSAP  "  Zs/<*’2 

Zs  =  Rs+jujLs 

zn  -  zP  +  z0 

Aj  1  -  [1  +  ZpY(~)]  /$> 

Zj  7  =  -4>Zq 

A12  =  ^[zSap  +  zp  +  zsapzpyo] 

Z21  =  12 

A2i  =  Yq/45 

Zto  =  -[4>-ZQ  +  Zg] 

A-»o  =  <£[1  +  ZSAP  Y0^ 
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TYPE  21 :  SERIES  IMPEDANCE 


Input  Parameters: 

R<5  =  Resistance  (ohms) 
Lg  =  Inductance  (H) 

C§  =  Capacitance  (F) 
f  =  Frequency  (Hz) 
cc  =  27rf 


(2) 


o 


Computed  Parameters: 


Let 


ZS  -  RS+J^LS 


and 


in 

Cg  =  10  for  no  capacitor, 
Lg  =  0  for  no  inductor. 


Rg  =  0  for  no  resistor. 

Then 

An  "  1- 

A1 2  ZS- 
A21  = 

and 

An  —  1  , 

where 


X 

A1  1 

> 

i  j 

_ i 

A. 

ri 

< 

_ i 

A-)oJ 

Eo 


IA 


NOTE:  THE  Z-FORM  DOES  NOT  EXIST 
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TYPE  22:  PARALLEL  IMPEDANCE 


In  put  Parameters: 


f 


CJ 


Resistance  (ohms) 
Inductance  (H) 
Capacitance  (F) 
Frequency  (Hz) 
27rf 


(2) 

— o 

+ 


-o 


Computed  Parameters: 


Let 

Cg  =  1 0^®  for  no  capacitor, 
Lg  =  0  for  no  inductor, 


and 

Rg  =  0  for  no  resistor. 

Then 


A11 

=  1, 

Z11  =  ZS’ 

A1  2 

=  o, 

Z12  =  _ZS 

A21 

=  }/Z  S’ 

Z21  =  ZS’ 

A^2 

=  1, 

Z22  “  _ZS 

where 


V 

rAn 

> 

to 

_ 1 

i 

tn 

to 

_ 1 

n 

<  j 

Ann 

hJ 

P'l 

'zn 

Z1  2 

v 

1 

r  4 

Jij 

_Z21 

1 

ro 

n 

N 

In 

— 

I  A  I  =  1 
IZI  -  0 


NOTE:  THE  Y-FORM  DOES  NOT  EXIST 
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TYPE  24:  SERIES  ADMITTANCE 


Input  Parameters: 

Rp  =  Resistance  (ohms) 
Lp  =  Inductance  (H) 

Cp  =  Capacitance  (F) 
f  =  Frequency  (Hz) 
co  =  2nf 


Computed  Parameters: 


Let 


Rp  =  10^0  for  no  resistor, 
Lp  =  1 0^®  for  no  inductor, 

and 

Cp  =  0  for  no  capacitor. 

Then, 


NOTE:  Z-FORM  DOES  NOT  EXIST 


I A  I  =  1 
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TYPE  25:  PARALLEL  ADMITTANCE 


Input  Parameters: 

Rp  =  Resistance  (ohms) 
Lp  =  Inductance  (H) 

Cp  =  Capacitance  (F) 
f  =  frequency  (Hz) 
cu  =  27rf 


Computed  Parameters: 


YP  = 

-  +j/cuCp-  ] 

Rp  y  cuLp  J 

Let 

Rp  = 

1 030  for  no  resistor. 

Lp  = 

1030  for  no  inductor. 

and 

n 

*0 

n 

0  for  no  capacitor. 

Then, 

A11 

-  1, 

A1  2 

-  0, 

A21 

=  YP< 

and 

A-n 

=  1. 

where 

11 

A1  2 

21 

a22 

NOTE:  Y-FORM  DOES  NOT  EXIST 
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TYPE  26:  LINEAR  MODEL  OF  A  DIFFERENTIAL  INPUT  OPERATIONAL  AMPLIFIER 


Input  Parameters: 

Rj  =  Input  resistance  (ohms) 

Cj  =  Input  capacitance  (F) 

Ro  =  Input  resistance  (ohms) 

CA  =  Input  capacitance  (F) 

Rq  =  Differential  input  resistance  (ohms) 

Cp  =  Differential  input  capacitance  (F) 

R3  =  Output  resistance  (ohms) 

L3  =  Output  inductance  (H) 

Aq  =  dc  open  loop  gain 

CO]  ,cot.co3  =  First  three  pole  frequencies  of  the  open  loop  transfer  function 
(radians/sec) 

f  =  frequency  (Hz) 
co  =  27rf 
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Computed  Parameters: 


Yj  =  1/Rj  +  jcoCj 
Y  =  1 /R-i  +  jcuCo 
Y3  =  1  /( R3  +  jcuL^) 

Yd  -  1  /Rd  +  jt^CD 


A  -  Aq/( ( 1  +jcu/c<j|)(l  +jcj/cj2)(l  +jto/cj3)) 
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APPENDIX  C 


THE  CONTROL  ROUTINE  USING  A 
TRANSDUCER  DESIGN  EXAMPLE 

This  appendix  presents  a  discussion  of  the  operations  that  are  written  tor  a  particu¬ 
lar  transducer  design  control  routine.  In  this  case,  the  specific  problem  dictates  part  ot  the 
routine;  however,  certain  operations  must  be  included  in  all  control  routines. 


The  first  step  in  writing  a  control  routine  is  to  define  the  amount  ot  storage  to  be 
used  in  the  master  storage  area  or  blank  common.  The  common  block,  DTA  TOT,  contains 
the  variable,  MX  DATA,  which  is  to  be  set  by  the  user  to  the  maximum  number  ot  storage 
cells  to  be  used  in  blank  common  (i.e.,  the  dimension  size  ot  the  blank  common  array , 
DATA). 


Secondly,  the  user  must  supply  the  port  table  description.  The  port  table  describes 
to  the  program  the  actual  interconnection  topology,  series  or  parallel  connection,  and  sign 
conventions.  This  operation  is  accomplished  through  subroutine  LD  PORT,  which  is  called 
once  for  each  accessible  port  of  the  network  (Note:  It  is  possible  to  define  several  networks 
at  this  time  for  use  later  in  the  program). 


Following  the  description  of  the  port  table,  each  basic  network  that  is  referred  to  in 
the  table  must  be  described  through  a  network  description  routine.  A  library  of  network 
description  routines  is  available  to  include  various  types  of  networks,  and  may  be  easily 
expanded  to  handle  any  special  types  desired  by  the  user.  The  purpose  of  these  routines  is 
to  supply  the  data  required  to  calculate  the  impedance  matrix  tor  the  network  (e.g.,  length, 
area,  velocity  of  sound,  etc.). 
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The  control  routine's  major  function  is  to  perform  the  necessary  calculations  (or,  as 
in  most  cases,  the  necessary  calls  to  subroutines  that  perform  the  calculations)  for  the  solu¬ 
tion  of  the  user's  specific  problem.  Since  the  specific  problems  vary  widely,  a  general  struc¬ 
ture  for  this  portion  cannot  be  provided.  This  approach  offers  the  user  maximum  generality 
and  flexibility  for  the  solution  of  his  problems. 

The  user  must  control  through  the  routine  the  printed  output  of  the  problem.  Sev¬ 
eral  print  flags  can  be  turned  on  by  the  routine  to  provide  some  printed  output;  however, 
the  user  will  have  to  supply  some  code  to  output  the  desired  information. 

The  control  routine’s  basic  functions  are  illustrated  by  a  routine  written  to  investi¬ 
gate  the  changes  in  transducer  performance  caused  by  varying  the  transducer's  tail  length 
and  ceramic  area.  Figure  C.l  shows  a  diagram  of  the  transducer.  Figure  C.2  shows  the 
listing  and  associated  printed  output  for  the  control  routine.  The  problem  in  this  case  is  to 
minimize  |Zecl.  the  voltage-velocity  control  impedance  of  the  radiating  face  (or,  equivalently, 
flee  =  ^ec  +  ^selfl  where  ^self  =  constant),  as  a  function  of  both  the  number  of  ceramic 
pieces  and  the  length  of  each  ceramic  piece.  Lines  8-22  of  figure  C.2  are  the  common 
blocks  needed  for  this  run.  Line  17  is  a  common  block  used  only  for  this  type  of  design 
procedure  and  is  used  primarily  for  transfer  of  information  to  RLIMMG,  a  subroutine  that 
determines  the  real  and  imaginary  parts  as  well  as  the  magnitude  of  a  complex  function 
evaluated  at  a  particular  point.  In  this  case,  the  function  to  be  minimized  is  programmed 
into  the  subroutine. 

Lines  27-33  contain  pertinent  data  such  as  the  names  of  networks  to  be  involved  in 
total  length  calculations,  ceramic  areas,  and  tail  lengths. 

Lines  46-59  define  to  the  routine  the  interconnection  table  and  sign  conventions 
utilized  by  the  user.  Each  of  the  basic  networks  that  appears  in  the  interconnection  table  is 
then  described  through  particular  network  type  calculations  in  lines  61-81. 
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Three  major  loops  begin  on  lines  108  (tail  lengths),  1  14  (ceramic  areas),  and  1  26 
(number  of  ceramic  pieces).  The  computations  that  are  a  function  of  these  variables  must 
be  defined  within  these  loops.  The  problem  is  to  determine  the  number  of  pieces  required 
to  minimize  I7  |  and.  using  that  number  of  ceramic  pieces,  to  refine  the  estimate  of  the 
minimum  by  varying  the  ceramic  piece  length.  The  refined  estimate  is  computed  by  both  a 
pole-zero  technique  and  a  parabolic  fit. 

The  subroutine  RLIMMG  is  presented  on  pages  8,  9  and  10  of  the  computer  print¬ 
out.  Computer  routine  loading  information  appears  on  pages  1  1  through  22,  while  pages 
through  30  are  given  to  computer  run  output. 
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I  J  . 


Figure  C.  2 . 

Listing  and  Output  for  Sample  Control  Routine  for 
Transducer  Design.  This  routine  investigates  the  effects 
on  the  performance  of  a  transducer  element  of  varying 
the  tail  length  and  ceramic  area. 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued). 


C-29 


iEALAJCEK  Kill.  atlOJa  PON  TP-229  APPEML-U  C  DATE  101972  PAGE 


■S’ 

Cu 


a  3  o 
0  3  0 
3C3EO 


d 

b  s 

TO  o  © 
<r  cj  <  c 

tL  05  O 

or  ^  *-*  a 

<t  r-r  « 

•  u 


O  _J  o 

0-3  0 

Q  o  2  o 

•-«  o  • 

=3° 


r-  e>  o 

«  t—  «.  o 
U  <V  2.  O 
X  rH  O 

4  H  • 

-  <_> 


3  .;  O 

r-  ta  o 
<i  r-  c  c 
i.  tv  2  c 

X  *-»  -4  O 
<  ^  . 

•  U 


o  ►—  o 
o  JO 
o  DC( 
a  o  x  o 
«-»  o  - 
o  u 


i  a 

*  d 

vo  o  a 

«.  O'  <  a 
i i-i  i  & 
X  .f>  «-*  d 
a  c\j  ^ 


oh  3 

0  3  0 
O  3  Oi 
3  0X0 

*-«  ©  • 

o  u  . 


h  iDC 
<r  ^  <  ol 
tu  tf  5 

a:  d  >-h  o 
<  to 

•  o 


0  1-3 
0  3  0 
D  3  o 
3  0X0 

#-h  o  * 

o  o 


•-<  O  c 
«r  «-(  <t  o 

u:  a  5  c 
a:  :*■  -•  o 

<  TO  • 

*  o 


O  JO 
*  <  ■ 

O  OJ  3  o 
X  a  c:  r~~ 
a.  O' 
bj  U  a 


<0  3  0 

•  a  • 

ONlilO 
X  o  c  o 
X  P"  <V 
cvi  o  m 


Figure  C.2.  (Continued). 
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Figure  C.2.  (Continued) 
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APPENDIX  D 


SEARCH  ROUTINE 


INTRODUCTION 


This  appendix  presents  a  description  of  the  search  procedure  employed  by  the  pro¬ 
gram  for  longitudinal  vibrators.  A  detailed  summary  of  naming  conventions,  equations,  and 
user  conventions  presently  existent  in  the  search  package  is  presented  as  well  as  a  sample 
run.  The  procedure  includes  a  decompaction  process  for  a  3-3  mode  ceramic  stack. 

The  first  step  of  this  search  procedure  is  the  optional  computation  of  a  velocity 
control  reactor  with  specified  properties  which  maximizes  the  velocity  control  of  the  trans¬ 
ducer  element  at  a  specified  frequency.  Following  the  loading  of  this  reactor,  the  trans¬ 
ducer  element  performance  is  computed  at  a  number  of  frequencies  over  a  band  centered 
about  a  specified  center  frequency.  1'q. 

At  each  frequency  within  the  frequency  sweep,  element  performance  is  computed 
for  each  of  a  number  of  radiation  boundary  conditions  by  means  of  a  boundary  condition 
loop.  Each  performance  variable  is  checked  and  only  the  maximum  and  minimum  of  each 
is  stored.  Several  computed  quantities  are  summed  over  the  boundary  condition  loop, 
providing  totals  and  averages.  Following  completion  of  the  boundary  condition  loop,  sever¬ 
al  quantities  are  computed  for  use  later  in  computing  two  new  quantities  as  a  function  of 
bandwidth  about  the  center  frequency.  Also,  selected  performance  variables  are  output  at 
this  point  by  printer,  and  either  card  punch  or  magnetic  tape. 

The  process  is  then  repeated  at  a  new  frequency  until  the  entire  frequency  sweep  is 
complete.  Following  completion  of  the  frequency  sweep,  the  two  “FBAND"  variables  are 
computed  as  a  function  of  bandwidth  and  output  by  specified  means. 
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MECHANICAL  NETWORK  CONVENTIONS 


All  longitudinal  resonator  transducers  using  3-3  mode  ceramic  may  be  compacted 
into  three  networks,  as  shown  in  figure  D.l .  The  double  arrow  (=>)  defines  the  direction  of 
positive  power  flow,  and  the  current  (velocity)  is  always  in  the  clockwise  direction.  The 
following  are  descriptions  of  the  basic  networks  and  their  ports: 

Network  A5S  is  defined  as  the  compacted  interconnection  of  networks  B4S,  C3S, 
and  T1S.  The  standard  interconnection  table  is 

1  S  1 A5S,  3C3S 

2  S  2A5S.3B4S 

3  S  3A5S.4B4S.  1T1S 

4  S  4A5S,  1B4S, -+1C3S 

5  S  5A5S,  2B4S,  -+2C3S  . 

Network  C3S  is  defined  as  the  compacted  (type  1  1  )  3-3  mode  ceramic  stack,  where, 

by  convention, 

1C3S  is  the  ceramic  mechanical  port  at  the  tail  end, 

2C3S  is  the  ceramic  mechanical  port  at  the  head  end, 
and  3C3S  is  the  ceramic  electrical  port. 

Network  B4S  is  defined  as  the  compacted  and  reduced  interconnection  of  the  re¬ 
mainder  of  the  transducer,  where,  by  convention, 

f> 


4A5S 


Figure  D.l .  Three  Basic  Networks  for  All  Longitudinal 
Resonator  Transducers. 
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1  B4S  is  the  mechanical  port  at  the  tail  end  of  C3S, 

2B4S  is  the  mechanical  port  at  the  head  end  of  C3S, 

3B4S  is  the  mechanical  port  at  the  radiating  face, 
and  4B4S  is  the  mechanical  port  at  the  tail  mount. 

Network  T 1 S  is  defined  as  the  compacted  and  reduced  interconnection  of  the  tail 
mount  assembly  and  is  of  the  form  of  a  complex  passive  tail  impedance,  Zj. 
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ELECTRICAL  NETWORK  CONVENTIONS 

The  electrical  network  naming  and  sign  conventions  assumed  by  the  procedure  are 
presented  in  Figure  D.2.  Also  shown  are  particular  voltages  (or  forces)  and  currents  (or 
velocities).  The  electrical  networks  are  described  as 

Network  S2S,  the  power  amplifier; 

Network  A2S,  the  transducer  element; 

Network  E2S.  the  velocity  control  network; 

Network  D2S  and  F2S,  the  matching  networks; 

Network  G2S,  [AG2S]  -  [AS2S]  *  [AD2S]  ; 

Network  H2S,  [AH2S]  =  [AF2S]  *  [AA2S]  ; 
and  Network  M2S,  [AM2S1  =  [ AD2S]  *  [AE2S]  *  [AF2S]  ; 
where  [AX]  represents  the  A-form  for  a  two-port  network  named  X.  and 


Figure  D.2.  Electrical  Networks  Assumed  in  the  Search  Procedure. 
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Similarly,  the  Z-form  t'or  a  two-port  network  is  defined  by 
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PORT  VARIABLE  DEFINITIONS 


The  following  symbol  definitions  and  relationships  for  port  variables  are  made 

1.  f£i’  v£i  =  force  and  velocity  at  the  left  (tail)  end  of  ceramic  piece  #i. 

2.  fri,  vri  =  force  and  velocity  at  the  right  (head)  end  of  ceramic  piece  #i. 

3.  Ec,  ici  =  voltage  and  current  at  the  electrical  port  of  ceramic  piece  #i. 

4.  FL,VL  =  force  and  velocity  at  Port  #1C3S  and  Port  #1B4S. 

5.  Fr.  VR  =  force  and  velocity  at  Port  #2C3S  and  Port  #2B4S. 

6.  Ec-  Ic  =  voltage  and  current  at  Port  #3C3S  and  Port  #1A5S. 

7.  Fj,  Vj  =  force  and  velocity  at  Port  #4B4S. 

8.  -Fj,  Vj  =  force  and  velocity  at  Port  #  1 T 1  S. 

9.  Erad’ ^rad  =  force  and  velocity  at  Port  #3B4S  and  Port  #2A5S. 

10.  Zrad  =  Fra(j/Vracj  =  radiation  impedance  connected  at  2A5S. 

11.  Zj  =  -Ft/Vt  =  tail  impedance  =  ZT1S;(Zy  =  0  for  air  termination). 

12.  Eg,  Ig  =  voltage  and  current  into  power  amplifier  at  1  S2S. 

13.  =  voltage  and  current  into  matched  transducer  at  1M2S. 


4. 

'Es“ 

d 

AS2S 

em 

Js_ 

5. 

_em 

d 

AM2S 

i 

m 

o 

_ i 

i 

3f 

1 _ 

_ 1 

Erad 
Vrad. 

17.  [AM2R]  £  [AM2S]  *  [AA2S] 

18.  [AS2C]  =  [AS2S]  *  [AM2S] 

19  d  [ AS2Sj  .[aM2Sj  .  [aA2Sj 


16. 


L‘cJ 


.  A2S 
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20. 


21 . 


23. 


24. 


f£i 

“V 

f  • 

lie. 

ZC 

vri 

ri 

ri 

E 

i  • 

c  _ 

_  Cl 

d  S2R,  S2R 

ec  Ai?  /A]  1 


7  d  S2R,  S2R 

Fec  =  1 /A^"^  *  (a  given  voltage) 
Fjc  =  1/A^J^  *  (a  given  current) 


— 1 

t- 

J 

vL“ 

fr 

d 

ZC3S 

VR 

1 — 

m 

o 

1 _ 

_lC  _ 

'fl  " 

_VL  ‘ 

fr 

d 

ZB4S 

VR 

Frad 

^rad 

i 

H 

U< 

_ 1 

vT 
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PERFORMANCE  VARIABLE  DEFINITIONS 


3. 

4. 

5. 

6. 

7. 

8. 
9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 


The  performance  variables  are  denoted  by  the  lollowing: 

T^  =  stress  in  the  \tjl  ceramic  piece  at  distance  (x)  from  the  left  (tail)  end  of  the 
piece. 

S*  =  strain  in  the  i ^  ceramic  piece  at  distance  (x)  from  the  lett  (tail)  end  ol  the 
33 

piece. 

P1  =  power  loss  in  the  \^}  ceramic  piece  (numbering  trom  lett  (tail)  to  right 
LC  (head)). 

77^  =  efficiency  of  the  irfl  ceramic  piece. 

E  =  electric  field  on  the  ceramic. 

PLC  =  total  power  loss  in  the  ceramic  stack. 
tjq  =  total  ceramic  stack  efficiency . 

Pld  -  total  power  loss  in  the  electromechanical  portion  oi  the  transducer. 

Pq  =  total  transducer  electromechanical  efficiency . 

Pet  =  total  power  loss  excluding  only  the  amplifier  (including  Pyj_). 

rtj  =  total  efficiency  excluding  only  the  amplifier  (including  Pjj_). 

PM  =  power  into  the  tuned  transducer  (amplifier  output). 

Pc  =  power  into  the  electromechanical  portion  of  the  transducer  (into  the 
ceramic). 

PjL  =  power  output  from  the  tail  into  the  tail  impedance. 

Prad  =  power  output  from  the  head  into  the  medium  (water). 

Zec  =  voltage-velocity  control  impedance  seen  looking  into  the  radiating  face  with 
port  1  S2S  short-circuited. 
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17. 


18. 

19. 

20. 


■>  l 


23. 

24. 

25. 


26. 


27. 

28. 

29. 


Zj .  =  current-velocity  control  impedance  seen  looking  into  the  radiating  lace  with 
port  1  S2S  open-circuited. 

Zq  -  impedance  seen  looking  into  the  ceramic  port  (Port  3C3S). 

Zj yj  =  impedance  seen  looking  into  the  matching  network  (Port  1  M 2 S ) . 

Yq  =  admittance  equivalent  to  Zq. 

Yjyj  =  admittance  equivalent  to  Z^. 

^Arad  =  tota*  array  output  power. 

Pma  =  total  array  input  power. 

77-p^  =  total  array  efficiency  excluding  only  the  amplifiers. 

X-  =  the  location  of  maximum  stress  magnitude  in  the  1  ceramic  piece  expressed 
as  a  distance  from  the  left  end. 

S  th 

X-  =  the  location  of  maximum  strain  magnitude  in  the  i  ceramic  piece  expressed 

as  a  distance  from  the  left  end. 

=  reactive  power  output  by  the  amplifier. 

^LM  =  Power  l°ss  in  ^e  electrical  matching  network.  M2S. 

PM  =  efficiency  of  network  M2S. 
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30.  |SMI  =  complex  power  magnitude  at  port  1M2S  (out  of  amplifier). 

SM  =  eM  X  lM*  =  Re  (EM  X  (EM  X  =  PM  +  jQM 

31.  ZM  =  reactor  inserted  to  extremize  l7ecl  or  l7ic|  in  midband  logic. 
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USER  INSTRUCTIONS 


1 .  Specify  interconnection  table  defining  the  BnS  network  which  will  be  reduced 
to  the  B4S  network. 

2.  Specify  interconnection  table  defining  the  TmS  network  which  will  be  reduced 
to  the  T1  S  network. 

3.  Specify  all  network  parameters  and  types,  including  those  for  the  Type  1  1  C3S 
network. 

4.  Specify  [  AF2S]  ,  [AD2S]  ,  [AS2S]  and  Zself. 

5.  Specify  the  midband  logic  parameters  to  be  used  in  computing  velocity  control 
network  matrix,  [AF2^]  (optional). 

6.  Specify  the  array  of  radiation  impedance  and  velocity  boundary  conditions  to 
be  searched  at  each  frequency. 

7.  Specify  frequency  range  and  frequency  increment  to  be  investigated. 

8.  Specify  print,  punch,  and  tape  output  options. 
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SEARCH  PROCEDURE 


The  following  paragraphs  provide  the  basic  equations  and  conventions  found  in  the 
search  procedure.  The  five  basic  operations  are  summarized  below. 

(Set  frequency  of  maximum  velocity  control) 

Midband  logic  procedure  (optional) 

(Enter  frequency  sweep) 

Compaction 

(Enter  boundary  condition  loop) 

Port  variable  computation 
Performance  variable  computation 

(Exit  boundary  condition  loop  -  output  results) 

(Exit  frequency  sweep) 

FBAND ZM 
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Midband  Logic  Procedure  (Optional) 

Single-reactor  midband  logic  computes  the  series  or  parallel  reactor  where 
=  lmag  (Z^f/Real  (Z^)  and  =  Real  (Z^*).  required  to  extremize  iycJ  or  b-J  at  a 

IT  ~>C 

specified  frequency.  If  this  subroutine  is  not  called.  [  A  ]  is  automatically  set  to  a  unit 
matrix  and  the  following  procedure  is  ignored. 

r  "ic  u  -ic 

1.  Compute  [A  -°]  and  [A  ]  from  available  information. 


N0  =  Kkl  ZSELF  +  Kk2 
N1  =  KC1  ZSELF  +  K£2 
D0  =  Kkl 
D]  =  Kgl 


D- 15 


3.  Compute  reactor  constants  (caution:  is  always  series). 

M0  -  Nq  +  RM  N  j 

Mj  =  jN,  +Nj/QM 
c0  =  d0  +  rmd, 

C,  =  jD,  +D1,/QM 

ao  =  |M0|2 

“1  =  -(M0rMlr  +  M0iMliJ 
a2  =  |M  j 

Po  =  ico|2 

h  =  ^c0rclr  +  c0lcH) 

=  1C]  I2 

S0  =  (al  %  "  a0  ^1 } 

6j  =  2  (ai  (3q  -  aQ  j32) 

6o  —  (a-i  |3]  -ckj  /3-> ) 

e0  =  ^0  61  "  2  01  50> 

el  =*--^0^2_^l^l~  ^@2  «0> 

€~>  —  (—3(ii  6  j ) 

e2  —  (-2(3')  5  -i ) 
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4. 


Compute  reactor  variables: 


=- 


25- 


-5  j  ±  ^6  j  -  -  45-i  5q  j 


1/2 


O  2 

en  +  e  i  X~  +  6t(X“)  +e->(X~) 
D±  =  - ! - = - - 


0o  +  0,  x-  +  p2  t'X-)" 


If  D1  <  0.  X1  maximizes. 
IfD1  >  0,  X111  minimizes. 


5.  Select  x  =  x±  to  maximize  (or  minimize). 
If  x  >  0  set  REACTOR  FLAG  =  L  , 
and  REACTOR  =  x/co. 

If  x  <  0  set  REACTOR  FLAG  -  C, 

and  REACTOR  =  -  —  . 

cox 


p  -)  c 

6.  Load  A  from  following  equations: 


.  E2S  .  E2S  .  ,.n 

Aj  j  =  Ani  =  1  +j  0 


If  REACTOR  FLAG  =  C,  then  ZM  =  RM  + 


j co  REACTOR  I  QM 


1  + 


If  REACTOR  FLAG  -  L,  then  ZM  =  RM +jco  REACTOR  (\  --J- 


If  parallel  reactor,  then 


aE2S  =  1/ZM 


and 


E°S 

A1  2  -  0 . 
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If  series  reactor,  then 


Compaction 


The  following  procedure  is  performed  once  at  each  new  frequency  within  the  fre¬ 
quency  sweep. 


1.  Compute 

[Z«S], 

[ZB4S] 

.  and  [ZT1S] 

2.  Compute 

[ZA5S]  . 

.  1  zA-S  j 

.  and  [AA“S 

DCJ 

3.  Compute  left  boundary  condition  matrix.  [K  ] .  where 


[K 


BCLi  3 


ASS 

'21 

?A5S 

Z23 

A>S 

7  1 

24 

7  A5S 
25 

ASS 

'31 

7  A5S 
33 

7  A5S 
34 

7  ASS 
Z35 

ASS 

7  ASS 

7ASS 

7  A5S 

'41 

Z43 

Z44 

Z45 

ASS 

7  ASS 

7A5S 

7  A5S 

'51 

Z53 

Z54 

zss 

4.  Compute  the  piece  matrix.  [Z^  ] ,  where 


=  -Z^  =  -jpc.A  cot  (kC ) 


Zp  =  -Z^  =  jpcAcsc(kC). 


,C 


,c 


c 


,c 


Zj3  =  Z, ,  =  Z^  =  -Z^  =  -j 


23 


(f) 


?33 


and 


7C  _  .  pcE  r  ,.D  I  " 
Z33  J  kA  Lg33  T  S33/e33_ 


Note:  These  equations  are  valid  only  for  type  1  1  ceramic. 
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r  r 

5.  Compute  right-handed  decompaction  matrix,  [K.1^] ,  where 


4C  =  zn/D 


Compute  [AS-SI .  [AD2S] ,  I AF^]  .  |AE2S,  .  and  |aM2Sj 


7.  Compute  Zec  and  Zjc  from  definitions.  Then, 
Tec  =  ZSELF  +  Zec  ’ 

Tic  =  ZSELF  +  Zic  * 

where  Zgppp  is  self  impedance. 
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Port  Variable  Computation 


The  following  procedure  is  used  in  computing  the  port  variables  for  every  network 
involved  in  the  interconnection  (decompaction)  and  is  performed  for  a  single  radiation 
boundary  condition  within  the  boundary  condition  loop. 

1 .  Compute  [K®^®]  ,  the  right  boundary  condition  vector,  where 


4.  Compute  Fj  =  -Z^*®Vy. 

5.  Compute  [£]  -  I  A***]  [J] 
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6.  Compute 


EM 


7.  For  each  ceramic  piece  in  the  stack,  compute 


V 

f£i 

vri 

-  [K^] 

v£i 

_  ^ci_ 

Er 

_  c  _ 

for  the  i^1 


ceramic  piece. 


where 


fri 

d 

fe.  i+i 

vri 

vC.i+l 

L  -I 

r 

ff«i 

I - 

< 

— i 

> 

_ i 
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Performance  Variable  Computation 

The  following  procedure  is  used  to  compute  the  performance  variables  from  a  single 
set  of  port  variables. 

All  performance  variables  are  searched  for  both  maxima  and  minima  over  the  speci¬ 
fied  boundary  conditions  at  each  frequency  except  for  the  cumulative  array  performance 
variables  of  Step  1  2. 

1.  Compute 


where  N  is  total  number  of  ceramic  pieces. 
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and 


7.  Search  over  all  ceramic  pieces  for  these  ceramic  variables. 


TM  =  MAX 


IT*  I  .i=1 
max 


■] 


=  MAX 


77^  =  MIN 


SM  =  MAX  IS*  |  ,  i=l ,  .  .  .  ,  N 

|_  max  J 

[<4  . N] 

-i=t . nJ 


i£?  =  MAX  [ici,  i=l . N]  . 


Compute  exterior  ceramic  stack  performance  variables. 

PLC  =  Re  (flvl  +  ec1c"fRvr) 

pc  =  Re  (ecI*) 

E  =  Eq/Q.  ;  £  =  ceramic  piece  length 

zc  =  EC/IC 
Yc  =  kite 
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9.  Compute  exterior  electromechanical  performance  variables. 
PTL  =  ~Re  (ZTVTVt) 

Prad  _  Re  ^Zrad  Vrad  Vrad^ 

1 0.  Compute  amplifier  boundary  condition  performance  variables. 
PM  _  Re  (  EM  ) 

Qm  ~  (em  ) 

/  -I  t\1/2 

|sm>  ~  (pm“  +  Qm  / 

ZM  =  EM^M 

ym  =  !m/£m 

1 .  Compute  transducer  performance  variables. 

PLT  =  Re  (EMlM)'Re  (ZradVradV*ad)  =  PM  ~  Prad 
t?x  =  100% .*(1  -pLT/pM> 

V/£s  =  Vrad/Es  Prad/lEsl2 

"  Yrad^s  Prad'/^s^ 
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Cumulative  array  performance  variables  are  these 


PMA  "  J  PM  ' 

PArad  "  zl  Prad  ■ 

r?TA  =  100%  *  ( P Arad/PM  A )  ■ 
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FBAND  ZM 


Following  completion  of  the  frequency  sweep,  the  ratio  of  ZM  MAX/MIN  is  com¬ 
puted  (magnitude  and  angle)  as  a  function  of  bandwidth  (in  Hz)  centered  about  a  specified 
center  lrequency.  The  following  procedure  is  used. 

1.  FBAND  |ZM| 

This  performance  parameter  is  defined  in  the  following  manner: 

Given:  fg  =  center  frequency  of  the  u,  •„  rating  band. 

f()  +  fj ,  ff)  _  f  i  -  edge  frequencies  of  band  BW  j . 

IZj^jl  maxj  =  max  mag  ZM  in  band  BW "j  . 

IZjyjIminj  =  min  mag  ZM  in  band  BW/ ] . 
d  |ZM|max, 

Thus,  FBAND  MAG  ZM  = -prn — — Lat  BW  =  BW,  =  2f, 

|ZM|minj  1  1 

Similarly  for  BW;  between  f0  +  fj  and  f0  -  fj,  FBAND  MAG  ZMj  is  computed, 
fj  is  stepped  from  zero  to  one-half  of  operating  band. 

FBAND  MAG  ZM  is  to  be  plotted  against  BW  and  not  FREQ. 

2.  FBANDZ.Zm 

This  performance  parameter  is  defined  in  similar  manner  to  the  foregoing: 

Given:  fg,  1q  +  f  j .  fg  -  fj ,  BW,  (ZLZj^j  )  max  j ,  (ZZ^j)  minj  as  before. 

Thus,  FBAND  ANG  ZM  =  (ZZM)maxj  -(ZZM)minj  at  BW  =  BWj. 

FBAND  ANG  ZM  is  plotted  against  BW  also;  but  the  units  are  “degrees”  in¬ 
stead  of  “unitless  ratio”  as  was  the  case  above. 
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Ceramic  Stress  And  Strain 


The  following  development  provides  a  detailed  documentation  of  the  methods  re¬ 
quired  in  D.6.5,  steps  2  and  3  to  compute  stress  and  strain  existent  at  the  interior  of  a  ce¬ 
ramic  piece.  The  symbols  used  here  are  defined  in  Appendix  B. 

1 .  The  stress  in  the  \th  ceramic  piece  at  a  distance  x  from  the  left  edge  of  the 
piece  is  given  by  this  expression: 

T^3  =  K  j  +  K'2  cos  kx  +  K3  sin  kx  , 

where 


K, 


.  £33  ‘ci 

J  A  D 

cc A 


K'2  =  -jpc  [VC  j+1  esc  kC  -  Vgi  cot  kfi]  , 
K3  =  jpc  VCj  , 


and  k  -  co/C  . 

Assuming  that  higher  order  terms  may  be  neglected  (i.e.,  small  angle  approxi¬ 
mation),  then 


T^  -  kJ  +  k[X  +  kTx2, 


where 

,T 


^0  =  K 1  +  K  2  = 


T  T 


K!  =  kK3  = 
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and 


and 


,T 


T  ,  :  „T 


Ki  =  -1/2  k-  K'-,  =  K  *  +  j  K,.  . 


dr  J  di 


Differentiating  IT^f ,  equating  to  zero,  and  solving  for  x  we  find  that 


x  =  xl  =  -i/: 


M 


T 


M 


T  ’ 


where 


T  T  - 
m‘  =  IKq I  , 


«•  =  2(KSXr  +  KJiKn)' 


T  T  - 
M*  =  IK  |  +2 


(-IAAA) 


The  strain  in  the  \^1  ceramic  piece  at  distance  x  from  the  left  edge  of  the  piece 
is  given  by 

S^3  =  (Kt  cos  kx  + K3  sin  kx) 

By  again  making  the  small  angle  approximation,  differentiating  IS33 1~ ,  and 
equating  to  zero,  one  finds 


M 


-1/2 

Mi 


where 
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and 


=  2  (k^  Kr  +  K^  K 


S  S  S  S 
n  K,  +  K^.  K,. 
Or  lr  Oi  h 


■  ) 


M 


2  lK?':  +  2(KorK:r  +  Ko,K:,) 


KS 

ko 


=  Ko  S 


D 

33 


S  S 

KOr+J4 


Kr  = 


,,T  D 

kl  S33 


KS  +}Kb 
lr  1  i 


S33  =  K^+jK^.. 


3.  Search  procedure  for  maximum  and  minimum  values  of  IT33I  and  IS3 3 1 .  If  x 
y  P 

or  falls  outside  of  [0,  £]  inclusive,  then  discard  that  particular  value.  Other¬ 
wise,  compute 


IT33'x=0=  K)I,J 
is«'x=o  =  (Mo)1/2 

lT33l  =  IK']  +  (K'o  cos  kC  +  K3  sin  kfi)| 


IS33I 


X=C 


(K'o  cos  kC  +  K3  sin  k£)| 


lT33!x=xT  =  IK1  +  K7  cos  ( kXj)  +  K3  sin  (k  )  I  if  0  <  xj  <  £ 
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SUMMARY 


To  summarize  the  requirements  of  the  search  procedure,  the  user  must  specify 

(1)  the  interconnection  table  that  defines  the  BnS  network,  which  will  be  reduced 
to  the  B4S  network, 

(2)  the  interconnection  table  that  defines  TmS  network,  which  will  be  reduced  to 
the  T1  S  network. 

(3)  [AF2S]  ,  [  AD“S]  ,  and  [AS~S]  .  which  can  be  done  either  through  interconnec¬ 
tion  tables  or  constants,  and 

(4)  the  frequency  range  and  frequency  increment  to  be  investigated. 

He  must  also  input  an  array  of  radiation  impedances. 

SAMPLE  SEARCH  ROUTINE 

Figure  D.3  shows  a  diagram  of  the  transducer  element  used  in  the  search  example:  a 
listing  of  the  control  program  and  associated  printed  output  is  shown  in  figure  D.4. 
Computer-generated  plots  from  the  search  procedure  are  shown  in  figure  D.5. 
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Figure  D.4. 

Listing  and  Output  of  the  Control  Routine  Used 
in  a  Search  Procedure 
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Figure  D.4. 


UjF  UR  #  I S  C  OLTPL .  COn  T  Rl 

f ok  ouio-io/?o/72-iu:l;>:-*5 


cm 


u. 

o 

< 

a 


«  >~ 

G  <r  < 

aJ  IX  i/)  U 

i  <wo 

o  uiri 

D  tr 


(\J  >C  OH 
O 

O  O  O  O 

o  o  o  o 

c  c  c  c 

o  o  ©  o 

C-hQ; 

fO  O  O  O 

O  C  O  Ol 

o  c  o  o 

G  O  O  O  ■ 


<x 

Q 


o  z  >- 

Z  Z  <t 
<x  i/it-  cr 
cr  cm  z  cr 

U.  o  o  < 
a  Z 


CN  G  IT) 

□  ONH 
C  G  C  C 

o  o  c  o 

G  G  O  O 
o  o  o  o 


O  G  O  O 
O  G  O  O 

3  3  a  G 
o  o  o  o 


h-  s  z  z  3-  2»!X  —I  O  ** 

<i  <i  •-»  >-  cu  fsij  Q  xm_i2q 


_j  *-«  cr  ^  cm  in  s 

hSa|aoQaogouaH<  JzsaLua 


COCQCZ-JO  ^  O  Ll1  Ctrl  cr  hUICCLH 
3  J  _|  J 


IT.  -£MC 
O  O  O  -H 
ooojoosaoo^ 
oooiooooooj 


cuKicririvOr'O-HWKjj 
•h  ^  oj  c\»  cj  oi 
o  o  o  o  o  ol 
o  o  o  o  o  o 


3-  in  o  -( 

cm  oj  cm  <\j  rO'fT 
3  0,0000 
O  O  O  Ol  o  o 


it  a  qh 
a  *-  *  i/r 

z  Z  z  3 


oj  fO  ?  ifl 
KJ  K»  rO  iO 
o  a  o  ol 
o  o  o  o 


ZL  Z  O  N 
u.  z  o  xju 
DO*-*  3Ejl^ 

O  OJ  o|k) 
O  o  CM  O 

o  o  c  o  c 

3  3  3  0  3 
O  O  O  O  O 

O  OO  o  o 
Q  Q  ~  ~  a: 

■ro  0  0  3-10 
O  O  O  O  IO 
3  3  3  013 

o  o  c  c  o 


z  z 

z  <1 

c  >  Z 
I—  Z  N!  U! 

<020; 
a  o*-h  _j 

3  3"  r-  o; 
o  o  <«h  o 

3  0  0  3' 
3  3  0  0 
O  O  O  O 
3  3  3  0 

a  a  h  h| 

ai  o;oh 
•o  oo  «-< 
o  o:c  o 
iC  c.o  o 


:lo  q.  —*  *— •] 

»0  3-  'sO  CM 
c\j  oo  cm 
o  O'O  o 
o  o  o  o 
C  O’O  o 
3  33  3 

a - 1 

o  to  o  o 

3  3  0  O 
O  OO  O 
0  3  0  0 


Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued) 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued) 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 
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Figure  D.4.  (Continued). 


D-69 


mom 

vD  GD  IT) 
CM  CQ  f- 


CVJ  "O 

K)  ro  m 


o  co  -£> 

CM  CM 


in  ir>  cm 


m  in  in 


in  in  i/> 
emu 


m  j  *o 


n  mm 


i 


in  in  i/i 


in:  d  rv 

K)  K)  rO 


COMMON  CATALb 


Lu 

O 

< 

CL 


CV 

r~ 

o 

CV 

o 


< 

Q 


S)  A  N  D:- t(OJ>h-J'-irO^)r~7'^d‘OCOOCV.n^ 

tvcvfOfOj-j^-j-if^nj^irirtni/'ivOvrl^cvor^r^r-r-i 


,  oj  m  ro  *  n  o 1  -  -  •  ~  ‘  —  -  • 

cr  cc  cc  cc  cc  cc 


O  3  3  3  »-4  Oi 


O  l-j  O  3  3  *-♦  • 


o  o  o  a  — *  ca 


la 


cO  -H  lD  (M  tv  cd 


i/>  u)  Si  <  <  trj 
33  H>  <  «i  <  jd 


-« <g  »o  sr  si  >q 


r-  o  r-*  o  r-  o 


l  O  <H  3  -H  3j 


HQHOHd 


cvi  tv  tv  cv  tv  cvb 


□j  uJ  lu  x  >  <aq 
n  lu  uj  :>  :>  H 


r-  cn  O'  o  h 


r-r-  ?  o  p-  coi  v  j  h  (\j  cr  , 
COffO'CC^t-'^WtL  f\J  i 


p-p-or-t-r^'or-  sn  3  io; 


*-H  *-4  O  •“<  W  *H 


AJ  <\l  lO  K>  -O  (O 


<\l  o  ^  fV  rq 


«  <1  C  cc  cr  £T 
h-  h;  oc  oq  x  a 

»o  ^  in  Oj  r-  33 


o  — •  ou  -<  J")  -h 


rO  rO  KJ  fO  O  O 


3  3  3  3  O  3.  O  O 


3  3  30  33  3  0  <\ 


0  0  3  0  0  31  30330  3;  3030  X 


3  0  300  0  3  0  3  3  0  3 


o^cvrOoo'oooocojoooooa 


OJ  CV  CV  rv  ro 


c\j  cv  cv  tv  <m  cv 


tn  u  o  u 
<  p  nj  as 


isiNifsjfM»-t/iLni/Ti/)t/tLnLr  s  Si  s  si  s>  s 
Nj.'jNirNirMciu.pL/iouJi:  2:  <a  x  p  <  <j 

O'  3  -t  tV!  fO  Si  D>flODW(\JfOj'All 
•-«  CVJ  CV  CV^AJ  C\J  OU  CV'  CV  C\J  V  ■OrOifOfO'OKIrO'O'O^ 


(V  &■  O  3  OJ  3 


0  0  3  3 


O  O  O  O  S) 


O  ^  lO  <0,  Ld 


r«-  ay  O'  o;  o 


Figure  D.4. 


(Continued) . 


D-70 


O' 

ir 


Figure  D.4.  (Continued) 


D-71 


Figure  D. 5 . 

Computer- generated  Plots  from  Search  Procedu 


hue  run  me  run  no  i  omt  10/20/72  h.u.c.  oissrw  version 


Figure  D.5.  (Continued). 


D-74 


31S0.0  3162.5  3175.0  3187.5  3500.0 

FREQUENCY  'Hz) 


unt  m*  Hi6  run  h o?  date  10/20/72  h.u.c.  disstlp  version 
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Figure  D.5.  (Continued). 
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Figure  D.5.  (Continued). 
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Figure  D.5.  (Continued). 
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APPENDIX  E 


FIND  PARAMETERS  EXAMPLE 


The  linear  approximation  iteration  technique  used  to  compute  the  set  of  ceramic 
parameters  is  described  in  this  appendix.  The  technique  designated  “find  parameters”  com¬ 
putes  the  ceramic  parameters  s|^3,  §333,  E^,  5^33  1,  §333  I,  eJ3  ij  from  the  experi¬ 
mentally  measured  set  of  quantities  [fm.  fn,  C  (f^X)).  lYml*.  lYnl,  and  D  -  where 

fm  =  the  frequency  of  maximum  admittance  magnitude, 

fn  =  the  frequency  of  minimum  admittance  magnitude, 

|Ym|  =  the  electrical  admittance  magnitude  at  frequency  fm. 


|Y n |  =  the  electrical  admittance  magnitude  at  frequency  fn, 

C  =  the  electrical  input  parallel  capacitance  at  a  frequency  f^Q  which  is  below  fm 
and  near  the  lower  frequency  limit  of  the  actual  transponder  operating 
frequency  band  (called  CLO  in  program). 


and  D  =  the  dissipation  factor  at  the  frequency  fj^Q  (called  CLOM  in  program ). 


*Thc  absolute  value  signs  will  be  omitted  from  Ym  and  Yn  henceforth  in  this  section  even  though  absolute  values  are 
used. 


E-l 


For  simplicity  of  explanation  of  the  method  employed,  the  subscripts  and  superscripts  on 
the  ceramic  parameters  will  be  omitted  in  the  remainder  of  this  presentation,  except  to 
designate  which  set  is  being  discussed.  That  is.  tire  set  of  ceramic  parameters  w  ill  be  |s.  g.  E. 
si.  gl.  El]  . 

The  equations  used  to  calculate  the  '‘measured"  quantities  from  the  parameters  are 
dependent  upon  the  exact  configuration  of  the  stack  under  investigation.  The  complete 
mathematical  model  of  the  stack  will  dictate  the  equations  to  be  used.  However,  in  general 
each  of  the  “measured"  quantities  is  a  function  of  every  material  parameter,  as  implied  by 
the  following  set  of  equations: 


fm  =  fm  is.  g.  E.  si.  gl,  El) 
In  =  fn  (s,  g.  E.  si.  gl.  El ) 


CHlo)  =  C  (s.  g,  E,  si,  gl.  El) 


imac  part  of  Yf 

fLO 


2x  7T  xfLO 


Ym  =  Ym  (s,  g.  E,  si.  gl.  El ) 
Yn  =  Yn  (s.  g.  E.  si,  gl.  El) 


D  (f 


LO 


=  D  ( s.  g.  E.  si.  gl,  El )  = 


real  part  of  Yp 


LO 


imae  part  of  Yf 

*LO 


) 


A  useful  relationship  can  be  set  up  between  these  quantities  by  taking  the  totai 
differentials,  where  f^  q  is  considered  to  be  a  constant  for  any  particular  case.  This  yields 
the  following  expression: 
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3  t  m 
oE 


"1 

dfm 

3fm 

3  s 

dfn 

3  fn 

3s 

dC 

3C 

3s 

dYm 

3Ym 

3s 

dYn 

3Yn 

3s 

dD 

dD 
_  3s 

3ln, 


o '  I  i 


3g 


3  !  m 


Although  equation  (2)  is  developed  tor  infinitesimals,  it  should  have  approximate  validity 


for  finite  differences.  One  estimates  a  set  of  material  parameters  [s.  g.  E.  si.  gl.  El]  .  From 

o  o  o  o  o  o 

this  initial  estimate  one  calculates  a  set  of  quantities,  [fm.  fn.  C.  Ym.  Yn,  D]  ,  from  the 

o  o  oo  o  o 


equations  of  the  complete  mathematical  muuel. 


A  set  of  increments  (As.  Ag,  AE,  Asl.  Agl,  AE1]  is  prescribed.  A  1  Off  change  in  s.  g, 
o  o  o  o  o  o 

E,  si.  gl,  and  El  is  used.  One  augments  each  of  the  initially  estimated  material  parameters 

one  at  a  time  leaving  the  others  alone.  A  new  set,  [fm,  fn.  C.  Ym,  Yn.  D]  ,  is  then  calculated 

for  each  of  the  six  augmentations.  For  example,  when  As  is  chosen  to  augment  s.  and  no 

o  o 

other  material  parameters  are  augmented,  one  obtains  the  following  set  of  values: 
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fm  (s  +  As,  g,  E,  si,  gl,  El)  , 
o  o  o  o  o  o  o 


fn  (s  +  As,  g,  E,  si,  gl,  El)  , 
o  o  o  o  o  o  o 


C  (s  +  As,  g,  E,  si,  gl,  El)  , 
o  o  o  o  o  o  o 

(3) 


Ym  (s  +  As,  g,  E,  si,  gl,  El)  , 
o  o  o  o  o  o  o 


Yn  (s  +  As,  g,  E,  si,  gl,  El )  , 
o  o  o  o  o  o  o 


and  D  (s  +  As,  g,  E,  si,  gl,  El) 
o  o  o  o  o  o  o 


If  we  subtract  the  quantities  computed  prior  to  augmentation,  expand  in  a  Taylor 
series  and  ignore  the  higher  terms,  we  have 


fm  (s  +  As,  g,  .  .  .)  -  fm 
o  o  o  o 


3fm 

3s 


As,  etc. 
o 


Thus,  the  following  linear  approximations  are  obtained: 

fm  (s  +  As,  g,  .  .  .)  -  fm  fm  (s,  .  .  .  El  +  AEI)  -  fm 

3fm  .  o  o  3fm  o  o  o 

3s  _  As  5  •  •  '  3EI  =  AEI 

o  o 


fn  (s  +  As.  g,  .  .  .)  -  fn 

3fn  ^  o  o  o _ 

3s  -  As 

o 


D  (s  +  As,  g,  .  .  .)  -  D 

3D  ^  o  o  o _ 

3s  -  As 

o 


3D 

3E1 


D  <s,  ...  El  +AED-D 
o  o  o 
AEI 
o 


(4) 
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By  this  procedure,  a  numerical  value  is  obtained  for  each  of  the  entries  of  the  partial  differ¬ 
ential  matrix  described  above  in  equation  (1).  After  this  matrix  is  obtained,  its  inverse 
matrix.  A,  is  found.  Now  one  can  write  the  finite  differences  of  the  set  [s,  g,  E,  si,  gl,  El] 
in  terms  of  the  finite  differences  of  the  set  [fm,  fn,  C,  Ym,  Yn,  D]  .  This  equation  has  the 
form 


As 

A11  A1 2  A1 3  a14  A1 5  A]  6 

Afm 

Ag 

A21  A22  A23  A24  a">5 

Afn 

AE 

A3  1 

AC 

ASI 

A4 1 

AYm 

Agl 

A5  1 

AYn 

AEI 

_A61  A66_ 

AD 

(5) 


One  now  chooses  Afm  to  be  the  difference  between  the  experimentally  measured  value, 

fmex  and  the  value  fm  calculated  from  the  estimated  ceramic  parameters,  and  likewise  for 
o 

the  other  quantities.  Thus 


Afm  =  fmex  -  fm  , 
o 

Afn  =  fnex  -  fn  , 
o 

AC  =  Cex  -  C, 
o 

AYm  =  Ymex  -  Ym, 
o 

AYn  =  Yn£x  -  Yn  , 
o 

and  AD  =  Dex  -  D  . 

o 


(6) 
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After  the  matrix  multiplication  indicated  in  equation  (5)  is  completed,  the  result  mu  changes 
in  the  material  parameters  [As.  Ag,  AE.  Asl.  Agl,  AEI]  are  added  to  the  initial  climates  [s. 


g,  E.  SI,  gE  El]  to  obtain  the  first  corrected  set  of  material  parameters  [s,  g.  E.  si.  gl.  El] . 

If  the  equations  were  nothing  more  than  a  linear  set  of  equations,  this  corrected  -ot  ot  mate¬ 


rial  parameters  would  be  the  exact  solution  which  goes  with  the  experimentally  measured 

quantities.  But  since  this  is  not  the  case,  this  corrected  set  of  material  parameters  [  s.  g.  E, 

si.  gl,  El]  is  substituted  into  equation  ( 1  )  to  obtain  a  second  calculated  set  [fm.  In. 

Yn,  D]  .  To  insure  a  refinement,  a  new  set  of  increments  [As.  Ag.  AE,  Asl.  Agl.  AEI]  is  also 

j|  lit*  1  i 

formed,  each  member  of  which  is  formed  from  the  old  increments  by  multiplying  such 
increments  by  0.1  :  As  =  As  .1 .  etc.  Thus,  starting  with  the  procedure  described  ndow 
equation  (2),  this  iteration  may  be  repeated  as  many  times  as  desired.  The  iteratn  n  is  ternti 
nated  if  each  of  the  following  terms  is  less  than  the  corresponding  member  of  some  preset 
list  of  increments. 


fmex  ~  fm 
fmex 


<  5fm 


fnex  ~  fn 
^nex 


<  6fn 


C  -  C 
'“ex 


C 


ex 


<  SC 


Ymex  ~  Ym 
Ymex 


<  SYm 


(7) 


Ynex~Yn 
^  nex 


<  SYn 


Dex-D 


D 


ex 


<  SD 
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Tins  test  is  made  just  following  the  calculation  of  each  corrected  set  of  material  parameters. 
With  5fm  =  Sfn  =  SC  =  .3*  10“^  and  6Ym  =  5Yn  =  SD  =  10"^.  not  more  than  two  iterations 
are  in  general  necessary  if  one  uses  reasonably  good  initial  estimates  of  the  material 
parameters. 

A  word  of  caution  is  appropriate  here.  The  program  can  converge  to  impossible 
parameters  if  there  are  errors  in  the  transducer  input  data  or  experimental  data.  Therefore, 
close  inspection  of  the  results  should  be  made  to  insure  validity. 

Figure  E.l  shows  a  diagram  of  the  transducer  used  in  the  find  parameters  example, 
and  the  associated  interconnection  schematic.  Figure  E.2  shows  a  listing  of  the  required 
control  routine  and  printed  output. 
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Figure  E.  2 . 

Listing  and  Output  of  the  Control 
Routine  Used  in  Find  Parameters. 
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Figure  E . 2 .  ( Continued ) . 
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Figure  E.2.  (Continued) 


ru 


U1 

<£ 

< 

a 


C\J 

h- 

o 

OvJ 

o 


Z  Z 
o  o 
o  u 


fO 


o 

•M 

fS 

-a 

o 

o 


z  \ 

O  o 

o  rvi 


O  at 

u.  o 


O  3 
3  3 
D  O 
O  3 
O  O 


o 

T» 

z 


z 

z 

o 


X  t~ 
L>  o 

»-*'  f- 
LA  < 
<•  K 

n  z> 


PTJ  J- 
o  o 
a  o 


I 


I 

i 


SC< 

(_>: 

o 

J: 

» 


ui 

o 

z: 

Ub 

tr: 

Ui 

li.; 

LlJ 

cr; 


C 

z> 

U> 

; 

Xi 

U4 


o  a: 

c  *- 

q_  _n 


JT  -«  CVJ  -* 

O  b~!  ^  *■”!  UJ  QJ 

cl  o  o  o  z  g 

33S9r^ 


_J  o  kp 

m  ja  w^q. 

hhKI-MO 
Q.UZKOI- 
a  cm  qimit 

a.  cl  z  z  z  z 


iftNDh 
o  o  c 
o  o  a 

DOC 


O^OJKl^U' 

o  a  o  o  o  c 

o  o  o  O  O  C 


XNOH 
OJ  rvJ 


OJ  IT 
od  e\ 
o  a 
o  a 


-O  o 
o  o 

,ga 

■si 


O  J- 

■  o  a 
;  o  o 
.  o  o 


G_.< 

d\<J> 


53. 

ujIo  b 


t-  *  3 
<  U; 

O  qK  - 

Z  Z  CD  I 
O  O  ; 

LA  lA!.m 

-  m.:j>  ■ 

o  o[  X  I 

Ld  UJ 

X  a:, 

D-  £L  K 


05  -3Dj 
D  ID 
O  O 


s  \ 

•4 

Z  Z  »- 

■  O  O  <L 

:  z.  z  o 
z  z 

I  O  O  X 

I  J  J  £ 


■o  X  X)  43, f-  T3 


3  h- 

X  ID 

o 

cc  *: 

UJ  o 
U  UJ 
D  X 
n  o 
z  <r 

*-*  UJ  X 

_J  tA  t— 

uj 

O  or  -  • 

<  »-  Lf)  — 

CL  LA  t- 

K  -1 

_»  -J  Z  Z 

_ J  I  • — *  LTC 

<  <nr  o 
J  o  x  u 


<3  0  0 

:C0  x>  co 
Id  o  cd 

VI  (\|lO 


D  O  O 

■CL.  X  X 


1  -I  -I  -I 
_)  — 1  — l 
I  «r  c  <t 


7>  3,  h  M  A  tl 


ru  rO  iD  in  vO  v£>  r-o 


Figure  E.2.  (Continued). 
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Figure  E.2.  (Continued). 
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Figure  E.2.  (Continued). 
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Figure  E.2.  (Continued). 
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Figure  E.2.  (Continued). 
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Figure  E.2.  (Continued). 
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Figure  E.2.  (Continued). 
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APPENDIX  F 


ACTIVE  BAND-PASS  FILTER  ANALYSIS  EXAMPLE 


An  active  band-pass  filter  analysis  is  obtained  by  cascading  an  eighth-order  Chebyshev 
low-pass  filter  (figures  F.l  and  F.5)  with  an  eighth-order  Chebyshev  high-pass  filter  (figures 
F.2  and  F.5).  The  high-pass  and  low-pass  filters  are  each  composed  of  four  individual  second- 
order  active  stages  (figure  F.3)  in  cascade.  Each  stage  is  modeled  using  Type  26  Operational 
Amplifier  (figure  F.4)  in  conjunction  with  Type  21  Series  Impedance  and  Type  22  Parallel 
Impedance  (figure  F.5). 

The  voltage  gain  of  the  composite  band-pass  filter  is  computed  as  a  function  of  fre¬ 
quency  and  is  plotted  in  dB  (referenced  to  unity  gain)  along  with  the  measured  response  of 
the  prototype  circuit  being  modeled  (figure  F.6). 


F-l 


Figure  F .  1 .  Low  Pass  Filter. 
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Figure  F.2.  High-Pass  Filter. 
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Figure  F.3.  Individual  Second-Order  Active  Filter  Stage  Showing 

Sign  Convention  Diagram  and  Corresponding  Computer  Input. 
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Figure  F.4.  Type  26:  Linear  Model  Of  a  Differential  Input  Operational 
Amplifier  With  Active  Filter  Input  Values. 
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Figure  F.5.  Series  And  Parallel  Impedance  Input  For  Individual  Active  Filter  Stages. 
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Figure  F.6.  Plot  of  Voltage  Gain  as  a  Function  of  Frequency 


Figure  F. 7 . 


Listing  and  Output  of  the  Control  Routine 
for  an  Active  Band-Pass  Filter  Analysis. 
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Figure  F.7. 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F. 7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued) 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued) 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F. 7 . 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F.7.  (Continued). 
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Figure  F. 7.  (Continued). 
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Figure  F.7.  (Continued). 
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